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It is well known that the process of ion implantation revolutionized the semiconduc- 
tor industry. In case of metals, ion implantation has been proved to be very effective 
as a surface treatment technique which reduces wear and friction of the implanted 
surfaces and improves corrosion resistance and adhesion properties. These properties 
are the consequences of the stable/metastable/amorphous alloys produced in the near 
surface region due to ion implantation. However, the major disadvantage of ion im- 
plantation in case of metals has been the long irradiation times required to achieve 
sufficiently high implant concentrations resulting in high cost of the process. Further, 
the achievable concentrations are limited by the sputtering of surface atoms. Therefore, 
the process of ion implantation had a limited success in case of metals as compared to 
semiconductors. 

During the past decade new techniques have been developed which combine thin film 
deposition and heavy ion irradiation processes to synthesize novel and technologically 
important phases or materials, thereby reducing the process time and cost by an order 
of magnitude. These techniques utilize the inherent capability of heavy ion irradiation 
process of displacing the atoms from their normal sites either during thin film deposition 
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(ion-beam assisted deposition) or after irradiation ( ion-beam mixing). 

In case of ion-beajn mixing (IM), the elements to be alloyed are first sequentially 
deposited on the metal or semiconductor or inert surface and subsequently irradiated 
by energetic heavy ions. Ions of inert heavy gases are generally utilized in order to 
avoid any chemical reaction with the alloy to be formed. The sample configurations 
which are generally used for IM studies are i) marker ii) bilayer and iii) multilayer. 
In marker configuration one employs a thin layer ( ~ 1 nm) of an element in a host 
matrix to determine the effective diffusion coefficient of impurity resulting from the 
collision-induced atomic migration. In another configuration a marker layer placed 
across the interface is used to monitor the atomic movement during the formation of 
phase. The marker experiments mainly address the basic aspects of ion mixing process. 
The bilayer configuration involves a film of element A of thickness ranging from a few 
tens to a few hundreds of nanometer deposited on a bulk material B or on an inert 
substrate coated with a film of element B. The interest here is focused around the 
formation of binary phases of different compositions in addition to the understanding 
of the basic processes involved at the interface. Multilayer configurations are mainly 
important from the point of view of formation of stable/metastable/amorphous alloys 
of desired composition having larger thicknos.sc.s. 

For ion-beam mixing studies, inert gas ions having energy of few hundred kilo 
electron-volts are mostly used. This is partly due to the wide availability of the ion 
implanters at these energies. However, these ion energies limit the maximum thickness 
of the alloyed region to about 100 nm. The use of ion-beams of MeV energies is now 
gradually increasing for synthesizing thicker surface alloys. 

In order to make ion-beam mixing a viable technique for producing desired phases, 
an important prerequisite is to decipher the latent facet of the physics involved in 
the basic mechanisms under different processing conditions. Several ion-beam mixing 
studies have been reported for metal/metal and metal/Si systems which reveal the 
character of the mixing mechanism under various combinations of the processing pa- 
rameters. The details of these studies on phase formation and mechanisms of mixing 
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have been the subject of some excellent review articles by Cheng[l], Was[2], Kelly and 
Miotello[3], Nastasi and Mayer[4]. 

It is well known that metal silicides are important as contact and interconnect 
materials in semiconductor device technology. The interest in the synthesis of silicides 
by different methods grew in 1970s when it became clear that with decreasing dimension 
the resistance of the conductors fabricated from polycrystalline silicon would become 
too high for effective use. The work done in the last decade on metal silicide thin films 
revealed a large number of interesting properties related to resistivity, electromigration, 
diffusivity, corrosion resistance.etc. All these properties are important from the point 
of view of semiconductor and microelectronis industry which require such materials in 
modern integrated circuits because the packing density, speed, etc. are usually limited 
by interconnects and contact requirements rather than by decreasing active device 
dimension. For example, if the circuit size is scaled by a factor k, the interconnect 
resistance increases by a factor 1/k while the contact resistance by 1/k^. For the 
formation of stable and metastable alloys of metal-semiconductor systems ion-beam 
mixing has been widely used for its several advantages over solid state reaction such as 
lower processing temperature and higher spatial selectivity. In the last decade a large 
number of studies on the formation of different near-nobel and refractory silicides have 
been reported. 

In contrast to the silicide systems, which have been widely studied using both 
thermal and ion-beam processes, the reports on the synthesis of germanides are very 
few in the literature. Recently the interest in the germanide systems has been triggered 
by the work of Krusin-Elbaum and Aboelfotoh[5] who report that the ei- C uzGe phase 
possesses lower resistivity (few pO — cm) than many silicides {CoSi 2 i NiSi 2 ,Gtc.). 
Also it possesses most of the desirable properties required for an interconnect /contact 
material in the existing GaAs and Si device technology. This copper germanide phase 
which is formed by thermal annealing at about 680 K, is found to be stable and 
reproducible. The studies on Cu/Ge system indicates that other germanides which are 
mostly unexplored might emerge as promising future materials for the semiconductor 
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industry. 

The interesting and technologically important properties of Cu^Ge, mentioned 
above, motivated the author to undertake the study of MeV ion-beam mixing in 
metal/Ge systems. Apart from Cu/Ge, two other systems viz. Ni/Ge and Co/Ge 
have also been selected for these studies. These systems have been chosen for the fol- 
lowing reasons. Firstly, the atomic masses of Cu, Ni and Co being similar it is expected 
that the ballistic effects produced under ion irradiation in all these three cases would 
be similar. Also, the equilibrium phase diagrams of Cu/Ge, Ni/Ge and Co/Ge are 
similar. The main objectives of this work are (i) to explore the directed energy depo- 
sition process of ion mixing for producing germanide phases in bilayer and miiltilayer 
thin film structures of Cu/Ge, Ni/Ge and Co/Ge systems at different ambient temper- 
atures, using Kr and Ar ions of mega electron-volt energy, (ii) to obtain quantitative; 
information of the parameters viz. mixing rates, mixing efficiencies which are used to 
understand the basic mixing processes and the activation energies from temperature 
dependence of mixing (the Q curve) and (iii) to develop a phenomenological model to 
explain the mixing parameters. 

The thesis is divided into seven chapters. The importance of metal/Ge systems and 
the application of ion-beam mixing techniques are brought out in Chapter 1. 

A brief discussion of mechanisms of ion-beam mixing and various phenomenological 
models for calculating mixing parameters viz. mixing rate, mixing efRciencies,etc. are 
presented in Chapter 2. A compilation of experimental ion mixing results relevant to 
this work is also presented in this chapter. 

Chapter 3 deals with the methodology of sample preparation along with a descrip- 
tion of experimental techniques utilized in this work. High purity (99.99%) Ge, Cu, 
Ni, Co materials were used for preparing multilayer and bilayer films of M/Ge config- 
uration on quartz substrates where M is one of the above mentioned metallic species. 
•For marker experiments high purity Au was utilized in M/Au/Ge configuration. The 

multi-hearth electron gun evaporation system which was used for the deposition of thin 
films has been described. 
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For ion irradiations, the 2 MeV Van de Graaff accelerator at IIT, Kanpur was 
utilized. For naixing experiments, an ion irradiation chamber to work in the range of 77 
K to 500 K was designed and fabricated along with other ancillaries. Thermal annealing 
was carried out at various temperatures in the range of 400 - 700 K in vacuum. A brief 
description of these facilities has been described. Most of the quantitative information 
viz. layer thickness, composition of the mixed regions and depth profiling were obtained 
by Rutherford backseat tering spectrometry (RBS)^sing 1.1-1. 5 MeV He"*" ions from 
the same accelerator. The method of backscattering analysis has been explained by 
giving an example of Cu/Ge system. 

The experimental results and analysis on multilayer, marker and bilayer configu- 
rations are presented in Chapters 4, 5, and 6 respectively. These chapters deal with 
the production of specific germanide phases, movement of the atomic species during 
mixing and the mixing mechanisms respectively. 

Chapter 4 describes results of MeV Kr and Ar ion induced mixing effects [4], leading 
to compound formation in multilayer structures of Cu/Ge, Ni/Ge and Co/Ge samples. 
The ion mixing was done at various substrate temperatures ranging from 100 to 450 K. 
Irradiations were performed for several doses in the range of (2 — 20) x 10^® ions j cm?. 
In all the cases films of the metal M and Ge were sequentially deposited to produce 
multilayer structure of 8 to 10 layers having total thickness in the range of 90 - 300 
nm. The thickness of the individual layers of multilayer samples were adjusted such 
that the average film compositions range from about 60 at.% to 80 at.% of metal M. 

The results of ion-beam mixing experiments have been compared with those of 
thermal annealing experiments performed at various temperatures in the range of 400 
- 700 K. The characterization of all treated samples were done by RBS and X-ray 
diffraction measurements. The surface morphologies of the samples were studied by 
scanning electron microscope before and after mixing. The resistivity measurements 
were done by Van der Pauw four-probe method. 

It has been observed that MeV heavy ion irradiation is a very efficient low temper- 
ature process for producing stable phases in metal/Ge systems namely Cu/Ge, Ni/Ge 
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and Co/Ge at different temperatures. Cu^Ge, Ni^Ge and Co 2 Ge have been observed as 
“first phase” in ion mixed as well as in thermally annealed multilayer Cu/Ge, Ni/Ge 
and Co/Ge systems respectively. In case of Cu/Ge system another Cu-rich phase 
(Cu^Ge) has also been observed only at higher irradiation doses which is not observed 
in thermally annealed samples. 

In case of bilayer samples, the first phase produced by thermal annealing are Cu^Ge, 
Ni 2 Ge and CoaG'e in Cu/Ge, Ni/Ge and Co/Ge systems respectively. At higher an- 
nealing temperatures NiGe and CoGe phases have been observed as second phase in 
Ni/Ge and Co/Ge systems. The present studies show that the first phases formed by 
thermal annealing and ion-beam mixing are identical. 

The formation of compound phases has been analyzed on the basis of various ex- 
isting empirical rules. The composition of the first and subsequent phases have been 
explained in the light of the effective heat of formation rule which considers both kinetic 
and thermodynamic aspects. 

Chapter 5 deals with marker experiments which give information on the dominantly 
moving atomic species during ion mixing and thermal annealing. The bilayer struc- 
tures of metal M and Ge deposited on quartz substrate along with a thin marker layer 
embedded at the M/Ge interface were used for these studies. The thickness of indi- 
vidual layers was chosen such that it yielded the composition of about MsqGc^q. The 
total thickness of bilayer films was in the range of 100-200 nm. These experimental 
results are discussed on the basis of the present knowledge of the atomic movements 
in metal/semi conductor systems. 

It has been observed that during the formation of Cu^Ge phase by ion-beam mixing 
at room temperature, both Cu and Ge are mobile but Cu is the dominant moving 
spedes aoross the interface while for Ni/Ge system both Ni and Ge are equally mobile 
during the formation JVijffe phase. However, tor Co/Ge system only Ge atoms are 
mobile to produce Co2Ge phase. 

In Chapter 6 the values of different mixing parameters viz. mixing rates, mixing 
efficiencies, critical temperature for initiating thermally activated process etc. obtained 
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from bilayer experiments are described. The bilayer configuration was used since accu- 
rate quantitative information on these parameters can be obtained in this configuration 
only. Ion mixing of bilayer samples also leads to the formation of the phases described 
in Chapter 4 for multilayer samples. The room-temperature mixing rates for Ar ion 
irradiation in case of Cu/Ge, Ni/Ge and Co/Ge systems are found to be 5.4 nm'*, 5.6 
nw^ and 1.6 respectively. Similarly for Kr ion irradiations the mixing rates are 
20.6 nm^, 18.5 and 7.7 respectively. The mixing efficiencies are independent 
of the deposited energy density in each case and the experimental values for Cu/Ge, 
Ni/Ge and Co/Ge are obtained to be 13.1 nm^fkeV, 10 nm^fkeV and 5 nm^/keV 
respectively. 

The temperature dependence of ion mixing in all the three systems show that the 
critical temperature To, above which the mixing becomes temperature dependent, is 
just above room temperature. The value of Tc is about 320 K for both Cu/Ge and 
Ni/Ge while for Co/Ge system it is about 347 K. There is no distinction between 
critical and equivalent temperature for crystallization in these systems. Compound 
formation in all the cases occur around the critical temperature Tc. Empirical models 
which estimate the critical temperatures have also been discussed. 

In this chapter we have also presented a new phenomenological model[7] based on 
chemical kinetics applicable to reactive systems (such as metal silicides or germanides) 
to explain the experimentally obtained mixing rates since the existing models underes- 
timate the mixing rates for reactive systems. In this model the mixing rate is enhanced 
by the compound formation factor retaining the linear dependence on the deposited 
energy density. The predictions of this model are very close to the experimental re- 
sults presented in this work. It is also able to explain observed mixing rate of Ni/Si 
system reported in the literature. Hence it is demonstrated that this model would be 
applicable to the important class of silicide and germanide systems. 

The last chapter. Chapter 7 presents a summary of the work followed by important 
conclusions and scope of the future studies. 

To conclude, the synthesis and characterization of technologically important ger- 
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manicle phases by ion-beam mixing and tliermal annealing along with the experimen- 
tally determined mixing parameters (to understand mechanisms of mixing) and their 
successful interpretation with the help of a new phenomenological model, are the main 
outcome of this work which accomplishes the objectives mentioned above. 
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Chapter 1 
Introduction 


1.1 Preamble 

Prior to 1960 the research interests in ion-matter interaction were confined to the 
radiation damage studies of the surfaces of reactor materials exposed to neutrons and 
energetic fission fragments. The understanding of the impact phenomena improved 
when Bohr's theory of energy loss of energetic ions in solids was modified by Lindhard, 
Scharff and Schiott ( LSS)[1] in early sixties. The LSS theory presented the energ}' loss 
of an ion in a solid by universal stopping curves. These curves helped in delineating the 
energy regions useful for the purpose of materials modification and materials analysis 
as illustrated in Fig 1.1. The revolution in the field of ion-matter interaction occurred 
when it was realized that ion beams could be used for doping of semiconductors and 
could modify the tribological properties of metals{2] and that the light ions viz. ife+, 
etc. could be used for materials analysis[3]. Subsequently, a dramatic growth in 
the use of energetic ion beams took place for studying materials science and solid state 
physics problems[4-8]. In fact, ion beams of almost all stable elements of energies 
ranging from a few keV to a few MeV have been put to use. 
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Figure 1.1: Universal stopping power vs. ion energy in solid in Thomson-Perini (T.F.) 
unit as predicted by LSS theory. Sn and Se are the nuclear and electronic stopping 
powers respectively. 
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It is well known that the process of ion implantation revolutionized the semicon- 
ductor industry [4. 9]. In case of metals, ion implantation has been proved to be very 
effective as a surface treatment technique which reduces wear and friction of the im- 
planted surfaces and improves corrosion resistance and adhesion properties[10, 11]. 
These properties are the consequences of the stable/metastable/amorphous alloys pro- 
duced on the surface due to ion implantation. However, the major disadvantage of ion 
implantation in case of metals has been the long irradiation times required to achieve 
sufficiently high implant concentrations resulting in high process cost [6]. Further, the 
achievable concentrations are limited by the sputtering of surface atoms. Therefore, 
the process of ion implantation had a limited success in case of metals as compared to 
semiconductors. 

During the past decade new techniques have been developed which combine thin film 
deposition and heavy ion irradiation processes to synthesize novel and technologically 
important phases or materials, thereby reducing the process time and cost by an order 
of magnitude[6. 12-17]. These techniques utilize the inherent capability of heavy ion 
irradiation process of displacing the atoms from their normal sites either during thin 
film deposition (ion-beam assisted deposition) or after irradiation ( ion-beam mixing). 

1.2 Ion-Beam Mixing 

In case of ion-beam mixing {IM)[17], the elements to be alloyed are first sequentially 
deposited on the surface of a material and subsequently irradiated by energetic heavy 
ions. Ions of inert heavy gases axe generally utilized in order to avoid any chemical 
reactions with the alloy to be formed. The sample configurations which are generally 
used for ESI studies are shown in Fig 1.2. Multilayer configurations shown in Fig 1.2a are 
mainly important from the point of view of formation of stable/metastable/amorphous 
alloys of desired composition. The bilayer configuration involves a film of element 
A of thickness ranging from a few tens to few hundred nanometer deposited on the 
bulk material B (Fig 1.2b) or on an inert substrate coated with a film of element 
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B (Fig 1.2c). The interest here is focused around the forrnatioii of binary phases of 
different compositions in addition to the understanding of the ba.sic processes involved 
at the interface. In marker configuration (Fig 1.2d) one uses a thin layer ( ~ 1 nm) 
of an element in a host matrix. This configuration is used to determine the effective 
diffusion coefficient of impurity resulting from the collision induced atomic migration. 
In another configuration (Fig 1.2e), a marker layer placed across the interface is used 
to monitor the atomic movement during the formation of phases. These studies mainly 
address the basic aspects of ion-beam mixing process. 

For ion-beam mixing studies, inert gas ions having energy of few hundred kilo 
electron-volt are mostly used[ll]. This is partly due to the wide availability of ion 
implanters at these energies. However, these ion energies limit the maximum thickness 
of the alloyed region to about 100 nm. The use of ion- beams of MeV energies is now 
gradually increasing for synthesizing thicker surface aIloys[llj. 

In order to make ion-beam mixing a viable technique for production of different 
phases, an important prerequisite is to decipher many a latent facet of the phy.sics in- 
volved in the basic mechanisms under different processing conditions. Several ion- beam 
mi.xing studies have been reported for metal/metal{ 17-26], met al/Si [27-31] systems 
which reveal the character of the mixing mechanism under various combinations of 
the processing parameters. The details of these studies on phase formation and mech- 
anisms of mixing have been the subject of some excellent review articles by Was[6]. 
Cheng[14], Nastasi and Mayer[15], Kelly and Miotello[26]. 


1.3 Synthesis of Silicides 

It is well known that metal silicides are important as contact and interconnect materials 
in semiconductor device technology. The interest in the synthesis of silicides by different 
methods grew in 70s[32] when it became clear that with decreasing dimension the 
resistance of the conductors fabricated from poh’crystalline silicon would become too 
high for effective use. The work done in the last decade on metal silicide-s thin films 
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Figure 1.2: Schematic diagram of sample configurations utilized for ion-beam mixing 
and. thermal annealing studies. 
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revealed a large number of interesting properties related to resistivity, electromigration, 
diffusivity, corrosion resistance etc. [33, 34]. All these properties are important from 
the point of view of semiconductor and microelectronis industry which require such 
materials in modern integrated circuits because the packing density, speed etc. are 
usually limited by interconnects and contact requirements rather than by decreasing 
active device dimension[35, 36]. For example, if the circuit size is scaled by a factor 
k. the interconnect resistance increases by a factor 1/k while the contact resistance by 
1/A’^[36]. For the formation of stable and metastable alloys of metal-semiconductor 
systems ion-beam mixing has been widely used for its several advantages over solid 
state reaction such as lower processing temperature and higher spatial selectivity[5, 
6. 13, 15]. In the last decade a large number of studies on the formation of different 
near-nobel and refractory silicides have been reported[15, 17, 29. 30, 37-39]. 


1.4 Synthesis of Germanides 

.As compared to the silicides systems, which have been widely studied using both ther- 
mal and ion-beam processes[37, 40-43], The reports on the synthesis of germanides are 
very few in the literature. Recently, the germanide systems have attracted attention of 
many researchers [44-48]. This interest in the germanide systems has been triggered by 
the work of Krusen-Elbaum and Aboelfotoh[49] who report that the €i- CuzGe phase 
possesses lower resistivity (few /iO-cm) than many silicides {CoSh, NiSi-i etc.). Also 
it possesses most of the desirable properties required for an interconnect/contact ma- 
terial in the existing GaAs and Si device technology [44, 49], This copper germanide 
phase which is formed by thermal annealing at about 680 K, is found to be stable and 
reproducible. The studies on Cu/ Ge system indicates that other germanides which are 

mostly unexplored might emerge as promising future materials for the semiconductor 
industry. 
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1.5 Motivation and objectives of this work 

The interesting and technologically important properties[44, 49] of CuzGt. mentioned 
above, motivated the author to undertake the study of MeV ion-beam mixing in 
metal/Ge systems. Apart from Cu/Ge, two other systems viz., Ni/Ge and Co/Ge 
have also been selected for these studies. These systems have been chosen for the fol- 
lowing reasons. The atomic masses of Cu, Ni and Co being similar, it is expected that 
the ballistic effects[14] produced under ion irradiation in all these three cases would be 
almost same. Also, the equilibrium phase diagrams of Cu/Ge, Ni/Ge and Co/Ge axe 
similar. The main objectives of this work are summarized below : 

• To explore the directed energy deposition process of ion mixing for producing 
germanide phases in bilayer and multilayer thin film structures of Cu/Ge, Ni/Ge 
and Co/Ge systems at different ambient temperatures, using Kr and Ar ions of 
mega electron-volt energies. 

• To obtain quantitative information of the parameters viz. mixing rates, mixing 
efficiencies, activation energies from the dose dependence, temperature depen- 
dence of mixing (the Q curve) etc. for the above systems. 

• To compare the results of ion induced phase formation in the above metal/Ge 
systems with those obtained by solid state reactions at elevated temperatures, 
especially concerning the formation of the first phase. 

• To study the atomic movements of different species during ion- beam mixing using 
marker layer at the interface. 

• To understand the mixing mechanism responsible for the growth of the compound 
layer across the interface and to develop a formalism to explain the experimentally 
obtained mixing rates and efficiencies. 

For ion irradiations, the 2 MeV \'an de Graaff accelerator at IIT, Kanpur was uti- 
lized. Most of the quantitative information viz. layer thicknesses, composition of the 
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mixed regions and depth profile was obtained by Rutherford backscattering spectrome- 
try[7] using MeV ions from the same accelerator. X-ray diffraction technique was 
utilized to identify the phases produced by ion-beam mixing or by thermal treatment. 


1.6 Organization of Thesis 

A brief review of different aspects of ion-beam mixing is presented in Chapter 2. The 
details of sample preparation and the experimental techniques utilized in this work are 
described in Chapter 3. Chapter 4 describes the MeV Kr and Ar ion induced mixing 
effects, leading to the compound formation in multilayer structures of Cu/Ge, Ni/Ge 
and Co/Ge samples held at various ambient temperatures. It also shows the compar- 
ison of these results with those of thermal annealing experiments. Also, the nature 
of the phases formed in these systems has been examined on the basis of the existing 
rules. Chapter 5 deals with the marker experiments which reveal the information on 
the dominantly moving atomic species during ion-beam mixing and thermal annealing. 
The mixing rates, mixing efficiencies etc. are determined using bilayer configurations 
in Chapter 6. In this chapter a new phenomenological model based on chemical ki- 
netics and applicable to reactive systems ( such as metal silicides, germanides) has 
been proposed to explain the experimentally obtained mixing rates. The work has 
been summarized in the last chapter giving summary, conclusions and scope for future 
studies. 



Chapter 2 


Review of ion- beam mixing 


2.1 Introduction 

In this chapter sorae basic aspects of ion-beam mixing which are relevant to this work 
will be reviewed. These include the influence of collisional and thermodynamic effects 
in the mixing process and temperatxire dependence of mixing. The various empirical 
rules reported in the literature for predicting the nature of mixed regions, will be 
mentioned . Also a brief description of the various phenomenological models used for 
quantitative calculation of different mixing parameters will be given. 

2.1.1 Basic mechanism of mixing 

The basic processes involved in ion-beam mixing are illustrated schematically in Fig 2.1. 
When the energetic heavy ion penetrates a top layer A to reach the bulk material B, it 
loses energy due to collision with target atoms. Due to these collisions the target atoms 
receive sufficient energy to get displaced from their original positions and thus produce 
a displacement cascade[5, 6]. The displacement of atoms occur near the interface of 
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layer A and the bulk material B resulting into a mixed region of A and B. I he effects 
of the vehement collisions of the incident ion with the taigct atoms have boon bioadly 
divided into two categories based on the time scales[14, 1 1] coinpan^d to the time taken 
by the incident ion to come to rest in the target which is of the order of pico-second. 
These are called as prompt and delayed processes. 

In the prompt process, the generation of primary recoils by the incident ion and 
the evolution of a cascade by the secondary collisions occur on a time scale of few pico- 
second. The prompt process tend to the mixing of A and B by direct recoil of the target 
atoms or by generation of cascade of secondary atoms called as cascade mixing. The 
direct recoil process takes the atom of layer A deep inside the bulk B (Fig 2.1b) and is 
called as recoil implantation[14, 15, 50]. The probability of direct recoil is very small 
and it cannot quantitatively explain the amount of experimentally observed mixing of 
A and B in most of the cases[14, 15). In case of cascade mixing[l4j, the secondary 
collisions along the track of the incident ion produce a collision cascade which involves 
many atoms with kinetic energies much smaller than the incident ion energy and the 
multiple relocations of the atoms occur resulting in the mixing of dissimilar atoms 
of A and B across the interface. Such collision cascades are also called as thermal 
spikes[5, 6, 14, 15]. Thus, each incident ion can be assumed to give rise to a small 
volume across the interface containing both atoms A and B (Fig 2.1b). The effect of 
the individual cascades overlap as a function of incident ion dose results in a continuous 
mixed region as shown in Fig 2.1c. Both processes namely the primary recoil and the 
collision cascade are essentially temperature independent. 

In the delayed process, further mixing of A a.ud B occurs due to the thermally acti- 
vated motion of defects produced after the cascade dies. The resulting atomic motion 
persists for a time period exceeding several nanosecond and depends on the tempera- 
ture of the target[14, 15, 21, 23]. Mixing occurring due to these delayed effects follows 
an Arrhenious behavior. It has been experimentally observed that for a binary system 
there exist two critical temperatures TepO, 51] and reij[51]. Above the mixing is 
temperature dependent and below it mixing is temperature independent. The critical 




(a) Initial stage 


(b) Intermidiate stage (c) Final stage 


Figure 2.1: Schematic diagram showing the ion-beam mixing processes. 
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temperature T^q, above which normally stable compound phases form, is known as 
the equilibrium temperature. It is always greater than or equal to 7;.[51]. However, 
there are many reports of equilibrium phase formation at relatively lower temperature 
around Tc in case of some metal silicides and aluminides systems[15. 25, .'H, 41, 42], 
For technological applications it is essential to know the nature of the final phase and 
its stability. In general, amorphous, crystalline metastable or stable phases form after 
ion-beam mixing. To predict the nature of the mixed region (amorphous/crystalline) 
several empirical rules[52-56] have been proposed which will be discussed in the fol- 
lowing section. 


2.1.2 Phase formation 

The phase formation in metal/semiconductor systems by ion-beam mixing can be di.s- 
cussed on the basis of metal/Si systems for which ion-beam mixing have been exten- 
sively studied using both bilayer (unlimited source) and multilayer (limited source) 
structures. The results[17, 31. 42, 57] of these are summarized in Table 2.1. The 
ion induced silicide formation characteristic can be compared with those obtained by 
steady-state thermal annealing. The results of silicide formation[33. 34, 57]) by ther- 
mal annealing are also given in Table 2.1, It can be clearly seen from this table that 
the first silicide phase induced by ion-beam mixing in unlimited supply (bilayer) sam- 
ples is same as that obtained by thermal annealing. Although some of the layers that 
were ion mixed at room temperature do not e.xhibit a distinct crystal structure due 
to low intensity X-ray diffraction pattems[57], the composition of the mixed layers is 
invariably same as that obtained by thermal annealing. This correlation concerning 
the first phase formation by ion-beam mixing vis a vis the first ph<^e formation by 
thermal treatment does not hold true for metal/metal systems[15]. 

In solid state reactions, the kinetic constraints play a dominant role, whereas in 
ion-beam mixing it is normally assumed that the kinetic restrictions are relaxed[27]. 
Thus, both the reaction kinetics as well as the thermodynamic driving forces play 
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Table 2.1: Comparison of first phase formation by ion-beam mi.xing and thermal an- 
nealing in metal/semiconductor systems 


Metal/Semi conductor 

system 

Phase formed by IBM 

Phase formed by 

thermal annealing 

composition 

phase 

Xi/Si 

XhSi 

NkSi 

NkSi 

Co/Si 

C02Si 

C02Si 

C02Si 

Pd/Ge 

Pd2Ge 

Pd2Ge 

Pd2Ge 

Pd/Si 

Pd2Si 

Pd2Ge 

Pd2Ge 

Cr/Si 

^ CrSi2 

- 

CoSk 

V/Si 


- 

VSk 

Ti/Si 

TiSk 

- 

TiSk 

Xb/Si 

NbSk 

NbSh 

NbSk 

Pt/Si 

PtiSi 

Pt2Si 

Pt2Si 


-weak x-ray reflection 
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aa active role during mixing. Extensive efforts have been made experimentally as 
well as theoretically to understand the effect of thermodynamic driving forces[15, 27. 
31]. mobilities[31, 56], composition[6, 13] etc. to predict phase formation during ion- 
beam mixing and thermal annealing in metal/metal aluminides[25, o8, 60, 61] and 
metal/Si[31. 57-60] thin films. 

Several empirical rules have been suggested to predict the formation of crystalline 
compound by ion-beam mixing[52-56]. Similar rules concerning the formation of amor- 
phous phases are also relevant to compound formation since the former can preclude the 
later[61]. Amongst these the most widely used rules are simple structure rule[52], mod- 
ified structure rule[53], negative heat of compound formation rule[62], Hagg’s rule[55] 
and Lau's rule[57]. All these rules have also been applied to metal/ssemiconductor 
systems by several authors. However, In many cases the formation of first and the 
successive phases in metal/Si systems cannot be predicted by the above mentioned 
rules[27]. Recently Pretorius et al.[56] have formulated effective heat of formation rule 
which successfully predicts the sequence of phase formation by thermal annealing and 
ion-beam mixing in many systems. There are also many other factors which influence 
the final phase formation which are not considered in these rules. In general, the factors 
which most likely influence the formation of phase are the structure of the terminal 
solution, complexity of the compound phase, width of phase field, mobility, effect of 
substrate temperature etc. These rules will be discussed in detail in the fdrth chapter 
while explaining the present studies of phase forniation in metal/Ge systems. 


2.2 Critical temperature 

The results of temperature dependents of mixing show that the mixing is divided in 
temperature independent and temperature dependent regimes[14, 15, 21. 23, 51]. As 
uieationed above in most of tlie systems it has been observed that, in general, there 
exist two critical temperatures and T,, with Tc < such that metastable phases 
are expected between Tc and^^Te, and equilibrium phases are expected above Ttq. de 
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Reus et al.[51] have correlated % and Teq with the energies required for creating * holes 
in the compound. The smaller hole formation energy is related with Tc while for the 
larger one with Teg. These temperatures are given by the following empirical relations. 


Tc = 4.0 X AHvsmall 


( 2 . 1 ) 


and 

Teg = 3.8 X AHviarge (2.2) 

where the temperatures are given in K and the hole formation enthalpies in kJ/mole. 
The plots of Tc or Teg (Fig 2.2) as a function of AHv gives a straight line. In several 
systems the difference between Tc and is not so pronounced and in those cases there 
is only one transition temperature[51]. Table 2.2 contains many reported data[51] for 
metal/metaJ and metal/semiconductor systems. The hole formation enthalpies AHv 
for an alloy Ai-xBj; are calculated from Miedemma’s model[63] and are given by 


AHf = caAMA- + (1 - Ci)(yjVBfl^AH^. (2.3) 


and 

AH^ = CBAHPy + (1 - CB)(VBlVj,fl^AHfy (2.4) 

where AHiy and AH^ are the mono vacancy formation enthalpies in pure A and B. 
The molar volume ratio Va/Vb can be regarded as a measure to which degree the atoms 
are surrounded by both constituents. The effective concentrations ca and cb depend 
on X. 

Cheng et al.[20] have correlated Tc with average cohesive energ}' of mixing elements. 
Since the cohesive energy is proportional to the melting temperature of solids, Rossi 
and Nastasi[23] have correlated Tc with the melting temperature. Cheng et al.[20] have 
plotted Tc values obtained experimentally against AHcok a.nd obtained a straight line 
whose slope is 100 K/eV. Thus 


Tc = IQOAHcoh 


( 2 . 5 ) 
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Figure 2.2: Plots showing variation of (a) i; as a function of AHv^maii and (b) Z 
a function of AHviarge (from ref 51 and data given in table 2.2). 


as 
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Table 2.2; Comparison of experimentally obtained critical temperatures with those 
calculated from the model of de Raus et al. (after ref 51) 


Bilayer 

system 

Critical temperature 

TcinK 

Eq. temperature 

Teq in K 

Exp. 

Calculated from 

Exp. 

Cal. from 



1) Al-Cu 

350 

284 

- 

353.4 

2) Al-Fe 

- 

356 


433.2 

3) Al-Mn 

350 

272 


326.8 

4) Al-Mo 

290 

516 

- 


5) Al-V 

- 

408 

420 

448.4 

6) Al-Zr 

- 

412 


497.8 

7) Au-Ni 

400 

384 

- 


8) Au-Si 

400 

292 

- 

326.8 

9) Au-Zr 

450 

504 


566.2 

10) Co-Si 

300 

356 



11) Cr-Si 

370 

! 

336 


368.6 

12) Cu-Ti 


428 

525 

547.2 

13) Fe-PT 

454 

528 

- 

638.4 

14) Fe-Ti 

480 

496 


634.6 

15) Fe-Zr 

480 

468 


744.8 

16) Mo-Ni 

- 

572 



17) Mo-Si 

630 

500 

- 



not available 


continue* • • 
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Cheng has shown that the Tc values of other systems specifically metal/metal sys- 
tems also follow the above relation. However, there are quite a few studies especially 
concerning metal/semiconductor systems, which report experimental values of Tc 
which are significantly different than the one predicted by the above relationship [64]. 


2.3 Phenomenological model 

For quantitative understanding of mixing it is important to estimate some of the useful 
mixing parameters like (i) mixing rate rate at which the squared thickness 

of the mixed region increases as a function of irradiation dose (#) at given temperature 
for a particular incident ion beam, (ii) mixing efficiency {X^I^Fd): normalized mixing 
rate with respect to the deposited energy density Fd- For theoretical prediction of 
these parameters, several phenomenological models have been proposed which depend 
on different experimental conditions viz. atomic number of incident ion and target 
atoms, energy deposited in collisions, nature of the thermal spike (whether local or 
global), thermodynamic parameters, nature of the final phase etc.[19, 22-24. 31]. 

The experimental data of mixing rates in various systems reported in the literature 
reveal that the mixing rate increases linearly with Fd up to certain critical damage 
energy (F^*)[24, 65] above which it increases much more rapidly exhibiting a square 
dependence on Fd . Obviously, the mixing efficiency is constant in the linear region. 
For example[66], in the case of Co/Ti system, the mixing rate vari^ linearly with Fd up 
to 2.25 keV/nm after which it becomes quadratic. For Pd/Si system (shown in Fig 2.3), 
this linearity is seen for Fd values up to 1.5 keV/nm[31]. For systems with Zave > 20 
this dependence comes from the nature of the thermal spikes. For systems having Zave 
< 20 the concept of thermal spike is not applicable[14]. Different phenomenological 
models have been developed to explain the mixing rates in the linear and quadratic 
region of the mixing rate vs. Fd curve. In the next section we shall describe some of 
the models frequently used for the calculation of mixing fates and efficiencies. 
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2.3.1 Ballistic Model 

Sigmund and Gras Marti[50] studied the effect of ion bombardment in a dilute solution 
(< 1%) of an impurity in a homogeneous matrix at low temperature and low current 
density limit. This implies that the mixing occurs due to linear superposition of ele- 
mentary events which are associated with individual ion impacts. This approach has 
been followed by several authors to model mixing rate. The formulae developed by 
Sigmund and Gras Maxti to evaluate the effects of recoil mixing and cascade mixing 
are given below. 

Recoil Mixing 

The effect of recoil mixing can be identified in the RBS spectrum as a tail in the low 
energy side of the impurity peak[50]. The probability of such an event is given by 
relocation cross-section da[Eo{x),T\ where Ec{x) is the energy of the iiTcoming ion at 
a depth x and T < Eq. Where m and M are the masses of incident and target 

atoms. 

The probability density P(z) for animp'urity at depth x to get relocated at depth z 
is given by 


F(Ax)d(Az) = ^da-(x,z) 

where $ is the irradiation dose. The long-range recoil implantation is mainly due to 

single ion-impurity knockon events whose cross section dcr is of the order of the square 

of the atomic radius, which is about 10“^® cm^. Thus long-range recoil implantation 

at typical dose $ < 10^® ions/cm^ is a rare event[14, 50]. 

\ 

Cascade Mixing 

As mentioned above, a cascade is formed when a large number of atoms are displaced 
from their lattice positions and move very small distances in random directions. This 
motion occurs due to the secondary collisions between impurity and matrix atoms. If 
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the impurity concentration is very low (< 1%) or if the deposited energy density (Fd) 
is low then the spread in the concentration profile is given by [14, 50], 

( 2 . 6 ) 

where F = 0.608, and M 2 being the masses of the two target 

atoms, C = average of atomic density, Ec = threshold or cut off energy (8 ev)[14, 15], 
Rc = mean square range at Ec- 

For stable Frenkel pair formation Ec should be replaced[14] by the displacement 
energy F?(f[67] and Rc should be equal to the minimum distance A (~ 1 nm) where 
recombination does not occur. Setting (F = 4Dt, the effective diffusion coefficient can 
be written as 


D 


IdCF- Ir-„,W2A2 


4 


3C 


E, 


(2.7) 


From the above relation we get oc i.e. oc which is similar to the thermal 
annealing process in which 0 a The ballistic mixing rate for a mixed layer of 
thickness X is given by 




( 2 . 8 ) 


The main limitation of this theory is that no general result can be derived from the 
relocation profile without approximation[68, 69]. Also the numerical solutions of the 
transport equations are very complicated. However, this kinetic view point of cascade 
mixing provides a good description of broadening and shift of the concentration profile 
of dilute impurity ( < 1%) in solid matrix. 


2.3.2 Thermodynamic effects 

The above mentioned ballistic approach ignores the thermodynamic effects which are 
important for large atomic concentration. Recently, Miotello and Kelly[26] incorpo- 
rated chemically guided defect motion in the ballistic approach of mixing. Due to this 
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reason ballistic model described above could not properly describe ion-beam mixing es- 
pecially the mixing rates and efficiencies in many systems[14, 15, 70]. For example[Tl], 
if only the ballistic effects are considered the Hf-Ni and Hf-Ti systems should show 
similar mixing behavior. Contrarily, it has been found that the Hf-Ni system mixes 
much more efficiently than Hf-Ti system. For equal atomic concentrations, the heat of 
mixing for Hf-Ni is -62 kJ/g at while for Hf-Ti it is 0 kJ/g at. The difference in the 
mixing was believed to result from the large difference in the heat of mixing between 
these two systems. The large ion-beam mixing effects in the miscible Au-Cu and the 
negligible mixing of immiscible W-Cu system is another such example[14]. It has also 
been realized that thermodynamic driving forces also play a dominant role along with 
the ballistic effects in the formation of various phases during mixing. [14, 15, 70]. On 
the basis of studies reported for various systems the dependence of mixing on AHmix 
and AHcoh can be summarized as follows. 

• The Mixing rate varies linearly with the heat of mixing AJif,n«x[14, 70]. The slope 
of this curve gives the effective temperature of mixing and hence the kinetic 
energy of the moving atom which is found to be of the order of 1 ev/atom. 

• The mixing rates and the mixing efficiencies are very low, moderate and high for 
binary systems having high positive, zero and negative values of heat of mixing 
respectively[14. 15]. 

• The mixing rate varies inversely as the square of Ai/coAllS], 

Several phenomenological models[14, 15, 19, 22, 24, 31] which consider the effect of 
the chemical driving force, nature (overlapping and non-overlapping subcascades) and 
geometry (cylindrical or spherical) of the thermal spike have been developed to explain 
above experimental results. Although there is no direct experimental evidence of a 
thermal spike, the concept of which has been found useful in explaining the ion-beam 
mixing effects. In the following we shcill describe the characteristics of a thermal spike 
followed by brief discussion of different models utilized for calculating mixing rates. 
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2.4 Thermal Spike 

A spike is defined as a limited volume inside a solid within which the majority of atoms 
are temporarily in motion[14]. A thermal spike is a spike where atomic motion reaches 
a state of quasi-equilibrium, so that the law of equipartition of energy holds in the 
spike volume and the velocity distribution can be described by Maxwell-Boltzmann 
statistics. The characteristics of such spike are 

• The average energy density is very high and exceeds the melting temperature. 

• The local temperature is as high as 1000-3000 K. 

• The cooling time of the spike is 0.1 psec to 10 psec. 

The thermal spikes axe distributed at random inside the cascade volume. The 
cascade displacement energy, that is, the number of displacement per unit volume is 
not homogeneous in space. The density is locally very high in the spike and low in 
the other part of the cascade. The density of spikes depends on the amount of energy 
deposition which in turn depends on mass and energy of the incident ion and the 
ballistic properties of target atoms. Depending on the stopping power of the target 
and ion energy, the distance between the primary collision can be larger than size of 
the cascade initiated by the secondary recoils[14] . It is to be noted that a cascade may 
contain several well separate subcascades each of which may have a separate spike or 
local spike. The ideas about spike have been recently quantified by Cheng[14]. On 
the basis of fractal geometry approach Cheng has argued that space filling fractals is 
a necessary condition for spike formation. He has shown that an atom will initiate a 
space filfing cascade or spike, only if its kinetic energy is smaller than E^. This critical 
energy Ec is given by 

ilc = 3.929 X eu (2.9) 

Cheng also calculated that spikes cannot form in a matrix consisting of elements 
with Z values less than 20[14]. He has proposed the following conceptual picture for 
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spike formation: 

1. When jEc is smaller than the threshold displacement energy Ed, a collision cascade 
will stop before reaching the space-filling condition. The formation of spike is thus 
impossible. 

2. When (Ec > Eo > Ed), where Eo is the incident ion energy, the space-filling 
condition is satisfied. Here, the formation of spike is possible. 

3. When {Eo > Ec> Ed), the initial kinetic energy of the incident ion Eo is shared 
by the displaced atoms. Each displaced atom is capable of generating a space- 
filling sub-cascade, or a local spike, if this displaced atom has a kinetic energy 
just below Ec- A global spike will form if local spikes overlap. 

It is to be noted that the formation of spikes takes place only by higher generation 
recoils and is determined by Ec, not the incident ion energj' Eo- 


2.4.1 Global or overlapping spike model 

Johnson et al.[19] have developed this model under the condition of overlapping of the 
subcascades producing a single cascade or global spike for higher deposition energies. 
The development of this model is based on Vineyard’s spike model[72}. In Vineyard’s 
model it is assumed that the energy is deposited instantaneously in a very small vol- 
ume, producing a localized increase of temperature which spreads according to the 
classical laws of heat conduction in a homogeneous medium. Vineyard has calculated 
the jumping rate t] for the atoms in a spike of cylindrical geometry as 


AK^ FI 


( 2 . 10 ) 


Here the heat capacity /c and thermal conductivity C are independent of tempera- 
ture, Q is the internal energy' per unit volume and A is the proportionality constant. 
Johnson et al.{19] modified this formula, assuming that the mixed region is in the form 
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of a regular solution, by empirically including in it the heat of mixing and the cohesive 
energ^^ (A^coa) which is related with Q[18]. The jump frequency t] is then given by[18, 
19] 


V 


AK^Fl 


1 - 




Q 


( 2 . 11 ) 


SttkCQ^ [ 
where AHmix is the heat of mixing. 

Considering the spikes to be independent of each other Johnson et al.[l9] have 
postulated the following equation for calculating the mixing rate. 


_ diAJDt) _ FI r A^l 
$ Wo ^ ^'AifcoA- 


( 2 . 12 ) 


where Ki and are constants, C is the average atomic density. Using this for- 
mula, Cheng[14] has plotted the quantity [d(iDt)/d^] x [C^^^{AHcoh)^ / F^] versus 
AHmix! ^Scoh which is shown in Fig 2.4. Here Cheng has used experimentally mea- 
sured values of [d(4£>t)/d$] and the tabulated values for AHmix and AHcoh reported 
by Miedema[63] and Kittel[73] respectively. From this curve the values of K\ and 
are found to be 0.035 A and 27.5 respectively. The above relation is valid for infinitely 
thick layers where the composition of the intermixed region is not influenced by the 
limited supply of one of the constituents[14, 74]. The experimental values of the mixing 
rates in the region where it exhibits a quadratic dependence on Fd can be explained 
fairly well by this model. 


2.4.2 Local or non-oyerlapping spike model 

The hnear dependence of mixing rate on damage energy cannot be explained by John- 
son et al. model. Alford et al.[74, 75] have found that the mixing rate and efficiency 
do not possess the same energy dependence as predicted by Johnson et al. model and 
IS higher by a factor of about 5. The constancy of the mixing efficiency or the linear 
dependence of the mixing rate on Fd has been explained by Borgesen et al.[22] and 
Bolse[24] by considering well separated local spikes. The higher degree of mixing in 
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Figure 2.4; Plot of ld{4Dt)id^ x [C^f%AHcoh.y / F^] ss a function of AHmixIAHcoh 
(from ref 14). 
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these cases has been taken into account as follows. A local thermal spike would cool 
much slowly as compared to the global spike because the local spikes are isolated from 
each other in a cold lattice. A larger cooling time implies a greater degree of atomic 
motion and hence a larger amount of mixing. 

According to Borgesen et al.[22], for medium Z materials the subcascades and hence 
local spikes are well separated to prevent mutual overlapping. Assuming a spherical 
non-overlapping spike they deduced the mixin g rate as, 




\BoTgesen 

loc 


KiZli^Fp r 5 Ag„i. i 


(2.13) 


where Kb is 0.034 nm-eV^/^ and K 2 = 27.4 are used[75]. Using this relation they have 


explained the experimental mixing rates of many medium Z (20 < < 40) bilayer 

systems. 

However, Bolse[24] has argued that the local spikes are cylindrical in nature and 

not spherical as assumed by Borgesen et al.[22]. According to Boise, the mixing rate 
is given by 


\Boise_ ^bZIi:^‘Fd 
$ i/oc — 


1 + if; 




AR 


coh 


(2.14) 


where Kp = 0.0163 nm^ eV and A', = 41 are empirical constants obtained from 
experimental results[76]. 

Boise has also estimated the critical damage[24] energy density Fff, for a given 
ion-target combination, above which the nonlinear dependence of mixing rate on Fd 
begins. In other wards the subcascades begin to overlap. The Fg is given by 


^ = Kr^C 


ZiZt ] 

m{ 1 

^P.23^20.23 

Im,- -f AfJ 


(2.15) 


* ui tae mciuent ion, Zj and Mt = 

atomic nnmber and mas, of the targe, atom and = 2.5 keV „m=[24, 651 is the 

proportionality constant. 
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All the above mentioned models satisfactorily predict the mixing rates for most of 
the metal/metal systems. However, many metal/semiconductor systems where growth 
of a compound occurs by ion-beam mixing exhibit high value of mixing rates which 
could not be predicted by the above formulas (Eq 2.12 or Eq 2.13 or Eq 2.14). For 
example, the recent results of Pd/Si bilayer system (shown in Fig 2.5) cannot be 
explained by the above models. These metal/semiconductor systems are generally 
considered as a highly reactive systems. The models which explain the formation of 
phases or compound axe described in the next section. 


2.5 Model for Compound Formation 

Desimoni and Traverse[31] have proposed an ion-beam mixing model that accounts for 
compound formation at a boundary between two materials during ion irradiation. Their 
model is based on Fick’s law which incorporates a chemical driving force to simulate 
the chemical reaction at the boundary. Ihe model repuduces the two fluence regimes, 
linear and quadratic found in several ion beam experiments. When the square4 of 
the mixed thickness (X^) is plotted as a function of the irradiation fluence ($) the 
plot exhibits a quadratic variation of X^ until a critical dose $c is reacted above 
which it shows a linear dependence. This is illustrated in Fig 2.5 for Pd/Si systems. 
The critical fluence in this case is 1.5 x 10^® ionsjcm^. According to Compound 
formation model of Desimoni and Traverse the critical dose of $c is proportional to 
where k is reaction rate and D is diffusion coefficient under ion irradiation. For 
highly reactive systems for which k is very large . The critical fluence $c is quite low 
(< 10^^ ionsfcm^ ) . Therefore in these systems the linear dependence of X'^ on $ 
is only seen. On the other hand, for systems for which k is low the $c is quite higl^ 
therefore the quadratic dependence is only observed. For example in Cr-Si case the 
k rate is 1 x 10“^ nm/s and the value turns out to be 8 x 10^® ionsfcm^. This 
quadratic regime is called reaction control process while the linear regime is called 
diffusion control process. When diffusion is the limiting process, D is enhanced by a 
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Figure 2.5: Squared thickness of the mixed layer vs the fluence in Pd/Si irradiated at 
80 K with 730 keV Au ions (after ref 31). 
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factor which accounts for the formation of a compound in the mixed layer. Desimoni 
and Traverse have incorporated this enhancement of D in the ballistic and global spike 
models of mixing mentioned above and reported the following relations for calculating 
mixing rates. 
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Here C is the atomic density of compound AaBh formed across the interface. Ca sJ^d 
Cb are the atomic densities of layers A and B respectively. 


2.6 Formalism based on kinetics 

Miotello and kelly[26] have followed a different approach to address the compound 
formation in ion-beam mixing. The phase formation has been viewed as a kinetic 
rather than thermodynamic phenomenon. They have developed a general formalism 
to describe ion-beam mixing by introducing a chemically guided defect motion to obtain 
the diffusion flux J,- of a component i in a bilayer system. The diffusion equation 

^ = A'fr 

dt ox^ 

where D\ = (Z)* -f- -f- Z?f), can be numerically solved for a bilayer geometry. 

Here {D\, Z?”^ and Z?f) are diffusivities due to ballistic mixing, chemical non-guided 
and guided mixing. In order to reproduce the different region of the experimentally 
obtained concentration profile the diffusivity values are adjusted according to the binary 
phase diagram. The diffusivity must be high when AZf^ is strongly negative, and low 
when strongly positive. Also the value of D\ changes significantly at the boundaries of 
the phase. The growth of a phase is ^ven by, 
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ft’here a* is the concentrate of the phase and Df~^ and indicate diffusivity at 


the two sides of the growing layer. This is a generalization of the model presented by 
Desimoni and Traverse which is described above. 



Chapter 3 


Experimental 


3.1 Introduction 

The ion-beam mixing studies presented in this thesis typically involved the following 
experimental steps : (i) preparation of clean substrates (ii) vacuum deposition of ele- 
ments to be mixed in the desired geometry (iii) irradiation of the sample kept at a fixed 
temperature, in the range of 77K to 500 K, with MeV Ar or Kr ions, (iv) compositional 
analysis by Rutherford backscattering spectrometry (RBS) (v) surface morphological 
study by scanning electron microscopy (SEM) and (vi) identification of phases by X- 
ray diffraction technique. In some cases, to study the thermal stability of the mixed 
region, the ion-beam noixed samples were annealed in vacuum and steps (iv) to (vi) 
were carried out. On the other hand, for comparing the effects of ion-beam mixing 
with thermally induced reactions, step (iii) was replaced by thermal anne^ng of the 
sample in vacuum. 

The 2 MeV Van de Graaff accelerator laboratory at IIT, Kanpur was used for RBS 
measurements using ion as well as ion-beam mixing experiments using Ar and 
Kr ions. For performing mixing experiments in the temperature range of 77-500 K, 
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an ion irradiation set-np was designed and fabricated along with other ancillaries. A 
water cooled multi-hearth electron gun evaporation system was used for the deposition 
of thin films and the heat treatments were performed in high vacuum. 

This chapter describes the experimental steps mentioned above along with the de- 
scription of the accelerator facility. 


3.2 Sample preparation 

3.2.1 Sample cleaning 

The detailed studies of phase formation and the mechanism of ion-beam mixing requires 
mainly three kinds of sample configurations (described in Chapter 1) namely multilayer 
(Fig l.la)^bilayer (Fig 1.1c) and marker (Fig Lie) thin films of two elements on inert 
substrate. In the case of ion-beam mixing studies of metal/Ge systems, experiments 
were performed using all the three geometries Experiments involving movements of 
atomic species may be significantly influenced in the presence of foreign atoms and 
hence it is essential to prepare as clean a saimple as possible[77]. Substrates of size 10 
mm x 10 mm were cut from a high purity quartz plate of thickness 0.5 mm. Standard 
cleaning procedures were employed using high purity acetone, tricloroethylene and 
ethyl alcohol to de-grease the surface of the substrate using ultrasonic cleaner. 

3.2.2 Thin film deposition 

Thin films were deposited using a multi-heeurth electron gun evaporation unit shown 
in Fig 3.1. It is equipped with a liquid nitrogen trap to have a clean vacuum in the 
deposition chamber. The working principle of the electron-beam apparatus involves 
the generation of thermionic electrons from a tungsten filament and the bending of 
these electrons towards the water cooled graphite crucible which contains the material 
to be deposited. The thermionic electrons are generated by passing current through 



35 


a tungsten filament held at a negative potential of few kV. The emitted electrons are 
bent and focused with the help of a permanent magnet on to the material. The unit is 
also equipped with a quartz crystal controlled digital thickness monitor for in situ mea- 
surement of film thickness. The monitor was independently calibrated by measuring 
the thickness of several films by RBS measurements [7]. The major advantages of this 
deposition system are the accurate control of rate of deposition which is an important 
parameter associated with film adherence[77] and the ability to perform deposition of 
a number of layers of different materials without breaking the vacuum which provides 
a clean interface between- two layers. 

High purity (99.99%) material of Ge, Cu, Ni and Co were used for preparing mid- 
tilayer or bilayer films of M/Ge configuration on quartz substrate where M is one of 
the above mentioned metallic species. A thin Au layer was used as a marker at the in- 
terface of M/Ge bilayer sarnples. In order to make the Kr and Ar ions pass completely 
through bilayer or multilayer structures, the total thickness of the film was always kept 
in the range of 100-300 nm. which is much lower than the range of 1 MeV Ar or Kr ions 
in these materials [78]. A number of «f samples were prepared for the present studies 
of ion-beam mixing of metal/Ge as described in Table 3.1. 


3.3 Van de Graaff accelerator facility 

The Van de Graaff accelerator is the most widely used among the commonly available 
electrostatic accelerators. A schematic diagram of the accelerator system is shown in 
Fig 3.2, and photographs of the same are presented in Fig 3.3a and Fig 3.3b. The 
accelerator (Model AN-2000), High Voltage Engineering Corporation make USA, and 
is a precision 2-MeV high intensity source of positive ions. The accelerator produces 
a beam of ions which is homogeneous, stable and controllable over a wide range of 
energies. It is pro\dded with three gas bottles and presently they are filled with three 
different gases namely He, Ar and Kr. An r.f. ion source is used to produce positive 


ions. 
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Figure 3.1: Photograph of multi-hearth electron-beam 


gun evaporation system. 
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Table 3.1: Description of the samples used for present studies 


Geometry 

geometry 

Sample 

Thickness of individual 

layers 

No. of metal 

&: Ge layers each 

Metal/Ge 

multilayer 

films 

deposited 

on Si 02 

5i02/Cu/Ge 

each Ge layer thickness 11-15 nm 

each Cu layer thickness 23-27 nm 

5 

^iOs/Ni/Ge 

each Ge layer thickness 16-26 nm 

each Ni layer thickness 32-40 nm 

5 

5i02/Co/Ge 

each Ge layer thickness 22-27 nm 

each Co layer thickness 22-25 nm 

4 

Metal/Ge 

bilayer 

film 

with a Au 

layer as a 

marker 

5i02/Cu/ Au/ Ge 

Ge layer thickness 80-85 nm 

Cu layer thickness 23-27 nm 

2 

(Thickness of Au 
marker < 1 nm) 

SiO^I Ni/Au/Ge 

Ge layer thickness 120-130 nm 

Ni layer thickness 65-70 nm 

^iOa/Co/Au/Ge 

Ge layer thickness 95-100 nm 

Co layer thickness 55-60 nm 

Metal/Ge 

bilayer 

films 

diposited 

on SiO^ 

5i02/Cu/Ge 

Ge layer thickness 110-170 nm 

Cu layer thickness 50-100 nm 

2 

Si02l Ni/Ge 

Ge layer thickness 95-105 nm 

Ni layer thickness 55-70 nm 

^iOj/Co/Ge 

Ge layer thickness 110-120 nm 

Co layer thickness 50-60 nm 
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Figure 3.2; Lay out of the Van de Graaff accelerator facility at I.I.T Kanpur. 
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Figure 3.3: Photographs of (a) Vaa de graaff accelerator and (b) the analyzing magnet. 
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The mass and energy analysis of the accelerated ion-beam is done by an analyzing 
cum switching magnet which is a 16 KGauss electro- magnet having a d.c. regulated 
(0-30 V) power supply which has a capability of supplying current up to 250 A[79]. 
During operation the magnet is cooled to 10 "C, by circulating de-mineralized water 
through the coils. However to operate the magnet at currents of less than 140 A, air 
cooling is done by passing air instead of water through the tubular aluminum coils. 
The stainless steel switching chamber has five ports at angles 0° , ±20° and ±45° . The 
energy and mass analyzed beam can be directed to the experimental chamber through 
any one of the five ports. 


3.4 Ion irradiation 

The detailed physical processes occurring during irradiation have already been dis- 
cussed in the previous chapter. Here we describe the beam line containing important 
ancillaries such as sample manipulator, beam stopper, beam defining slits, beam pro- 
file monitor etc. which were developed and installed (for this work) for carrying out 
uniform irradiations on a sample in the temperature range of 77- 500 K. 


3.4.1 Heavy ion-beam line 

The accelerated ions obtained from the accelerator reach the experimental chamber 
through a clean and highly evacuated beam line as shown in Fig 3.4. The beam line 
has a specially adapted liquid nitrogen trap placed in between the accelerator and the 
analyzing magnet to achieve clean vacuum in that region. The ion- beam after leaving 
the analyzing magnet passes through beam stopper foUowed by a beam profile monitor. 
The beam stopper is made up of a ‘flat plunger’ which moves in or out of a coil when 
It IS energized or de-energized. The plunger can thus stop the beam from falling on the 
sample when a predetermined dose on the sample is achieved. The beam profile monitor 
(. Model 522, DANFYSIK) is used to monitor the position, shape and uniformity of 
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the beam. The beam then passes through a specially designed beam defining X-Y 
slit arrangement mounted on a double cross chamber (as shown in Fig 3.4a). This 
arrangements is used for the purpose of selecting a uniform portion of the beam of 
required size onto the target placed in the irradiation chamber. 


3.4.2 Heavy ion irradiation chamber 

‘The chamber used for irradiation is shown in Fig 3.5. The chamber is mounted on 
a turbo molecular pump (BALZERS) to obtain a hydrocarbon free clean vacuum. In 
the center of the chamber there is a movable sample holder as shown in Fig 3.5 on 
which a number of samples can be mounted and by moving the holder vertically by 
a predetermined distance, the required sample can be introduced in the path of the 
beam. The holder is electrically isolated from the rest of the chamber. A current 
integrator (ORTEC model 439) which measures the total charge accumulated during 
irradiation is connected to the sample holder. 

The sample holder is surrounded by a cylinder made up of aluminum sheet which 
has an opening for the passage of the heavy ion-beam. This cylinder is maintained 
at a poetical of nearly -260 volt relative to the grounded chamber. The arrangement, 
called as secondary electron suppressor, inhibits the interference of secondary electrons 
emitted from the sample surface during ion irradiation. 

The arrangement described above is suitable for room temperature irradiations. 
For performing irradiations at high and low temperatures, specially designed sample 
holders were fabricated. Figure 3.6a shows the photograph of the sample holder along 
with the heating arrangements for the high temperature irradiations. This sample 
holder is made up of high purity Copper for good thermal conduction. The heat loss 
by conduction via the sample holder is prevented by connecting the sample holder 
with a ceramic insulator. Temperature measurements are done by placing a chromel- 
alumel thermocouple junction as close to the sample as possible- With the help of a 
temperature controller the temperature of the sample can be maintained to within ±3 
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Figure 3.4: Photographs of beam line and ■ irradiation chamber. 
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K of the set temperature. Mica sheets are used to electrically insulate the thermocouple 
junction from the sample holder. The heating arrangement designed for heating the 
sample up to 550 K consists of a tungsten filament lamp (25 V, 10 A) which is placed 
on the sample holder to give uniform heating. A reflector made of stainless steel sheet 
is mounted to focus the light on to the sample. The reflector also serves the purpose 
of protecting the chamber walls from getting heated. 

The sample holder fabricated for low temperature experiments is shown in Fig .3.6b. 
It is made up of a rectangular copper block and a stainless steel container which store 
liquid nitrogen. The Copper block on which samples are mounted comes under direct 
contact with the liquid nitrogen stored in the stainless steel container. A copper- 
constantan thermocouple is attached with the copper block to measure the tempera- 
ture. The sample holder can reach a minimum temperature of 90-100 K and stays for 
about 2 hours once the trap is filled with liquid nitrogen. 


3.5 Thermal annealing 

A photograph of the annealing set-up is presented in Fig 3.7. The furnace used for 
annealing has a ceramic tube (i.d.=2 in.), wound with nichrome wire. A uniform 
temperature (maximum 1300 K) zone of 4 in. length can be obtained in the central 
portion of the ceramic tube. The furnace temperature is measured by a chromel- 
alumel thermocouple. The furnace temperature can be controlled to within ±2 K 
using a temperature controller (Indotherm Model 400). The furnace is kept on a 
horizontal movable platform. The sample to be annealed axe kept in the quartz tube 
(i.d.=l in.) which is continuously evacuated by a diffusion pump equipped with a 
liquid nitrogen trap. After reaching the desired temperature, furnace is moved, such 
that the quartz tube containing the samples enter in the central portion of the furnace. 
The temperature of the sample is measured with a chromel-alumel thermocouple kept 
directly m contact with the samples in the quartz tube. To achieve thermal equilibrium, 
on an average, it takes three minutes during heating and six to ten minutes while 
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Figure 3.6: Photographs of the sample holders used for (a) high temperature ion-beam 
mixing (b) low temperature ion-beam mixing. 
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cooling the sample. The thermal annealing was done at various temperature ranging 
from 400-700 K for duration of 30-300 minutes. 


3.6 Rutherford Backscattering Spectrometry 

Rutherford Backscattering Spectrometry (RBS) is a well established method for de- 
termining the composition at the interface of two layers, identifying the mass of the 
impurity and host atoms and thickness of the films etc. RBS is well suited for charac- 
terization of thin films or near surface region of bulk samples specially when the sample 
contain atoms of heavier species. It is based on the framework of the discoveries of 
Rutherford [1911] and Geiger and Marsden [1913]. In general, He'^ etc. ions in 
the energy range of 0.5 to 3 MeV are used for RBS. The well collimated ion beam 
of size 1 mm^ is made to impinge the surface of the sample making a suitable angle 
of the surface normal as shown in Fig 3.8. The backscattering particle are detected 
using a semiconductor detector subtanding a small solid angle. The signals and energy 
analized using conventional nuclear electronics. The frequency of occurrence of the 
scattering as a function of energy of the backscattered particles, known as "spectrum” 
is recorded by performing multi-channel analyzer. 

The basic physical phenomena occurring at the target are:, 

1. The process of energy* transfer from the particle to the target nucleus in an elastic 
collision. This leads to the concept of kinematic factor (A'e) and the capability 
to determine the mass of the target nucleus. 

2. The process of probability of occurteace of a two-body collision, which leads to 
the concept of scattering cross-section (<r) and the capability to quantitatively 
analyze the atomic composition of the target. 

3. The process of average energy loss of the projectile ions through the target, 

leading to the concept of stopping cross-section (e) and to the capability of depth 
perception. 
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Target atom (M) 



Figure 3.8: Schematic representation of an elastic coIUsion process in in laboratory 
frame. Taxget of mass M is initially at rest. 
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4. The process of statistical fluctuations in the energy loss of the projectile in the 
target. This phenomenon leads to the concept of energy straggling (il') and to a 
limitation in the ultimate mass and depth resolution of backscattering spectrom- 
etry. 

The RBS measurements are usually carried out in high vacuum (pressure in the 
range of 10“® to 10“® Torr). An ultra high vacuum may be needed if prolonged mea- 
surement on the same spot (to prevent carbon build up on the target due to the breaking 
up of residual hydrocarbons in the chamber under ion bombardment) is required. The 
method haa been described in detail by Chu et al.[7]. 


3.7 The RBS setup 

The complete backscattering spectrometry set-up used in the present study consists of 
the following main stages 

1. Production of ions and acceleration using the 2-MeV Van de Graaif accelerator 

2. Mass and energy analysis of the accelerated ions using an electromagnet. 

3. Focusing of the energy analyzed beam and separation the O"*"*" ions from the 
ions of the same energj’ 

4. Experimental setup for backscattering analysis. 

3.7.1 The beam line 

The backscattering set up is situated at the -b45° beam port. After deflection the 
beam passes through energj’ stabilizing slits and quadruple magnets. This focused 
and collimated ion-beam then passes through an uniform electric field generated by 
parallel plates (electrostatic beam separator) in order to separate out from the 
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beam[80] of the same energy. The ion-beam collimating unit placed between 
the RBS chamber and the ion-beam separator has horizontal and vertical slits, which 
define the beam size, so that a beam spot of required dimension (a minimum size of 
0.5 mm X 0.5 mm) be made incident on the target. Beam viewers and beam stoppers 
are placed strategically along the beam path to see the ion-beam for diagnosing the 
beam alignment problem and to block the ion-beam from entering the chamber as and 
when required. 

The backscattering chamber used in this work has six ports as shown in the Fig 3.9 

and Fig 3.10. The sample holder is momited on the top of a flange in such a way 

that the ion-beam entering from one of the ports is incident perpendicularly on the 

sample surface. The sample holder is a plate of dimensions 130 mm X 25 mm X 2.5 

mm. Five to six samples can be mounted on it at a time. It can be moved in the 

vertical direction, such that the region of interest of a sample can be brought in front 

of the incident beam for analysis, with an accuracy of ±0.5 mm. In this way several 

measurements on difilerent regions of the same sample and measurements on various 

samples are possible ^vithout breaking vacuum. An aluminum grid is placed around the 

holder (inside the chamber) and a negative voltage of -260 V is applied to it, to suppress 

the secondary electrons emitted from the samples due to the primary ion-beam. Beam 

currents of 5-10 nA were utilized for the present analysis. Currents more than 10 nA 

increased the detector dead-time and pulse pile-up effect, and hence it was maintained 

at all times below 10 nA. The target or beam current was measured and integrated 

to obtain the total fluence or charge (which relates the number of incident Hc+ ions 

incident on the target) using a current integrator (Ortec, Model-401). Backscattering 

spectra, in the present study, were obtained for a total charge collection of 2-10 /iC for 
each target. 

he partially passivated ion implanted or silicon surface barrier detector is placed 
at a scattering angle (9) erf 150" to the inddent beant. and ia fixed at one of the porta 
of the chamber. The e.xperimental geometry for RBS is shown in the inset of Fig 3.9. 
e detector surface area is decreased by about 80% by inserting a plate having a small 
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RBS 



Figure 3.9: Block diagram of RBS chamber along with target, detector and other 
accessories. 
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Figure 3.10: 



Photograph of RBS chamber along with target and other accessories. 
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aperture infront of the detector. The detector to sample distance was approximately 
10 cm. The solid angle of 2.57 mSr subtended by the detector was measured by taking 
several RBS spectra of of standrad samples. 


3.7.2 Pulse shaping electronics and multichannel analyzer 

The signals from the detector are initially amplified by a charge sensitive preamplifier 
[Ortec Model-142]. The bias circuit (in the preamplifier), for the detector is dc-coupled 
as it has very high resolution. The preamplifier is specially designed to accept the 
detector signal and amplify it with some shaping so as to preserve the maximum signal- 
to-noise ratio. In backscattering, the useful information in the preamplifier output 
signal is the amplitude of each pulse. The subsequent main amplifier [Ortec Model- 
572] amplifying this signal creates a suitable pulse shape (rise time of 1 [iC) thus 
optimizing resolution and count rate capability. The pulse shaping is done to avoid 
pulse pile-up and to enhance signal-to-noise ratio. Each individual detector signal has 
to be terminated in a shorter time compared to the time interval between two pulses. 
Otherwise pulses will overlap causing pulse pile-up, which leads to erroneous amplitude 
measurements. 

The CR-RC (differentiating-integrating) pulse shaping circuit of the amplifier op- 
erates with time constants much shorter than the decay of the preamplifier signal and 
much longer than its rise time. This effectively removes the slow component of the 
preamplifier signal preventing the pile-up error. It also removes the low and high fre- 
quency signals and noise components which significantly enhance the signal-to-noise 
ratio. Hence the amplifier finally produces individual pulses whose amplitudes convey 
the quantity of interest, i.e., the energy of the detected particles. The sigoal output 
from the main amplifier is then fed to a biased amplifier [Ortec Model-444]. With the 
help of this, the amplifier signals are compressed or stretched and the edges are cut off 
as desired, for more efficient analysis in the multichannel analyzer (MCA). The signals 
from the biased amplifier are then collected, either in the MCA [Model-ND65, Nuclear 
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Data Inc.], or a PC-based MCA. The MCA has two distinct data analysis modes: the 
pulse height analysis (PHA) mode and the multichannel scaling (MCS) mode. The 
PHA mode is exclusively used for RBS measurements. 


3.7.3 Calibration of the MCA 


The MCA is calibrated by determining the calibration factor (SF) from the chan- 
nel position of the elements from the experimental RBS spectrum. For this reason, 
backscattering spectra of pure elements (or calibration samples) were taken in every 
run. In each run three to four metal/Ge samples were analyzed along with the calibra- 
tion samples. A set of calibration spectra are depicted in Fig 3.11, for thin film of pure 
Au on pure A1 substrate and pure Al. In these spectra, the surface position of the Au 
is at channel number 494 and that of Al is at channel number 278. The backscattered 
energy for each of these elements (e) is calculated from the following equation 


£'i,e = A'i,e X 


(3.1) 


where Al, is the kinematic factor of the element e and F, is the incident energy of the 

ife+ beam (m this case 1.34 MeV). The calibration factor can be calculated from the 
following equation 


6F = keV/ch 

(cAicAa) 

considering two calibration samples (ej and ea) at a time. If we calculate for th. 
case the calibration factor is obtained as 2.212 keV/ch. 


(3.2) 

present 


3.7.4 Backscattering spectra analysis 

A backscattering spectrum obtained from the same run as that of calibration spectrum 
shown in Fig 3.U. for thin bilayer film of Cu and Ge on quartz substrate is shown in 
F.g 3 11 Th.s spectrum was obtained using He ions of 1.34 MeV at a scattering angle 
(9) of 150». the detector soUd angle (fl) equal to 2.57 mSr and for a beam dose of 5 
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[iC. This RBS parameter description is also given below the figure. It is useful to plot 
the backscattering yield by normalizing it with respect to some of the above mentioned 
parameters, by using the following equation 

A ormahzed y,dd = ^ ^ ^ ^ (3.3) 

where ’raw counts’ are the actual counts recorded; ’corr’ is the correction factor to 
compensate for inaccuracies during charge integration; Q is the beam dose in //C; dfl 
is the solid angle in milli-steridian (mSr); and 6E is the calibration factor of the MCA 
in keV/chaimei. The normalized yields are utilized for comparison of experimental 
spectra for which beam energy and scattering geometry are same but the dose solid 
angle and calibration values are different. The information about the thickness of the 
film and the depth distribution of the diffusing element (in the case where diffusion 
has taken place) can be obtained from such a spectrum. The backscattering spectra in 

the present study have been analyzed by RUMP simulation package[81] as described 
in the following. 

Recently, RBS simulation programs using computers have been developed at several 
places. One such simulation program developed at the Cornell University by Doolit- 
tle[811, has been used in the present studies tor the analysis of the backscattering 
spectra. The simulation program makes use of almost the same eqnations which have 
been described By Chu et al.[71. In simulation, a theoretical sample structure consist- 
ing of several layers of varying thickness is first prepared. The experimental parameters 
are also fed mto the program. A backscattering spectrum is then constructed using 
this sample structure and the experimental parameters. This theoretical profile is then 
compared with the experimentally obtained RBS spectrum. This procedure is utilised 
m an .terahve manner to krive at a best set of parameters (sample structure) using 
w 1C t e given sample spectrum compares quite well with the experimental spectrum. 

^ e WIU now dlustrate ae analysis of the RBS spectrum of Fig 3.12. by the simulation 

Th^ I ♦ ^ layer structure of the sample is given below the Fig 3.12. 

^he fat t« l^ers .mprising of Cn and Ge having thicknesses 660 nm and 1240 nm 
respectively and are further divided into snblayers of 2 nm each. The third layer is the 
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substrate{quartz of thickness 1000 nm, which is further divided into sublayers 

of 50 nm each. The program requires this information along with the experimental 
parameters (as defined below Fig 3.12). The simulated profile is then compared (as 
shown in Fig 3.12), with the experimental profile. One of the important features of 
this simulation programme is that it is an extremely powerful tool for understanding 
complex spectra involving overlap of numerous elements. It also takes into account 
the detector resolution as well as energy straggling of the ions. Thus the thick target 
yield for quartz has been simulated very correctly giving the analysis of the substrate 
as well. Such analysis of the multi-elemental target as mentioned above is not possible 
by hand calculations. 


3.8 X-ray measurements 

The nature of reacted layers after ion-beam mixing and thermal annealing were also 
investigated by the X-ray diffraction (XRD) method to supplement the findings and 
analysis of Rutherford backscattering spectroscopy. The XRD measurements were done 
using a X-ray powder diffractometer [Model ISO-Debyeflex-1001], generating Cu — Kq 
radiation along with a graphite monochromator attached to the goniometer. Some of 
the samples were measured by glancing angle x-ray diffraction technique in Seeman- 
Bohlin geometry using CuKa radiation from University of Poona. A diffractometer 
in general consists of the following main units: (i) power supply; (ii) X-ray generator; 
(iii) goniometer, with monochromator, scintillation counter or proportional counter 
and preamplifier; (iv) amplifier unit; (v) rate meter, and (vi) X-Y chart recorder. The 
spectral lines obtained in the X-ray spectra of the samples in units of 2^ is converted 
into d-values using the Bragg’s rule i.e., 

2dsm = nA (3.4) 

On comparing these values with standard tables of d-values of various alloys and phases 
available in literature, the presence of the suggested or probable phases were confirmed 
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in the samples. 


3.9 Scanning electron microscopy 

In Scanning Electron Microscopy (SEM) the area to be examined is irradiated with 
a focused electron beam, which is swept in a raster across the surface of the speci- 
men. The types of signals which are produced when the focused beam impinges on 
a specimen surface include secondary electrons, backscattered electrons, characteristic 
X-rays, Auger electrons ajid photons of various energies. They are obtained from spe- 
cific emission volumes within the sample and are used to measure many characteristics 
of the sample including surface topography. 

The Auto Scan model JEOL-JSM 840A having 10 nm resolution, about lOX to 
.300000X magnification and operating voltage up to 40 KV was used in secondary 
electron mode. 


3.10 Four probe resistivity measurement technique 

The four point probe comprises four equally spaced metal probes which are pressed 
onto the surface under measurement- The probes are usually set in-line or at the corner 
of a square* For the present work a square probe was used.The sheet resistivity of an 
infinite thin layer in which the current flow is parallel to the surface is [82] 

Ps'^^ = X {VzA/ha) (^■^) 

The resistivity measured by a square probe with the current I flowing between two 
probes and the voltage V measured across the other two probes is given by 

„ _ „» 2,34 Y 

P = Ps 

= 4,532(K3.4//i.2).(X) 


(3.6) 
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where X is the thickness of implanted layer. This formula can be used for calculating 
the resistivity of the modified layers because it meets the requirement of X/S < 0.4 for 
probe spacing S = 1.5 mm. 



Chapter 4 


Ion induced phase formation 


4.1 Introduction 

This chapter^ describes the mixing eflFects, leading to the formation of compounds, 
induced by Ar and Kr ions of megaelectron volt energy in multilayer structures of 
Cu/Ge, Ni/Ge and Co/Ge samples, held at various ambient temperatures in the range 
of 100 - 450 K. For comparing the ion induced mixing effects with those induced by 
thermal reactions, results of thermal annealing experiments performed at temperatures 
in the range of 400 - 700 K are also presented. We then examine the nature (amorphous/ 
crystalline) of the ion induced phases in these systems in the light of the existing rules, 
which axe based on composition and crystal structure of the reacting elements, ratio 
of the atomic radii, equilibrium phase diagram and heats of formation of the possible 
stable/metastable phases. 

^Part of this work is published in Appl. Phys. Lett. 67 (1995) 1700. 
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4.2 Experimental 

In preparing multilayer samples (as shown in Fig.l.la), alternate layers of metal M (M 
= Cu, Ni, Co) and Ge were deposited sequentially on clean quartz substrates with- 
out breaking vacuum by e-beam evaporation system equipped with a digital thickness 
monitor as described on chapter 3. The deposition of the films was controlled to 0.4 
- 0.6 nm/s to obtain uniform film of good quality. In all the cases four to five layers 
of metal M and Ge each were deposited. The thicknesses of the layers were adjusted 
such that the average film compositions range from about 60 at% to 80 at% of metal 
M. The average individual layer thickness of metal M was not more than 40 nm with 
total thickness of 90-300 nm. 

The bilayer structures of metal M and Ge were prepared on quartz substrates with 
Ge as the top layer. The thickness of individual layer was chosen such that it yielded 
the composition of about AfsoGeso. The total thickness of bilayer films was in the range 
of 100-200 nm. 

Ion-beam mixing was done at various substrate temperatures ranging from 100 to 
450 K using 1 MeV Kr and Ar ions. The temperature of the substrate was controlled 
to ±3 K during irradiation. Irradiations were performed for several doses in the range 
of (2 - 20) X 10^® ions/cm^. The ion flux was kept low ( < 1 x 10^^ ) 

to prevent the sample temperature rise during irradiation. A uniform portion of the 
beam was selected by a X-Y beam defining slit arrangement such that the samples got 
irradiated uniformly over an area of 5 mm x 5 mm. The beam position and uniformity 
of the beam were monitored in situ by a beam profile monitor. The projected range of 
the incident ions were calculated using TRIM program[78]. For the incident ion energj- 
of 1 MeV, most of the implanted ions pass through the multUayer and bilayer structures. 
The thermal annealing in vacuum was done for different annealing durations at various 
temperatures in the range of 400 - 700 K. The amount of mixing, the composition of 
the ion mixed or thermally reacted layers were determined by theoretically simulating 
(RUMP code)[Sl] the experimental Rutherford backscattering (RBS) spectra obtained 
using 1.1-1.5 MeV beam. The RBS spectra of the treated samples did not show 
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any visible oxygen or carbon concentrations. The phases were identified by standard 
X-ray diffraction (XRD) using & Cu — Ka source. A few X-ray measurements were 
also performed in the Seeman-Bohlin geometry using the glancing angle XRD facility 
at University of Poona[83]. The surface topography of the samples was studied under 
scanning electron microscope before and after mixing. The resistivity measurements 
were done by Van der Pouw four-probe method[82]. 


4.3 Results 

The results of ion-beam induced effects at different temperatures and those induced by 
thermal annealing alone are presented separately for each metal- Ge system in multilayer 
and bilayer configurations in the following subsections. 


4.3.1 Cu/Ge system 
Ion-beam induced effects 

A typical experimental RBS spectrum (open circles) of as-deposited Cu/Ge multilayer 
sample is shown in Fig 4.1. The spectrum of the as-deposited sample shows a multi- 
peak structure. The peaks corresponding to individual Cu and Ge layers are not com- 
pletely resolved because of the small mass difference between Cu and Ge atoms. The 
surface positions of elements are shown by arrows. The RBS spectra of all as-deposited 
samples were simulated using layer structures of the kind5i02/C'u/Ge/C'u/---C'u/Ge 
as the input to the simulation programme RUMP till best visual fits were obtained in 
each case. One such theoretical simulated spectrum is shown in Fig 4.1 by continuous 
line and the detailed simulated structure is given in Table 4.1. This sample consists 
of five pairs of copper and germanium layers with an overall thickness of Cu and Ge 
equal to 130 nm and 75 nm respectively. 
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Figure 4.1: The as-deposited RBS spectrum of Cu/Ge multilayer sample deposited on 
SiOr. The simulated spectrum was obtained using the sample structure given in Table 
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Table 4.1: Description of the layer structure used for the simulation of RBS spectrum 
of as-deposited Cu/Ge multilayer sample shown in Fig 4.1. 


Layer structure used for simulation 


Sample 

description 

Layer 

no. 

Thickness Composition 
(nm) 

Cu/Ge 

1 

15 

Ge 

1 

Multilayer 

2 

23 

Cu 

1 

films 

3 

15 

Ge 

1 

deposited 

4 

23 

Cu 

1 

on Si 02 

5 

11 

Ge 

1 


6 

23 

Cu 

1 


7 

17 

Ge 

1 


8 

25 

Cu 

1 


9 

17 

Ge 

1 


10 

37 

Cu 

1 


11 

1500 

Si 

1 0 2 
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The effect of ion-beam mixing at room temperature (RT), induced by Kr ions at 
doses of 4 X 10^= and 1 x 10^® Kr+lcm^, is depicted in Fig 4.2. It shows that 

the mixing progresses as a function of increasing dose and is indicated by a systematic 
decrease in the peak to trough ratios of the RBS signal. In order to determine the com- 
position of the mixed layer of irradiated samples, the experimental RBS spectra were 
simulated using structures of the kind 5f02/(7u/C'u3rGei_i/Ge/C'uir • • ■ fCu^Gei-xIGe 
as input to the simulation program. A typical result of such type of analysis for an 
irradiated sample is presented in Fig 4.3 for Cu/Ge multilayer sample irradiated at a 
dose of 8 X 10^® Kr'^/cm'^. The detailed layer structure used for simulation is given in 
Table 4.2 and the simulated spectrum is shown in Fig 4.3 by a continuous line. The 
complete structure can be visualized cis a layer of thickness 200 nm consisting of a film 
of composition GursGe 25 and unreacted Ge and Cu atoms. The simulated spectrum 
shown by continuous line in the Fig 4.3 is the the best visual fit that could be obtained. 
Mostly, the simulation analysis gives unique result and is significant as no two entirely 
different types of compositions and structures can give rise to same backscattering 
spectrum. 

The effect of Kr irradiation at a higher dose of 1 x 10^® Kr"^ f cm? is depicted by 
the experimental and theoretical RBS spectra in Fig 4.4. Here the wiggles in the ex- 
perimental backscattering spectrum in comparison to Fig 4.3 have almost disappeared 
indicating the formation of a homogeneous mixture of Cu and Ge. The peak at the 
surface position of Ge indicates that the surface Ge layer has not reacted completely. 
First, the simulation was tried using a single layer of thickness 190 nm having a com- 
position of CtirsGezs along with an overla\-er of 12 nm of Ge. This simulated profile is 
shown by the continuous fine in Fig 4.4. The fit appears reasonable except for the small 
hump around channel number 400. For fitting this hump a thin layer of composition 
Cu^Gc^r was introduced in the layer structure as shown in Table 4.3 which gives a 
better simulation fit shown by dashed line. 

It may be noted that the compositions of GursCejs and Cu«iGen used in simu- 
lation resemble the Cu^Gc and Cu^Ge phases shown in the Cu-Ge equilibrium phase 
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Figure 4.2: RBS spectra of Cu/Ge multilayer sample in the as-deposited condition and 
after ion-beam mixing at room temperature for two different doses. 
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Table 4.2: Description of the layer structure used for the simulation of RBS spectrum 
of ion mixed Cu/ Ge multilayer sample shown in Fig 4.3. 


Sample 

description 

Layer structure used for simulation 

Layer 

No. 

Thickness 

(nm) 


Composition 

Cu/Ge Multilayer 

1 

12 

Ge 

LOO 



films after ion 

2 

16 

Ge 

0.25 

Cu 

0.75 

mixing at RT for 

3 

7 

Cu 

1.00 



a dose of 

4 

13 

Ge 

0.25 

Cu 

0.75 

8 X 10^® Kr/cm? 

5 

4 

Ge 

1.00 




6 

17 

Ge 

0.25 

Cu 

0.75 


7 

2 

Cu 

1.00 




8 

20 

Ge 

0.25 

Cu 

0.75 


9 

1 

Ge 

1.00 




10 

15 

Ge 

0.25 

% 

Cu 

0.75 


11 

2 

Cu 

1.00 




12 

17 

Ge 

0.25 

Cu 

0.75 


13 

2 

Ge 

1.00 




14 

25 

Ge 

0.25 

Cu 

0.75 


15 

3 

Cu 

1.00 




16 

19 

Ge 

0.25 

Cu 

0.75 


17 

1 

Ge 

1.00 




18 

20 

Ge 

0.25 

Cu 

0.75 


19 

1000 

Si 

1.00 

0 

2.00 
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Figure 4.4: RBS spectum of ion-beam mixed Cu/Ge multlayer sample irradiated at 
RT for a dose of 1 x Kr+jcw? along with the simulated spectrum by considering 
formation of Cu^Ge and Cu^Ge phases. Another simulated spectrum obtained by 
considenng formation of only CuzGe layer is also shown for comparison. 
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Table 4.3: Description of the layer structure used for the simulation of RBS spectrum 
of ion mixed Cu/Ge multilayer sample shown in Fig 4.4. 


Layer structure used for simulation 


Sample 


Layer Thickness Composition 


description 

No. 

(nm) 





Cu/ Ge multilayer 

1 

12 

Ge 

1.00 



films after ion 

2 

20 

Ge 

0.25 

Cu 

0.75 

mixing at RT for 

3 

170 

Ge 

0.25 

Cu 

0.75 

a dose of 

3 

1500 

Si 

1.00 

0 

2.00 


1 X 10^® Kr/cm?' 
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diagram (see Fig 4.26) Therefore, one can predict from the RBS analysis that at room 
temperature, Kr ion irradiation up to dose of 8 x 10^® A'’r+/cm^ results in the forma- 
tion of a CuzGe phase and higher Kr doses result in the formation of both Cu^Ge and 
Cu^Ge phases. 

It is interesting to note that the ion-beam mixing carried out with Cu/ Ge multilayer 
samples at a temperature of 100 K ( near liquid nitrogen temperature) also shows the 
formation of a mixed layer having composition of GujiGe^z resembling the Cu^Ge 
phase. This is illustrated by the experimental and simulated RBS spectra presented 
in Fig 4.5. The irradiated area on this sample appears dark to naked eyes unlike the 
RT irradiations where the irradiated region appears reddish. The dark appearance of 
the irradiated region is due to the deposition of carbon film caused by the cracking 
of hydrocarbon on the surface. The carbon film thickness is negligible in case of RT 
irradiations while ion-beam mixing at 100 K, carbon deposition of as thick as 70 nm has 
been observed. This leads to the shift in the surface position of Ge as shown in the RBS 
spectrum of Fig 4.5. As mentioned in Chapter 3, the experimental chamber is equipped 
with a turbo molecular pump to get high vacuum which is almost free of hydrocarbons. 
Therefore there is almost negligible deposition of carbon at RT irradiation. It is felt 
that, when the sample is taken at 100 K, the small quantity of residual hydrocarbons 
which is still present in the chamber condenses on the sample at such low temperature 
and gets decomposed under irradiation fonning a carbon layer on the sample. 

For definite inference about the nature of phases formed, it is necessary to perform 
experiments which are sensitive to the atomic structure of the mixed region viz. X-rav 
diffraction (XRD), transmission electron microscopy (TEM) etc., the results of which 
can in turn be used to improve the quality of the simulation analysis by incorporating 
them m constructing the theoretical layer structures. The changes in surface topog- 
raphy, VIZ. formation of globules on the surface etc. can also affect the simulation 
analysis. The scanning electron microscopy (SEM) pictures taken for the as-deposited 
Gu/Ge multilayer sample and the one taken after irradiating the sample for a dose 
of 1 X Kr^cm^ are shown in Fig 4,6. Note that the surface topography of the 
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Figure 4.5: RBS spectum of ion mixed Cu/Ge multlayer sample at 100 K for a dose of 
8 X 10^® Kv^fcrn^ along with its simulation. 
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cLS-deposited sample is featureless aad it remains so after irradiation. 

The x-ray diffraction pattern for the as-deposited Cu/Ge multilayer structure sam- 
ple whose RBS spectrum is shown in Fig 4.1 is depicted in Fig 4.7a. The figure shows 
diffraction peaks for 26 values at 43.45° and 48.55° due to copper (111) and (200) planes 
respectively [84]. The Ge peaks are not seen in the entire diffractogram. This might 
be due to the amorphous nature of Ge layers which is caused by the fast quenching 
during the deposition process. The amorphous nature of Ge films prepared by vacuum 
deposition has been reported by Chopra et al.[85]. After irradiating the sample at RT 
by Kr ions for a dose of 8 x 10^® Kr^jcm?, the diffraction peaks belonging to copper 
disappear and the diffraction pattern (Fig 4.7b) exhibits an intense and sharp peak 
at 26 = 42.75° along with a weak peak at 26 = 45.47°. These peaks are found to 
correspond with the (020) and (012) planes of monoclinic Ci - Cu^Ge phase having 
a=0.2631 nm, b=0.42 nm, c=0.4568 nm and = 89.68° [47]. The observed diffraction 
peaks and the ones reported in the literature are compared in Table 4.4. 

Note that the RBS analysis of this sample presented in Fig 4.3 showed the formation 
of mixed layer of composition CuT 5 Ge 25 which is the composition of the compound 
CuzGe. Thus both RBS and XRD measurements confirm the synthesis of Cu^Ge 
phase at RT at this dose as well as for lower doses of 6 x 10^® Kr'^'fan^ and 4 x 10^® 
Kr+/cmK 

The sample irradiated at 100 K up to a dose of 8 x 10^® Kr+/cm^ does not show 
any clear peak (weak reflection) to confirm the formation of a compound phase. The 
XRD plot for the sample irradiated at room temperature for a higher dose of 1 x 
10^® Kr+lcm^ is shown in Fig 4.7c. Here, in addition to CusGe diffraction peaks, 
two distinct peaks at 2^ values of 40.40° and 43.30° also appear which belong to the 
hexagonal Cu^Ge phase also termed as ( phase having a=0.2582 nm and c/a = 0.1632 
nm[47]. The expected and observed diffraction lines are tabulated in Table 4.4. Thus 
the predictions of RBS analysis, that the Cu^Gt phase forms up to a dose of 8 x 10^® 

Kr+/cm^ and both Cu^Ge and CusGe phases form at higher doses are confirmed by 
the X-ray diffraction studies. 
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Figure 4.6: SEM photograph of Cu/Ge multilayer sample a) before and b) after ion- 
beam mixing at RT for a dose of 8 x 10^® Kr'^fcm^. 


76 



2 6 (degree) 


Figure 4-7: X-ray pattern of Cu/Ge multilayer sample taken before and after ion-beam 
mixing at RT (a) as-deposited (b) 8 x 10^® Kr^Jcrri^ (c) 1 x 10^® AV+/cm2. 
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Table 4.4: Comparison of experimentally obtained d-values of the observed copper 
germanide phases with the one reported in literature. 


Cu/Ge multilayer 
sample 

^Exp. ^Rep. 

(A) (A) 

(hkl) Phase 

As-Deposited 

2.08 2.08 

1.87 1.808 

(111) Cu 
(200) 

Ion mixed at RT, 

8 X 10^5 Krlcrr^ 

2.11 2.10 

1.99 2.01 

(020) CusGe 
(020) 

Ion mixed at RT, 

1 X 10^® Krlcm^ 

2.23 2.24 

2.11 2.10 

2.08 2.12 

1.99 2.01 

(100) CusGe 
(020) CusGe 
(002) CusGe 
(020) CusGe 
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The effect of ion-beam mixing in the Cu/Ge multilayers using 1 MeV Ar ions for 
a dose of 1 x 10^® Ar'^ jcrri^ is compared with the as-deposited sample as shown by 
the RBS spectra in Fig 4.8. It is clear that complete mixing has not occurred in this 
case. The simulation analysis shows the formation of a. mixed layers of composition 
C'u 7 sGe 25 across each Cu/Ge interface indicating synthesis of a CuzGe compound. The 
XRD analysis also confirms the formation of Cu^Ge compound. Although the mix i ng 
effects observed using Ar ions are qualitatively similar to the those seen in case of 
Kr irradiation, the rate of mixing is lower as compared to the Kr ion irradiation for 
the same dose and incident energy. The amount of mixing obtained with Ar ions at 
a dose of 1 x 10^® ions fern? is almost same as that obtained with Kr ions at a dose 
of 4 X 10^® ionsfern^. The formation of these copper germanide phases has also been 
observed in bilayer samples. The detailed comparison of mixing rates and efficiency 
values obtained using both these ions is dealt with in Chapter 6. The results of RBS 
and XRD analysis of these ion mixed multilayer samples are summarized in Table 4.5. 
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Figure 4.S: RBS spectra of Cu/Ge multilayer sample before and after ion-beam mixing 
at RT for a dose of 1 x 10*® Ar’^/arri^. The simulated curves are shown by dased and 
continuous line respectively. 
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Table 4.5: Summary of ion-beam mixing results in Cu/Ge multilayer systems using Kr 
and Ar ions at various doses and temperatures. 


Sample 

Ion 

Temperature 

Dose 

Phase observed after IM 

description 

species 

of Irrad. (K) 

ions j cm? 

composition 

Phase 

Cu/Ge 

Kr 

300 

(2 - 8) X 10^5 

C iiTsGezz 

C u^Gc 

multilayer 






sample 



(1 - 2) X 10^® 

C urzGtzz 

CuzGe 

deposited 




+ 

+ 

on SiOz 




C u^Gti^r 

CusGe 


Kr 

100 

(4 - 8) X 10^5 

C U7sGe25 

\ 

* 


Ar 

300 

(4 - 10) X lO^s 

C ujsGezs 

Cu^Gc 


*Weak X-ray reflection of CuzGe phase 
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Thermal annealing effects 

The thermal annealing experiments have been performed at 630 K for different an- 
necding durations from 30-180 min. After annealing, the peaks and troughs in the 
as-deposited spectrum start reducing just as in the case of ion mixed samples. Fig- 
ure 4.9 shows RBS spectra of a multilayer sample before and after annealing for 150 
min. which shows the formation of a homogeneous layer of Cu 75 Gt 25 . The RBS anal- 
ysis shows that the composition of the reacted region is Cu 7 ^Ge 2 z and is similar to 
the one obtained in ion-beam mixed samples. The simulated structures used for an- 
nealed sample is shown in Table 4.6. It is seen from the table that a single layer of 
thickness 190 nm of composition Cu 7 zGt 2 z was used- for the simulation of the reacted 
layer. Similar results were also found for other annealed samples. The X-ray studies of 
this multilayer samples annealed at 630 K (Fig 4.10) show diffraction peaks of CuzGe 
phase while it does not show any peak corresponding to the copper rich CusGe phase. 

Figure 4.11 shows a RBS spectrum of a bilayer sample before and after aimealing 
at 580 K for 60 min. The simulation analysis gives the stoichiometric composition 
of Cu 7 zGe 2 z- The X-ray studies also show diffraction peaks from Cu^Ge crystalline 
compound. The results of thermal annealing of multilayer and bilayer samples are 
summarized in Table 4.7. 
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Figure 4.9: RBS spectra of Cu/ Ge multilayer sample before and after thermal annealing 
at 630 K for 150 min along with their simulations. Annealing leads to the formation 
of homogeneous Cu^Ge phase. 
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Figure 4.10: X-ray pattern of Cu/Ge multilayer sample before and after thermal 
nealing (a) as-deposited (b) 630 K for 150 min. 
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Figure 4.11: RBS spectra of Cu/Ge bilayer sample 
K for 60 min along with simulations. 
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Table 4.6: Description of the layer structure used for the simulation of RBS spectrum 
of thermally annealed Cu/Ge multilayer sample shown in Fig 4.9. 


Sample 

description 

Layer structure used for simulation 

Layer Thickness Composition 

No. (iirn) 

Cu/Ge multilayer 
samples after 

annealing at 

630 K for 150 min. 

1 190 Ge 0.25 Cu 0.75 

2 1500 Si 1.00 0 2.00 

i 


4.3.2 Ni/Ge system 
Ion-beam induced effect 

For Ni/Ge system, the Kr and Ar ion-beam mixing experiments were performed at 
sample temperatures ranging from 100 K to 425 K. The RBS spectrum of as-deposited 
multilayer sample (open circle) is shown in Fig 4.12 along with the simulated spectrum. 
The sample consists of five pairs of nickel and germanium layers with an overall thick- 
ness of Ni and Ge equal to 170 nm and 120 nm respectively. The detailed simulated 
structure is shown in Table 4.8. 

The RBS spectrum after irradiation with Kr ions (triangle) for a dose of 6 x 10^® 
Kr'^fcm^ along with the simulated spectrum (continuous line) are also shown in 
Fig 4.12. This spectrum could be simulated using the structure shown in Table 4.8 
consisting of a single layer of thickness 240 nm and composition for the mixed 

region. This composition corresponds to Ni^Gt phase present in the Ni-Ge equilibrium 
phase diagram (shown in Fig 4,27). Irradiations up to a dose of 1 x 10'® Ar'*‘/cm* do 







Table 4.7: Summary of thermal annealing results of Cu/Ge sample 


Sample 

Geometry 

Summary of thermal annealing results in Cu/Ge sample 

Temperature Annealing Composition Phase identified 
in K time in min. given by RBS by X-ray 

Cu/Ge 

multilayer 

sample 

deposited 

on Si 02 

630 30 01*750625 Cu^Gc 

90 

120 

150 

180 

Cu/Ge Bilayer 
film deposited 

on 5*02 

580 30 Cur$Ge2s Cu^Ge 

60 




Normalized yield 


Energy (MeV) 
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Channel 


Figure 4.12: RBS spectra of the as-deposited and ion mixed (6 x 10^® Krfcm^, RT) 
Ni/Ge multilayer sample deposited on Si02’ The Continuous lines show simulations. 
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Table 4.8: Description of the layer structure used for the simulation of RBS spectra of 
as-deposited and ion mixed Ni/Ge multilayer sample shown in Fig 4.12. 


Layer structture used for simulation 

Sample Layer Thickness Composition 

description No. (nm) 

Ni/Ge multilayer 1 26 Ge 1 

films deposited 2 32 Ni 1 

on Si02 3 24 Ge 1 

4 34 Ni 1 

5 22 Ge 1 

6 34 Ni 1 

7 22 Ge 1 

8 35 Ni 1 

9 24 Ge 1 

i 10 35 Ni 1 

11 1000 Si 1 0 2 

Ni/Ge multilayer 1 10 Ge LOO 

films after ion 2 240 Ge 0.34 Ni 0.66 

mixing at RT for 3 15 Ni 1.00 

a dose of 4 1000 Si 1.00 0 2.00 

6 X 10^5 AV/cm2 
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not change the composition of the mixed regions. For irradiation performed at higher 
sample temperature of 390 K the mixing in Ni/Ge system proceeds rapidly with the 
formation of a uniform composition corresponding to the Ni 2 Ge phase. This is illus- 
trated by the experimental and simulated RBS spectra given in Fig 4.13 corresponding 
to a dose of dose of 4 x 10^® Krfcrn^. Ion-beam mixing at a sample temperature of 
100 K the mixed region leads to the same composition across the individual Ni/Ge 
interface. The temperature dependence of mixing and ion induced phase formation 
have been discussed in detail in Chapter 6. 

Conventional X-ray diffraction measurements as well as those performed in Seeman- 
Bohlin geometry did not show any peak for the as deposited Ni/Ge multilayer sample 
indicating the amorphous/highly disordered state of the individual Ni and Ge layers[85, 
86]. The Seeman-Bohlin X-ray diffraction spectra recorded after irradiating the sample 
at RT with a dose of 6 x 10^® Kr'^fcm^ is shown in Fig 4.14a, which shows diffraction 
peaks at 20 values 45.14° and 46.68° belonging to Ni^Ge phase. The observed and 
reported[84] d values for phase are given in Table 4.9. The irradiations at 

other doses also show the formation of Ni^Ge phase. No other Ni-Ge phase has been 
observed after ion-beam mixing. 

Irradiations with Ar ions also lead to the formation of Ni 2 Ge phase. The X- 
ray diffraction spectrum taken after irradiating the Ni/Ge multilayer sample at RT 
with Ar ions for a dose of 1 x 10^® ionslcm^ is shown in Fig 4.14b which also shows 
diffraction peaks at 20 = 47.6° and 49.68° belonging to Ni 2 Ge phase[84]. However, the 
sample irradiated at 100 K dose not show any prominent peak although the composition 
corresponding to ATiaGe phase has been obtained form RBS analysis. The experimental 
and reported d- values are also compared in Table 4.9. It should be noted that the effect 
produced due to RT irradiation for 1 x 10^® ions ( cm? of Ar is obtained at a much lower 
dose of 4 X 10^® ions fern? of Kr ions for RT irradiations and is clearly illustrated in 
Fig 4.15. The results of ion-beam mixing presented here are summarized in Table 4.10. 
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Figure 4.13: RBS spectra of Ni/Ge multilayer sample before and after ion-beam mixing 
at 390 K for a dose of 4 x 10^® Kr'^fcm? showing formation of a uniformly mixed region 
of composition Nie^Gesi. 
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Figure 4.14: X-ray pattern of Ni/Ge multilayer sample after ion-beam mixing at for 
irradiation fluence of (a) 6 x 10^® AV'^/cm^ (b) 1 x 10^® Ar'^'/cm^. 
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Table 4.9; Comparison of experimentally obtained d-values of the observed Ni 2 Gt 
phase with the one reported in literature. 



Experimental 


Reported 


Sample 

Dose 

dhkl 

dhkl 

(hkl) 

Phase 

description 

{ions ! cm?) 

(A) 

(A) 



X-ray spectrum of 

6 X 10^® Kr[cm? 

2.01 

2.04 

(211) 

NiaGe 

Ni/Ge multilayer 
sample after RT 


1.94 

1.91 

(020) , 


ion-beam mixing 

1 X 10^® At fen? 

1.91 

1.91 

(020) 

NhGe 



1.834 

1.836 

(212) 
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Figure 4.15: Comparison between the RBS spectra of Ni/Ge multilayer sample after 
ion-beam mixing at doses of 4 x 10^® A'r+/cm? and 1 x 10^® Ar'^fcm^. 
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Table 4.10: Summary of ion-beam mixing results in Ni/Ge multilayer systems using 
Kr and Ar ions at various doses and temperatures. 


Sample 

description 

Ion 

species 

Temperature 

during Irrad. 

Dose 

ions/cm^ 

Phase obtained by IM 

composition Phase 

Ni/Ge multilayer 

Kr 

300 

(2 - 10) X 10^® 

NiseGes4 

Ni^Ge 

film deposited 

on Si 02 

Kr 

100 

(6 - 10) X 10^5 

7 ? 



Kr 

390 

(2^ 10) X 10^5 

77 

NhGe 


Kr 

430 

1 X 10^® j 

77 

77 


Ar 

300 

(5 - 10) X 10^5 


NhGe 


At 

380 

1 X 10^® 

77 

77 


* Very low intensity x-ray peak. 
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Thermal annealing effects 

The RBS spectrum for a multilayer Ni/Ge sample annealed at 580 K for 60 min. is 
shown in Fig 4.16 along with the corresponding as-deposited spectrum. As seen in 
the spectrum a uniform homogeneous composition has been obtained after annealing. 
The spectrum of annealed sample was simulated using a composition of and 

some unreacted Ni and Ge. Similar results were obtained for other temperatures and 
annealing times. Figs 4.17a and 4.17b show the standard X-ray spectra for samples 
annealed at 580 K for 60 min. and 180 minutes. The analysis shows that the peaks 
correspond to the stable NizGe phase. The formation of Ni^Ge phase was detected 
after annealing at temperature as low as 480 K by normal and glancing angle X-ray 
techniques. 

Bilayer samples also show synthesis of the nickel germanide phases after annealing 
at various temperature. A typical RBS spectrum of an as-deposited bilayer sample 
is shown in Fig 4.18a. The simulated spectrum is shown by continuous line and the 
simulated structure is given in Table 4.11. The RBS spectrum after annealing at 480 
K for 180 min. is also shown in Fig 4.18a along with its simulated spectrum. The sim- 
ulation analysis shows the formation of two layers having compositions corresponding 
to NiGe and Ni^Ge phases. The simulated structure is shown in Table 4.11. The XRD 
spectrum for this annealed sample is shown in Fig 4.18b which confirms the formation 
of NiGe and Ni^Gt phases. Figure 4.19 shows the comparison of RBS spectra of three 
bilayer samples after annealing at 430 K for 90 min., 480 K for 60 min. and 530 K for 
180 min. respectively. The detailed analysis show that only Ni 2 Ge phase is produced 
at 430 K and both NiGe BudNiiGe are produced above 475 K. These results of RBS 
and XRD analysis of the thermally annealed multilayer and bilayer Ni/ Ge samples are 
summarized in Table 4.12. 
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Figure 4.16: RBS spectra of Ni/Ge muItOayer sample before and after annealing at 
580 K for 60 min. along with simulations. 



Figure 4.17: X-ray pattern of Ni/Ge 
K for (a) 60 min. (b) 180 min. 





98 


Energy (MeV ) 




Figure 4.18: (a)RBS spectra of Ni/Ge bilayer sample before and after annealing at 
480 K for 180 min. along with simulations (b) X-ray diffraction pattern confirms the 
formation of iVijGe and NiGe phases. 
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Table 4.11; Description of the layer structures used for the simulations of RBS spectra 
of as-deposited and thermally annealed bilayer Ni/Ge bilayer sample shown in Fig 4.18. 


Sample 

description 

Layer structure used for simulation 

Layer Thickness Composition 

No. (nm) 

Ni/Ge bilayer 
films deposited 

on Si 02 

1 106 Ge 1.00 

2 71 Ni 1.00 

3 1000 Si LOO 0 2.00 

After annealing 

at 480 K for 180 min 

1 2 Ge 1.00 

2 80 Ge 0.50 Ni 0.50 

3 30 Ge 0.34 Ni 0.66 

4 50 Ge 0.50 Ni 0.50 

5 1000 Si 1.00 0 2.00 
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Figure 4.19: Compaxison of RBS spectra of Ni/Ge bilayer samples annealed at three 
different temperatures (a) 430 K for 90 min. (b) 480 K for 60 min. (c) 530 K for 180 


mm. 
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Table 4.12: Summary of thermal annealing results of Ni/Ge sample 


Sample 

Annealing 

Annealing 

Phase obtained by annealing 

description 

temperature 

in K 

time 

(in min.) 

Composition 

Phase 

Ni/Ge multilayer 

480 

60-180 

NiesG6s4 

NhGe 

films deposited 

on SiO^ 

530 

180 

Nie6Ge34 

Ni2Gc 


580 

60-180 


NizGe 

Ni/Ge bilayer 

430 

60-180 

NieeGe^ 

NhGe 

film deposited 

on Si 02 

480 

60-180 

Nio&Gez^+NisoGe^o 

NhGe+KiGe 


530 

i 

90-180 

NimGez^+NisoGesQ 

NiiGe+KiGe 
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4.3.3 Co/Ge system 

Ion-beam induced effect 

The multilayer Co/Ge samples were made with four pairs of cobalt and germanium 
layers each having total thickness of 110 nm. The RBS spectra in Fig 4.20 compare the 
backscattering yields of as-deposited sample with the one irradiated by 1 MeV Kr'^ 
ions at a dose of 1 x 10^® ions j cm?. The analysis of the irradiated sample shows the 
formation of mixed regions having composition of at the interface of each 

layer along with unreacted layers of Co and Ge. The simulated curves are shown by 
continuous line for both the cases. It may be noted that unlike Cu/Ge and Ni/Ge 
system the multilayers of Co/Ge do not mix completely for the irradiation dose of 
1 X 10^® Kr'^ (cm?. The SEM studies of as-deposited and irradiated sample do not 
show any noticeable change in surface features. 

The effect of ion-beam mixing performed at elevated temperatures of .380 K and 420 
K are depicted in Fig 4.21 and Fig 4.22 respectively for an irradiation dose of 1 x 10^® 
Kr'^ (cm? along with as-deposited spectrum. A single layer of composition 
was used for simulation of the mixed region in both the cases. The layer structures used 
for simulation for the irradiated samples as shown in Fig 4.21 is given in Table 4.13. 
It is also clear from these figures that at these elevated temperatures the mixing is 
more uniform across the thickness, nevertheless, a small amount of topside Ge film and 
the bottom side Co film stiU remain unreacted. The X-ray diffraction patterns of the 
irradiated sample shows many new diffraction peaks as compared to the as-deposited 
spectrum. Unfortunately, x-ray data for C 02 G& is not available hence we could not 
identify x-ray diffraction peaks corresponding to this phase. However, the very fact 
that RBS gives uniform stoichiometry and that x-ray pattern indicate crystallinity of 
the film hence we conclude that the composition of the film belongs to the 

C 02 GC. phase and that the phase has formed after ion-beam mixing. The RBS results 
of ion mixing at different temperatures are summarized in Table 4.14. 
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Figure 4.20: RBS spectra of Co/Ge multilayer sample before and after ion-beam mixing 
at RT for a dose of 1 x 10*® Kr'^fcm^. The simulated curves in both cases axe shown 
by continuous line. 
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Figure 4.21: RBS spectra of Co/Ge multilayer sample before and after ion-beam mixing 
for a dose of 1 x 10^® A r"*"/ cw^ at 380 K. In this case mixing leads to the formation of 
CoiGe phases. The simulated curve is shown by continuous line. 
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Figure 4.22: RBS spectra of Co/Ge multilayer sample before aad after ion-beam mixing 
for a dose of 1 x 10^® Kr'^[an^ at 420 K. In this case mixing also leads to the formation 
of Co^Ge phases. 
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Table 4.13: Description of the layer structures used for the simulation of RBS spectrum 
of ion mixed Co/Ge multilayerr sample shown in Fig 4.21. 


Layer structure used for simulation 


Sample 

description 


Layer Thickness Composition 
No. (nm) 


Co/Ge multilayer 

1 

7 

Ge 

1.00 



sample after 

2 

130 

Ge 

0.34 

Co 

0.66 

ion-beam mixing at 

3 

10 

Co 

1.00 



380 K for a dose 

4 

1000 

Si 

1.00 

0 

2.00 


of 1 X 10^® Krjcrri^ 




Table 4.14: Summary of ion-beam mixing results of Co/Ge sample 


Sample 

Ions 

Temperature 

Dose 

Phase obtained by IBM 

description 


of Irrad (K). 

ions/cm^ 

Composition 

Phase 

Co/Ge 

Kr 

300 

(6-10) X 10^5 

CoseGc^^ 

CO 2 G 6 * 

multilayer 






sample 

Kr 

380 

1 X IQi® 

C 

Co^Ge 

deposited 






on SiOz 

Kr 

430 

1 X 10‘® 

CoeeGcz4 

Co^Ge- 


Ar 

300 

1 X 10‘® 

C 066^634 

CozGf 


Ar 

380 

1 X 10‘® 

CoeeGeu 

Co^Ge* 


* X-ray data for CozGe phase are not available 
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Thermal annealing effects 

In case of Co/Ge the thermal annealing experiments have been done for both bilayer 
and multilayer samples at 480 K to 680 K for different annealing durations. The RBS 
spectrum for a multilayer sample annealed at 530 K for 180 min. is shown in the 
Fig 4.23a along with the as-deposited spectrum. The spectrum of annealed sample 
can be simulated (shown by continuous line) by a layer of composition of Co-iGe and 
some unreacted Co and Ge. The multilayer sample whose RBS spectrum is shown in 
Fig 4.23b completely reacts after annealing at 680 K for 120 min. duration and the 
simulation analysis shows the formation of uniform layer of composition Co^Ge. The 
XRD analysis also shows a single sharp peak. The RBS spectra of Fig 4.24 compare the 
ion-beam mixing effect with the one induced by thermal annealing alone . It is seen 
form this figure that the sample annealed at 580 K for 60 min. matches quite well with 
the spectrum of the sample irradiated at 380 K for a dose of 6 x 10^® Kr'^/cm?. 

In case of bilayer samples, annealing at 480 K produces composition corresponding 
to Co^Ge phase. Whereas at 530 K both Co^Ge and CoGe phases are formed at Co rich 
and Ge rich sides respectively. At still higher temperature of 580 K both of the layers 
react almost completely forming a uniform composition of CoGe. Figure 4.25 shows 
the comparison of RBS spectra corresponding to a bilayer sample annealed at three 
different temperatures. The detailed results of these bilayer and multilayer samples are 
tabulated in Table 4.15. 


4.3.4 Resistivity measurements 

The resistivity of the films after ion-beam mixing was measured by four probe method 
as described in Chapter 3. The measured resistances are shown in Table 4.16 for. both 
ion mixed and thermally produced phases. In case of ion mixed phases, the resistivity 
values are reported for the minimum ion dose at which a uniform homogeneous layer 
is formed. The resistivity of the phases produced by ion-beam mixing are almost twice 
as compared to the phases produced by thermal annealing. This might be due to the 
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Figure 4.23: 
for 180. min. 


sp«t« of thermally annealed Co/Ge multilayer samples (a) 530 K 
fian K for 120 min. Simulated spectra are given y coo muou 


lines. 
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Energy (MeV) 



Figure 4.24: Comparison between RBS spectra of Co/Ge multilayer sample taken after 
ion-beam mixing at 380 K for a dose of 6 x 10^^ Kr'^fcm^ and after thermal annealing 
at 580 K for 60 min. 
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:e 4.25: Company, ot RBS sp«:tra of thermally annealed Co/Ge bilarer 

ree different temperaturee (a) 480 K for ISO min. (b) 530 K for 120 mm. (c) 580 


K for 120 min. 
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Table 4.15; Summary of thermal annealing results of Co/Ge sample 


Sample 

Annealing 

Annealing 

Phase obtained by annealing 

description 

temperature 

time 

Composition 

Phase 


inK 

in min* 



Co/Ge multilayer 

530 

60-180 

C oeeG^34 

C 02 G 6 * 

sample deposited 





on Si 02 

580 

60-120 

C 066^634 

Co2Ge* 


680 

120 

C 0seGe^4 

Co2Ge* 

Co/Ge bilayer 

480 

90-240 

CoseGtz4 

Co2Ge* 

sample deposited 





on Si 02 

530 

60-240 

C o^sGez4 

Co2Ge* 




+ 

+ 




Co^oGcso 

CoGe 


580 

60-120 

C 

CoGe 


* see text (p 102) 
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Table 4.16: Resistivities measured on Cu^Ge^ Ni^Gc and Co^Ge layers produced by 
ion-beam mixing and thermal annealing. 



Phase synthesized by 
ion mixing 

Phase prepared by 
thermal annealing 

Binary 

system 

Dose 

Krfcm? 

resistivity 

fiO. — cm 

Temperature and 
annealing time 

Resistivity 
yiQ, — cm 

Cu/Ge 

8 X 10^® 

20 

350 “G, 150 min. 

8.6 

Ni/Ge 

8 X 10^® 

50 

300 '’G, 180 min. 

23 

Co/Ge 

1 X 10*® 

74 

300 “G, 120 min. 

50 


presence of irradiation induced defects in the mixed regions. The resistivity of the 
CusGe phases compare reasonably well with the reported value of 6 /z-cm[49]. For the 
Ni-iGt and CotGe phases no resistivity data is available in the literature. Probably 
these are the first measurements to be reported. It is important to note that the 
resistivity values are lower as compared to metal silicides. Hence we believe that these 
germanide phases would become important for semiconductor industry. 
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4.4 Discussion 

Before starting the discussion let us summarize the main points of the ion-beam mixing 
and thermal annealing results of Cu/Ge, Ni/Ge and Co/Ge systems presented in Tables 
4.5, 4.7, 4.10, 4.12, 4.14 and 4.15 respectively. 

• Both Kr and Ar ion irradiation induces mixing in the multilayer metal/Ge films 
at 100 K, RT and above. In all the cases the metal-rich compounds viz Cu^Ge, 
Ni^Ge and Co^Ge form first in Cu/Ge, Ni/Ge and Co/Ge systems respectively 
after mixing at RT and above. Thermal annealing alone also produces the same 
phases which suggests that synthesis of the first phase is same for both ion-beam 
mixing and thermal reactions. 

• At higher doses the ion-beam mixing also produces Cu^Ge phase. In case of 
Ni/Ge and Co/Ge bilayer systems, higher annealing temperatures (above 475 K) 
produce NiGe and CoGe phases respectively. 

• .A.r ion induced mixing occurs at a much slower rate as compared to Kr irradiation. 


4.4.1 Phase formaton mechanism 

In general, ion-beam mixing refers to the homogenization across the interfaces of mul- 
tilayers or bilayers of two elements by energetic heavy ion beams. When the energetic 
ion passes through these layers (as described in Chapter 2) it looses most of its energy 
by soft atomic collisions and comes to rest within 10"^® s[14, 17]. During primary 
collisions the atoms of the target acquire sufficient energy and gets dislodgedfrom their 
original positions producing a number of secondary collisions. Thus within a short 
period of time, all the atoms in a small volume are set into random motion, which is 
commonly known as collisional cascade which lasts for about 10“^^ to 10“^^ s. During 
this short time period the displaced atoms interact in a collective fashion producing 
a significant amount of disorder. At the end of the collision cascade when the veloc- 
ities of the moving atoms assume Boltzmann distribution one can invoke the notion 
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of temperature. Subsequently a relaxation process starts in which the excited mixture 
transits to a state of lower free energ}^ This process lasts for another 10“^^ to 10-^° 
s and is called as thermal spike[14]. If we consider the temperature of the spike to be 
of the order of 1000 K the region cools down with rates of as high as 10^^ K/s. It has 
been commonly accepted that the ion mixed alloys form during the relaxation process 
and not in the collisional process because the atoms are always in dynamically excited 
states in the cascade. In principle, if the temperature and time conditions during re- 
laxation are sufficient for crystallization, the random mixture transforms to a stable or 
metastable crystalline state[13, 15]. However, if the conditions are not favorable then 
amorphous phase is formed. 

ThU general picture described above has been found to be true for most of the 
metal/metal systems which form crystalline phases at high irradiation temperature 
and amorphous phases around room temperature and below. The characteristics of 
these binary metal/metal systems are (i) the non-reactive nature of the elements and 
(ii) their binary phase diagrams show equilibrium compounds which can be formed by 
thermal treatment at very high temperatures. On the other hand, many highly reactive 
systems like metal/semiconductor[15, 17, 30, 31, 41, 42, 57] and metal/Al[25, 58, 87] 
form simple as well as complex compounds around room temperature or below. It is 
interesting to note that these systems form compound by thermal treatment around 
400-500 K. It should be mentioned over here that the thermal annealing at elevated 
temperatures also form the same phase as that formed by ion-beam mixing suggesting 
a common mechanism for the formation of the first phase in both the cases[57]. 

The emphasis here is to put the mixing processes on the scale of ’phase structure’. In 
this sense it is relevant to treat the relaxation process by invoking thermodynamics and 
kinetics of crystallization[13]. Apart from temperature, there are many other factors 
which influence the final phase formation. In general, the factors which are most likely 
to influence the formation of phase and its nature are the structure of the terminal 
solution[13], complexity of the compound phase[13, 56], width of phase fields[59, 60], 
atomic size[55, 62], atomic mobility[27. 56], heat of formation of various phases[56, 62], 



116 


the effect of substrate temperature[13] etc. All these factors act together in a complex 
manner. It has been mentioned earlier in Chapter 2 that several emperical rules[.52-56] 
have been proposed from the knowledge of the final phase formed in a large number 
of binary systems. In turn, these rules have been used for predicting the formation of 
a phase and its nature for an unknown system. The literature shows that these rules 
are very successful (within their limitations) in most of the metal/metal and metal/Si 
systems[6, 13]. 

Since the chemical nature of Si and Ge is quite similar, one expects a similarity 
of results in producing silicide and germanide phases by ion-beam mixing as well as 
by thermal annealing[48]. This is indeed true and can be seen from Table 4.17 which 
compares the results of ion-beam mixing in Metal/Si with the results of this work on 
metal/Ge systems. We now test the predictions of the existing rules regarding phase 
formation with the finding of the present work. 


4.4.2 Nature of Phase formation 

Here we investigate the rules which predict the possible nature of the phases formed 
under different conditions by ion-beam mixing. Table 4.18 gives a summary of some of 
the widely used rules and their predictions in case of metal/Ge systems studied here. 

The simple structure rule (column 1 of Table 4.18) proposed by Hung et al.[52] states 
that ion-beam mixing would produce crystalline compound having simple structure 
whenever the overall composition is close to that of a compound. Here simple structure 
means the fee or bcc or hep structures containing one or two atoms per unit cell. 
On the other hand, when the overall film composition is not close to an equilibrium 
alloy an amorphous phase would form. In the present case we have used samples 
having composition of one of the stable compounds present in the phase diagram. The 
equilibrium phase diagrams[91] of Cu/ Ge , Ni/Ge and Co/Ge binary systems axe shown 
in Figs 4.26, 4.27 and 4.28 respectively. The equilibrium phases present in the phase 
diagrams and other relevant information[63, 91-93] for these three systems are given 
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Table 4.17. Comparison of first phase formation in metal silicides and germanides 
by ion-beam mixing and thermal annealing. 


Metal / Semiconductor 

system 

Phase formed by IBM 

Phase formed by 
thermal annealing 

Ref. 

composition 

phase 

Cu/Ge 

CuzGe 

Cu^Gc 

CuzGe 

[this work] 

Cu/Si 

# 

# 

CuzSi 

[48] 

Ni/Ge 

Ni2Ge 

NkGe 

Ni2Ge 

[this work] 

Ni/Si 

Ni^Si 

Ni2Si 

NhSi 

[88] 

Co/Ge 

Co^Gt 

C 02 GC* 

Co2Ge 

[this work] 

Co/Si 

Co^Si 

C02Si 

Co2Si 

[89] 

Pd/Ge 

PdiGe 

Pd2Ge 

Pd2Ge 

[42] 

Pd/Si 

PdiSi 

Pd2Si 

Pd2Si 

[30, 90] 


# Not reported in the literature 


* X-ray data are not available. 
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Table 4.18; The nature of phases predicted by the existing empirical rules for the 
present case of metal germanide systems. 


Emperical 

Parameters 

Predictions for present systems 

Rule 

used 



Co/Ge 

(1) Simple 

Structure 

rule [52] 

Composition 

C 

c 

C 

(2) Hagg rule 
[55] 

atomic radii 

the elements 

C 

c 

C 

(3) Structure 

difference 

rule[53] 

Structure 

c 

c 

A 

(4) E.xtended 

structure 

ruIe[o4] 

Composition 

. k 

Phase diagram 

c 

! 

c 

C 


G = Crystalline phase 
A = Amorphous phase 
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Table 4.19. This table shows that the Cu/Ge system has two stable metal rich phases 
viz., CuzGt and CusGe. For this system we have used multilayer samples having 
composition close to the composition of CuzGt phase. Simple structure rule of Hung's 
et al. that the CuzGg phase should form, which is in accordance with our experimental 
results. The Ni/Ge system contains three stable phases NizGe, NiGe and NizGe in 
equilibrium phase diagram. In this case we have used a composition close to the Ni 2 Ge 
phase. This rule then predicts the formation of Ni^Ge compound which also agrees 
with our experimental results. Similarly, Co-Ge phase diagram has three stable phases 
namely Co 2 Ge, CoGe and CoGe 2 . Amongst these, we have chosen a composition of 
CoeeGeu corresponding to C 02 Ge phase. The prediction of this rule to form Co 2 Ge 
phase matches with our experimental findings. Table 4.19 shows that structure[92] of 
the CuzGe phase is CuzTi type whereas structures of the iVi^Ge and Co 2 Ge are that of 
Co 2 Si type. There are many reports of formation of this kind of structures by ion-beam 
mixing. For example, Pd 2 Si[o 6 , 94], €0281 which have complex structures, are formed 
by ion-beam mixing at room temperature. Similarly, in case of aluminide systems the 
formation of complex phases like AuA/'ifSS, 95], NizAl[25] at room temperature and 
below are reported in literature. 

Hagg[55] and Rauschenbach[96, 97] have proposed (Column 2 of Table 4.18) that 
complex crystalline structures can be formed in a binary system if the ratio (R) of the 
radii of metalloid and metal atoms happens to be greater than 0.88 and compounds 
having simple crystalline structure would form if R is less than 0.59. Between these 
two limits amorphous phase formation taJkes place. This is shown in Fig 4.29. It can • 
be seen from Table 4.19 and Fig 4.29 that in case of metal/Ge systems the ratio of 
radii[93] for all the three systems are greater than 0.88 thus suggesting the formation of 
compound phases in all the three cases. Thus both simple structure and radii ratio rule 
predict correctly the formation of compound at room temperature and above by ion- 
beam mixing. However, both these rules do not give any prerequisite of composition 
range under which these compound phases will be produced. 



120 


Atomic percent copper 


0 20 40 60 80 100 



Figure 4.26: Equilibrium phase diagram of Cu-Ge system. 
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Figure 4.27; Equilibrium phase diagram of Ni-Ge system. 




122 


Weight percentage germanium 



Figure 4.28: Equilibrium phase diagram of Co-Ge system. 




Table 4.19: The difFernet useful parameters of Cu/Ge, Ni/Ge and Co/Ge systems. 


Binary 

system 

A/B 

Lattice 

structure 

[92] 

atomic radii 

ratio R=rB/rA 
[93] 

Equilibrium 

phase 

[91] 

Width of 

phase field 
at%[91] 

kJ/mole 

[62] 

Crystal 

structure 

[92] 

Cu/Gc 

fee- fee 

0.95 

CuzGe 

4 

16 

CuzT i 




Cu^Gc 

8 

> 16 

Mg 

Ni/Ge 

fee- fee 

0.98 

iVijGc 

0 

37 

Co2Si 




NiGe 

0 

32 

MnP 




NizGe 

2.5 

25 

CuzAu 

Co/Ge 

hep-fee 

0.97 

C<hGt 

12 

15 

Co^Si 




CoGe 

3 

16 

FeSi 




CoGez 

3 

12 

PdSnz 
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The formation of phase can be also predicted by structure difference rule[53] of 
Liu (column 3 of Table 4.18 ). This rule states that an amorphous binary alloy will 
be formed by ion-beam mixing of multilayer samples when the two constituents are 
of different crystalline structures. Thus according to this rule in case of Co/Ge (hcp- 
fcc) one expects an amorphous phase but for Cu/Ge (fcc-fcc) and Ni/Ge (fcc-fcc) one 
expects a compound. Later Liu[54] proposed an extended structure difference rule 
bs-sed on the ecjiiilibrium phase diagram of the binary system under consideration 
(column 4 of Table 4.18). 

The extended structure rule states that (1) if the overall composition is in the two 
phase region of the phase diagram, an amorphous phase is most likely to form, (2) 
if the overall composition is in or near the single phase region of the phase diagram 
and structure of this phase is not simple, an amorphous phase is likely to be formed 
and finally, (3) if the overall composition is in or near the single phase region of the 
phase diagram and the structure of the phase is simple, a crystalline phase is formed. 
From the phase diagram (Fig 4.26-Fig 4.28) of these systems it is evident that the 
composition of all multilayer samples of Cu/Ge, Ni/Ge and Co/Ge are within the 
stable single phase region therefore according to third rule a compound phase should 
be formed in each case. 

Brimhall et al.[59] as well as others[60, 61] have pointed out that for better predic- 
tions, atomic mobilities and widths of the phase field width fAJC) have to be taken into 
consideration. It has been found that the situation for a phase having large phase field 
is different than a line compound in the phase diagram. For example, Ni/Al system[61] 
has three phases in the phase diagram NiAi {AX = 50 at.%), NiAk (line compound, 
A.Y =! 0 at.% ) and NisAl (small phase field, AX = 4 at.%). Room temperature 
ion-beam mixing in multilayer samples with composition 3:1 and 1:1 form a crystalline 
disordered cubic phase NizAl and crystalline NiAl phase respectively. But for 1:3 
composition an amorphous NiAlz phase is formed[25] as this is a line compound with 
complex (orthorombic) crystal structure. However, for Sb/Al systemsfST] the SbAl 
rhombohedral phase is formed in case of bilayer sample at 77 K whereas in multilayer 
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case an amorphous phase is formed. It may be noted that Sb/Al system has a single 
line compound SbAl in the phase diagram. Brimhall et al.[59] have pointed out that 
the phase field wider than 10 at.% are stable whereas they become unstable below a 
few atom percent. In our case Cu^Ge, Ni 2 Ge and Co^Ge phases have phase widths 
AA = 4 at.%, 0 at.% and 9 at.% respectively. We suggest that the formation of stable 
compound phases due to irradiation at room temperature and above are probably due 
to the high atomic mobilities of the atoms. There are many other reports of stable 
compound phase formation {Co^Si^ Pd 2 Ge, Pd 2 Si etc. ) for these kind of narrow 
phase fields. The lowest is the TiRu which has a phase, field width 0.8 at%[60]. 

In summary, we have shown that the formation of crystalline phases is possible 
in these systems depending on certain conditions viz., temperature, structure, atomic 
size, phase width etc. In the following section we shall discuss the formation of “first 
phase" by ion-beam mixing and thermal annealing with the help of existing rules. 


4.4.3 Formation of “First phase” 

In many metal/semiconductor systems[42, 57, 88-90], it has been found that when 
one studies the sequences of phase formation by thermal annealing and ion-beam mixing 
the phase which synthesizes first in both processes is same in most of the cases. We 
have sho^vn in Table 4.17 the comparison of the phases produced by ion- beam mixing 
and thermal annealing in metal/semiconductor systems. Our results for germanide 
systems (also presented in Table 4.17) clearly show that the “first phase” produced by 
thermal annealing and by ion-beam mixing are same for all the three systems. 

Similarity of ion-beam mixing and thermal annealing results in metal/semiconductor 
systems can be viewed through the phase diagram. The correlation between “first 
phase produced by thermal annealing in a bilayer sample and the phase diagrami was 
proposed by Waiser and Bene[9S]. Their rule states that the phase to grow first upon 
thermal annealing of a bilayer sample is the highest congruently melting compound 
next to the lowest melting eutectic one in the phase diagram. On a similar line, for 



127 


ion-beam mixing in metal/Si systems, Lau et al.[57] proposed that the “first phase” 
which would produce, is the highest congruently melting compound next to the lowest 
melting eutectic in the phase diagram provided the atoms have sufficient mobility to 
form compound. 

From Cu-Ge phase diagram (Fig 4.26) we find that the lowest temperature eutectic 
is at the position of 63 at.% of Cu shown by dashed vertical line. Therefore, it is 
expected that Cu^Ge compound should form first. This is consistent with our thermal 
annealing and ion mixing results. In case of Ni/Ge system, the eutectic point is at 67 
at.% of Ge so the first phase is expected to be NizGe^- But our experimental results of 
ion-beam mixing and thermal annealing show that it is the Ni^Ge phase which forms 
first. Similarly for Co/Ge system this rule predicts that the CoGe 2 phase should form 
first whereas we obtained Co 2 Ge phase in our experimental studies. Such violations 
of the prediction of Lau et al or Walser and Bene rule have been reported for other 
silicide systems viz. Co/Si[27] etc. 

Let us consider the sequence of phase formation on the basis of thermodynamic 
point of view according to which “the phase corresponding to highest negative heat of 
formation has the largest chemical driving force and hence should form first” [27]. From 
Table 4.19 the heat of formation of Cu^Ge phase is - 16 kj/mole of atom[99] while the 
value for 4 " phase is not available. From the phase diagram (Fig 4.26) it is seen that the 
melting temperature of ^ phase is slightly higher than that of CuzGe phase signifying 
that the heat of formation of ( phase should be higher than that of CuzGe phase. Thus, 
from the thermodynamic point of view the ( phase should form first. In case of Ni/Ge 
system the values of heat of formation of iViaGe, NiGe and NizGe phases axe -37, -32 
and -25 kJ/mole of atom respectively (Table 4.19). Therefore, in this case NizGe is 
expected to form first which is in accordance with the ion-beam mixing and thermal 
annealing experimental results. For Co/Ge system, the values of heat of formation for 
CoGe, CozGe and CoGez being -16, -15 and -12 kJ/mole of atom respectively, CoGe 
phase should form first. But experimental results show that CozGe has formed first. 
Similar deviations (like Cu/Ge and Co/Ge) from the heat of formation rule have been 
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reported in many metal/Si systems[27]. However, it is to be noted that the difference 
in energy values of the phases in all such cases being very close to each other, any 
phase can appear as a first phase. Thus the thermodynamics alone cannot predict the 
formation of phases in these cases. In short, none of the criteria mentioned above could 
simultaneously explain the formation of the “first phase” in all the three cases. 

The similarity of ion-beam mixing and thermal annealing in the formation of equi- 
librium phases suggests that both thermodynamic and kinetic considerations[100] have 
to be taken into account* since in case of thermal annealing, it is assumed that in 
determining the sequence of phase formation, kinetic factors dominate over the ther- 
modynamic driving forces[27] where as. in ion-beam mixing the kinetic restriction is 
normally relaxed. In the following, a rule proposed by Pratorius based on the cal- 
culation of effective heat of formation [56] which considers both thermodynamics and 
kinetics aspects will be discussed. 

The kinetic barrier and effective heat of formation are both related to the number of 
atoms taking part in the reaction which in turn are related to the atomic mobility and 
supply of atoms during phase formation. The effective concentration at the reacting 
interface could be influenced by many other factors such as diffusivity, temperature etc. 
Pretorius et aJ. assumed that the effective concentration at the reacting interface is 
the concentration of the lowest temperature eutectic. Based on this concept the rules 
for phase formation proposed by Pretorius et al. are 

1. The first compound phase to form during interaction of two elements is the phase 
having the largest negative effective heat of formation (AH'). 

2. The second phase to form at the interface between the phase and the remaining 
element is the next phase richer in the unreacted element which has the most 
negative heat of formation. 

The effective heat of formation AH' i.e., the total heat released in the reaction is 
equal to the product of the heat of formation and the total number of atoms used in 
the reaction. One can write AH' as [27] 
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T he compound concentration of the limiting element 
It can be seen from the above expression that for each phase, the most negative 
AH' and therefore the release of the highest energy from the system occurs, when 
the atomic concentration matches with a particular compound i.e. AH' = AH. The 
eifective heat of formation diagrams can be easily constructed by plotting the heat of 
formation AH and its compositional concentration and completing the triangulation 
by connecting those points to the end points of the concentration axis. 

Figure 4.30a shows the AH' diagram calculated for Cu/Ge system. The AH 
value for CuzGe is -16 kJ/mole. The maximum value of AH’ is thus -16 kJ/mole 
which is shown by point ‘a’ in the figure. The triangle is constructed by connecting 
this point to 0 and 100 on the concentration axis. In the absence of the value of 
AH for CuzGe phase we have assumed that the the difference of AH between Cu^Ge 
and CuzGt phases is small because of the relatively small change in composition and 
melting temperatures[47]. For the purpose of estimation, we have considered a slightly 
higher AH value for CusGe phase than that of CuzGe phase as shown by the point 
‘b’. The actual value might be different from this but certainly it is not below the 
value of CuzGe phase. It may be noted that the value of AH' for CuzGe phase is 
larger for Cu concentration below 76%. Therefore, according to the Pretonus’ rule, 
in this case CuzGe phase should form first as the effective heat of formation AH is 
larger for CuzGe at the eutectic concentration of 60% Cu. However, since the actual 
AH value of Cu^Ge phase is not correctly known and for Cu concentrations below 
76% the difference in the AH' values is small, it possible that either CuzGe or CuzGe 
could form at the interface. In reality, the actual concentration of the liimting etement 
available at the interface and the formation of first phase depends on several factors 

viz. diffusing species, their mobilities, presence of other impurities etc. It appears that, 

in the present case of Cu/Ge system these factors contribute in such a manner that 
the “available concentration of the limiting element” (i.e. Cu) favours the formation 
of CuzGt' phase at the interface. 
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The A//' diagram fi>r Ni-(5e system is shown in Fig 4.30b which shows that at 
the eutectic coiuxmtrat ioti ol titi'X Ge. NiGc phase has the largest negative heat of 
formation and hence* \vt>ul<i ft^rin first. Bui we have observed that Ni 2 Ge phase forms 
first in both ion-ixvini mixing or thermal annealing. The NiGe phase is produced 
subsequent to formation, litis can be explained from the AH diagram by 

considering the supply of atoms which actually determine the available concentration 
across the interface. It is s«*en form fig 1.30b that Ni^Ge has largest negative heat of 
formation for Ge conct*ntrati<m of less than 47%. Above 47% Ge concentration NiGe 
has largest negative h<*at of formation. .Vt 33% of Ge, Ni^Ge phase has highest effective 
heat of formation. In t his case. Ni and Ge are almost equally mobile (see Chapter 5) 
The occnrn'iice of Si /it as a first phase cati be seen from the following calculations. 
For illustration. w<‘ calculate the effective lurat of formation AH’ by considering that 
25% of total ( «e at«ims are available for reaction at the reacting interface. The effective 
heat of formation fc>r .N’i( le pluuie at the same interface is 

AH' =42 (25 %fWA>) = -i(» kJ/nio!eof atom. 

For Sidt'i pliiwe the effeetivt* heat of formation at the same interface is 

AH' =-25 (25 %/(>(>%)) = -9.5 k.l/mo!e of atom. 

Similarly, for Si/sf phase the effective heat of formation at the same interface is 

AH' =-37 (25 %/33%) = -2S kJ/mole of atom. 

The formation of Ni/it phase is followed by the formation NiGe phase (which has 
second highest effective heal of formation for Ge concentration less than 47%) at the 
NiiGe. and Ge interface because at this interface Ge supply will be more. 

For Co/Ge system, the AIV diagram is shown in Fig 4.30c. In this case CojGe 
phase has the largest AH' for Ge concentration less than 41%, CoGe 2 has most negative 
heat of formation for the Ge concentration greater than 65% and CoGe phase has most 
negative heat of formation between 42% and 64% of Ge. The heat of formation for 
these three phases are almost equal. In this case Ge is the only moving specie. On the 
basis of calculations similar to the one performed for Ni/Ge system it can be shown 
that in case of Co/Ge syHt«*in, CoaGc phase will appear first followed by the growth of 
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CoGe phase at the Ge and Co 2 Gt interface. 

To snmmarize, in this section it has been shown that the formation of first and 
subsequent phases in metal/Ge systems can be understood on the basis of effective 
heat of formation idea and the mobility of the species. 


4.5 Conclusions 

The conclusions of this chapter are given below. 

1. It has been shown for the first time that MeV ion-beam mixing is an efficient low 
temperature process for producing stable equilibrium phases in Cu/Ge, Ni/Ge 
and Co/Ge systems. 

2. In the case of multilayer samples, ion-beam mixing as well as thermal annealing 
lead to the formation of Cu^Ge Ni 2 G& and C 02 G& as ffirst phase”. In case of 
Cu/ Ge system, another Cu rich phase CusGe, has also been observed at higher 
doses. 

3. In case of bilayer samples, thermal annealing produces C uzGe^ NizGe and Co^Ge. 
phases as "first phase”. At higher annealing temperatures above 475 K NiGe and 
CoGe phases have been observed as “second phase” in Ni/Ge and Co/Ge systems. 

4. The first phase formed by thermal annealing and ion-beam mixing is same in 
all the cases. 

5 . The formation of crystalline phases by ion-beam mixing in these systems is in 
accordance with the prediction of the existing empirical rules. 

6. The first and subsequent phase formation can be explained on the basis of the 

effective heat of formation rule and by considering the mobility of the atomic 
species. 
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7. The resistivity values of CuzGe. NizGe and Co^Ge layers produced by ion-beam 
mixing are 20, 50 and 74 /xfl! — cm respectively. The resistivities of the phases 
produced by thermal annealing are 8.6, 23 and 50 /xfl — cm respectively. It is 
important to note that these resistivity values are lower as compared to metal sili- 
cides. Hence it is expected that these germanide phases would become important 
for semiconductor industry. 



Chapter 5 


Mobility of species during 
ion-beam mixing 


5.1 Introduction 

In the previous chapter we have shown the formation of compound phase by ion- 
beam mixing. It has also been argued that it is the dominant moving species that 
control the first and subsequent phase formation. In this chapter^ we will describe the 
experimental results of the atomic movements across the M/Ge interface during the 
sjTithesis of germanide phases by ion-beam mixing. 


^Part of this work is published in Semiconductor Devices Ed. Krishan Lai, Interna- 
tional conference on semiconductor devices, India, 1996 (Narosa Publishing House, New 
Delhi) p 502. 
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5.2 Experimental 

Bilayer samples {SiO^I M f Au{marker) f Ge) having a Au film of thickness less than 1 
nm acting as a marker for monitoring the atomic motions across the interface, were 
prepared by e-beam evaporation technique. In all these samples Ge was deposited as 
the top layer. 

The ion-beam mixing at different temperatures were performed by 1 MeV Kr"^ 
and Ar'^ ions for different doses in the range of (2 — 20) x 10^’ ions! cm?. At this 
energy most of the irradiating ions are transmitted through the bilayer structure. The 
deposited energy distribution was calculated using TRIM code[78]. The Rutherford 
back scattering (RBS) measurements with 1. 3-1.5 MeV He^ beam were carried out 
to study the composition of the mixed layers and movement of metal and Ge atoms 
across the interface by monitoring the position of the marker. The RUMP simulation 
code[81] was used to simulate all such RBS spectra. 


5.3 Experimental reisults 

Typical RBS spectra of a sample [Si02f C uf Au(^marker) f Ge] in the as prepared con- 
dition and after performing mixing at room temperature using Kr ions for a dose of 
1 X 10^® Kr^ foT? are shown in Fig 5.1. The simulation analysis shows the composition 
of the mixed region as C tiT 5 Ge 25 (or CuzGe). This type of growth of germanide phase 
requires a mass flow of one or both species across the marker leaving a different amount 
of material in front of the marker. This is equivalent to have a cover layer which results 
into an energy loss producing a shift in the position of the marker signal in the RBS 
spectrum. In this case Au marker has moved towards lower energy side after mixing. 
Also, there is a spreading in the Au peak which is due to the movement of Au atoms 
in Ge and Cu layers. However, the peak position is more important for evaluating the 
relative migration of the atoms of the two adjacent layersjlOlj. 

The movement of the marker signal is determined by measuring the energy dis- 
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Figure 5.1; The RBS sepctra of Cu/Ge sample embedded with Au marker layer at the 
interface. The continuous line shows RBS spectrum of as prepared sample while the 
dashed line spectrum is taken after ion-beam mixing. The Au marker shifts towards 
low energy side suggesting dominating motion of Cu atoms. 
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placement {AE) of marker peak position. For measuring AE we have followed a sign 
convention since it can shift in either lower or higher energy side after mixing. The 
negative sign will be used for the marker movement towards the lower energy side and 
positive sign for movement towards the higher energy side. Analysis of the spectra pre- 
sented in Fig 5.1 gives the energy shift of the marker position as -16 keV. The marker 
is located in the mixed region ( CuzGe ) of total thickness 25 nm of which 17 nm is on 
the higher energy side of the marker as shown in Fig 5.1. This gives the ratio of the 
flux of Cu and Ge across the interface as 7:1. Other samples irradiated for different 
doses using Kr and Ar ions, whose RBS spectra are not shown here*, also show marker 
movement towards the lower energy side. This shift increases with increasing ion dose. 
The shift in the marker position is plotted in Fig 5.2 as a function of the thickness of 
the germanide layer. The top(bottom) solid line is drawn from the calculated value for 
the extreme case where only Ge(Cu) atom is mobile. These results show that the Cu 
and Ge both are mobile but Cu is the dominant mobile species. d’Heurel et al.[47] have 
reported similar movements of Cu and Ge in thermal annealing experiments performed 
with Xe marker. 

In case of Ni/Ge system, the marker position remains stationary which can be 
seen from RBS spectra of Fig 5.3 for a sample before and after irradiation at room 
temperature for a dose of 6 x lO^^ Xrlcm^. This feature has been found for both Kr 
and Ar ion irradiations for all the doses used in this work. This implies that in case of 

bii/Ge system both Ni and Ge are almost equally mobile during the growth of NizGe 
phase. 

For Co/Ge system, the marker movement is just opposite to that of Cu/Ge system. 
The RBS spectrum of a Co/ Ge sample after irradiation for a dose of 1 x 10^® Kr /cm? is 
shown in Fig 5.4 along with the as-deposited spectrum. Here, the Au peak position has 
shifted towards higher energy side. The composition of the mixed region is CojGe.' The 
plot of vursus thickness of the mixed (cobalt germanide) layers is shown in Fig 5.5. 

In this case, the experimental points coincide with the theoretical line calculated by 
assuming the movement of Ge only. This suggests that for Co/Ge system, Ge atom is 
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Figure 5.2: The plot of energy displacement AE of the marker vursus thickness of 
the CuzGe layers. The upper (lower) continuous line is obtained from the theoretical 
calculations by considering that only Ge (Cu) atoms aire mobile. 
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Figure 5.3: The RBS sepctra taken before and after ion-beam mixing at RT for a dose 
of 6 X Kr^fcm^ of Ni/Ge sample with a Au marker. 
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the only moving species. Figure 5.6 shows the comparison of RBS spectra of thermally 
annealed marker samples at different temperatures. Annealing at 480 K for 60 min. 
shows very little marker movement towards the surface (higher energy side) but at a 
higher temperature of 580 K when both the layers are completely consumed (forming a 
CoGe phase) the marker is found on the surface of the bilayer suggesting the movement 
of Ge atoms only. 


5.4 Discussion 

We have mentioned earlier that the ion-beam mixing process can be divided into two 
regimes (1) prompt regime (2) delayed regime. In the prompt regime the incident ions 
transfer energy to the target atoms producing number of primary recoils in the forward 
direction and is considered to have a negligible effect to ion-beam mixing[14]. Mixing 
in the prompt region is mainly due to the secondary recoils which are isotropic. In a 
collision cascade the lighter target atoms always have a larger displacement as compared 
to the heavier ones. The experimental results of Tao et al.[102] on twelve binary 
systems show that the lighter species tend to migrate towards the heavier elements 
and is correlated with the damage energy via inverse Kirkendell effect. This motion 
can also be seen from the spike point of view[103]. In a spike the temperature is 
instantaneously very high for a short duration within a small volume. The reaction 
in the spike may be treated as a normal thermal process for the duration of that high 
temperature. This temperature is much higher than the substrate temperature so the 
variation of ambient temperature should not affect the atomic motion in the spike. 
This has been shown for Ni/Si and Pt/Si systems[28]. 

Apart from these collisional effects the movement of atoms is also influenced by 
the thermodynamic forces[28, 102] which are similar to those operating during ther- 
mal thermal annealing. However, the collisional effects are totally absent in the case 
of thermal annealhag. In ion-beam mixing if the two components namely collisional 
and thermodynamic operate in the same direction then only one element is observed 
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Figure 5.4: The RBS sepctra of Co/Ge sample with a Au marker taken before and 
after ion-beam mixing at RT for a dose of 1 x 10^® A'r+/cm^. 
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Cobalt germanide thickness (nm) 


Figure 5.5: The plot of energy displacement AE vs. thickness of the cobalt -^rmanide 
layers. The upper (lower) continuous line is obtained from the theoretical calcxilations 
by considering that only Ge(Co) atoms are mobile. 
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Figure 5.6: The compaxison of RBS sepctra of Co/Ge marker sample taken before and 
after thermal annealing at two different temperatures of 480 K and 580 K. 
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as the moving species and if they operate in opposite direction movement of both 
elements is expected [102]. This facts can be seen from the experimental data of 
metal/semi conductors systems given in Table 5.1. 

In order to see the effect the collisional process, the distribution of damage energie 
obtained from the TRIM calculation[79] is shown in Fig o.7a for Cu/Ge case. According 
to the inverse Kirkendall effect[102] it is expected that Ge atoms will move in the Cu 
layer. On the other hand, the results obtained from the thermal reactions[47] show 
that Cu is the only moving species for the CuzGe phase formation. Thus, for Cu/Ge 
the two factors namely collisional and thermodynamic are acting in opposite direction 
resulting into the motion of both Cu and Ge. Our calculation of the ratio of the flux of 
Cu to Ge which is 7:1 which shows that the thermodynamic processes dominate during 
mixing as par as the mobility of the atoms is concerned. The equal mobility of Ni and 
Ge for the case of Ni/Ge systems can be explained in analogous manner. 

In case of Co/Ge system Ge should be mobile under collisional effects as predicted by 
the damage distribution shown in Fig 5.7c. Motion of only Ge atoms has been revealed 
by thermally annealing experiments. These results suggest that both collisional and 
thermodynamic effects are acting in the same direction which result in the mobility of 
Ge atoms only which is in accordance to our observations. 


5.5 Conclusion 

In summary, we have shown that during the formation of CuzGe- phase by ion-beam 
mixing at room temperature, both Cu and Ge are mobile but Cu is the dominant 
moving species across the interface while for Ni/Ge system both Ni and Ge axe equally 
mobile for the formation Ni^Ge phase. However, for Co/Ge system only Ge atoms 
move to produce Co^Ge phase. 
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Table 5.1: Reported results of marker movements in metal/Si systems during thermal 
annealing and ion-beam mixing 


System 

Marker 

Moving species 

Ion-beam mixing 




and flux ratio 




formed 



under thermal 

Flux species 

Atoms whose 

Layer of 




annealing 

and moving 

motion is rela- 

lower damage 





ratio 

tive enhanced 

density 


Si/Ni 

W. Ag 

Ni only 

M ~ 1 - 2 

Si ~ 2 1 

Si 

Si 

NiiSi 

Si/Pt 

Mo 

^ 13 

Si — 1 

EL ~ 1-1 

St ~ 1 2 

Si 

Si 

Pt2Si 

Si/Pd 

W 

E£ > 1 

f <1 

Si 

Si 

Pd2Si 

Si/Mo 


Si only 

Si only 

Si only 

Si 

MoSi2 

Si/W 


Si only 

Si only 

Si only 

Si 

WSi2 

Si/Cr 

W 

Si only 

Si onlv 
* 

Si only 

Si 

CrSi2 

Si/Mg 

W 

Mg only 

Mg only 

Mg only 

Mg 

M g2Si 

Ge/Ti 

W 

Ge ~ 15 

Ti — 1 

Ss. Z 

Ti — X 

Ti 

Ti 

Ti ^ 6 

Ge - 5 

Ge/Mg 

W 

Mg only 

Mg only 

Mg only 

Mg 

Mg2Gt 

Cu/Ge {*) 

Au 

Cu only 

I 

Ge 1 

Ge 

Ge 

Cu^Gc 

Ni/Ge (*) 

Au 


Ni^i 

Ge 

Ge 

NiiGe 

Co/Ge (*) 

Au 

Ge only 

Ge only 

Ge 

Ge , 

Co2Gt 


* represents present work and other data from ref 101 
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Figure 5.T; The distribution of damage energr- due to 1 MeV Kr ion bomb^dment as 
a function of depth in (a) Cu/Ge and (b) Ni/Ge and (c) Co/Ge systems. Each point 


(x) corresponds to energy deposition in 10 nm 



Chapter 6 


Mechanism of ion-beam mixing 


6.1 Introduction 

The formation of crystalline compounds in multilayer configuration of M/Ge (M = 
Cu, Ni, Co) systems by ion mixing has been dealt in Chapter 4. Some basic issues 
concerning the migration of the different species during mixing have been discussed in 
the last chapter in the light of the experiments performed using marker configurations. 
In this chapter^ we discuss the mechanism of ion mixing in M/Ge systems on the basis 
of the experimentally determined parameters such as mixing rate, mixing efficiency etc. 
For this purpose bilayer configuration is the most suitable, since accurate quantification 
of these parameters is possible in this configuration only. To understand the basic 
processes of mixing mechanism, the prerequisite is to know the dependence of mixing 
on ion mass, ion energy, fluence, heat of mixing, cohesive energies, mobilities of the 
elements, composition of the mixed region etc. 

In this chapter the experimental results of ion mixing of M/Ge bilayer systems by 
Kr and Ar ions at different temperatures are first described. The critical temperature 


‘Part of this work is published in Phys. Eev. B 54 (1990) 5769. 
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Tc above which thermally activated processes set ia for these systems, have been de- 
termined. Empericcd models[20, 23, 51] which estimate the critical temperatures have 
been discussed. 

The experimentally obtained room temperature mixing rates are then compared 
with the predictions of the existing phenomenological models[19, 22, 24, 31, 50] which 
show that the existing models cannot explain our results. A new model for compound 
formation has been proposed which considers mixing to occur under non-overlapping 
subcascades. The predictions of this model are shown to be in good agreement for 
highly reactive medium Z bilayer compound forming systems (such as Cu/Ge, Ni/Ge, 
Co/ Ge etc. ) where mixing rate is proportional to the deposition energy Fd • In addition, 
it has been shown that this model also able to explain the observed mixing rate of Ni/Si 
system reported in the literature[31]. 


6.2 Experimental 

Metal/Ge samples having bilayer configuration (Fig Lid) were prepared by sequen- 
tially depositing high purity metal M (Cu, Ge, Co) and Ge films on clean fused quartz 
substrates, without breaking vacuum, using e-gun evaporation system. The total thick- 
ness of these bilayer samples was in the range of 150 nm to 270 nm with an overall 
composition in the vicinity of MsoGcsq. 

Ion-beam mixing was performed in the range of 100-450 K. using 1 MeV Kr+ and 
in the dose range of (2 — 20) x 10^® ions/cm-. The incident ion flux was kept 
low (< 5 X 10^2 ions/cm^s-'^) to minimize the rise in the sample temperature. It 
has been reported that MeV ion irradiations cause a typical rise of about 10 K of the 
film temperature and the temperature rise in Si 02 substrate was believed to be higher 
by another ~ 10 K, suggesting actual rise of sample temperature by 20 K[74]. The 
values of range and deposited energy densities by the irradiating ions in the target 
were calculated by Monte Carlo computer program TRIM[78]. In the present case, 
most of the irradiating Kr and Ar ions (1 MeV) are transmitted through the bilayer 
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structure. H.utherforcl backscattering (RBS) measurements with 1.2-L5 MeV 
beam were carried out to determine the composition and growth of the mixed region 
of irradiated samples. RLMP code[81] has been used to simulate all experimental 
RBS spectra. RBS spectra of all samples did not show any noticeable oxygen or 
carbon concentrations. X-ray diffraction measurements were performed using Cu — 
source. However, the intensity of the peaks were very weak due to small volume of 
the mixed region. Therefore, the composition obtained from RBS analysis of these 
bilayer samples after mixing arc used to infer about the nature of the phases formed by 
correlating with corresponding compounds (determined by X-ray analysis) obtained 
in these systems after ion mixing in multilayer configuration (presented in Chapter 
4). The surface morphology for all the samples has been studied by scanning electron 
microscopy whicli did not show any noticeable feature before and after ion mixing. 


6.3 Results 

This section describes temperature dependence of ion-beam mixing in Cu/Ge, Ni/Ge 
and Co/Ge systems in the range 100 K to 450 K. 


6.3.1 Cu/Ge system 

The RBS spectrum of an as-deposited Cu/Ge bilayer sample is shown in Fig 6.1 along 
with tile simulated spectrum. This sample consists of a top Ge layer of thickness 100 
nm and bottom layer of Cu having a thickness of 50 nm. The simulated structure 
used for this is shown in Table 6.1. Figure 6.1 also shows the RBS spectrum of the 
sample taken after ion-beam mixing at room temperature (RT) for a dose of 6 x 10 
Kr+fcni^. The simulated spectrum for the irradiated case is obtained by using a 
layer structure shown in Table 6.1. The simulation analysis shows a stoichiometnc 
composition of across the interface which is indicated by the step in the 

copper signal at channel number 365. The formation of the stoichiometric composition 
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is clearly seen in Fig 6.2 which shows the RBS spectrum of the sample irradiated at 
a higher dose of 8 x 10^® Kr'^fcm?. The simulated spectrum is obtained by using 
two adjacent layers of the mixed region one with an overall composition of Cur 5 Ge 25 
and the other of composition Cus^Geir. The layer structure utilized for simulation is 
given in Table 6.1. Similar results are observed for Ar ion irradiations performed at 
room temperature. However, the thickness of the mixed region in this case is lower as 
compared to the one obtained for Kr irradiated samples for the same dose. 

The squared thickness (X^) of the mixed layers of Gu^Ge phase produced at room 
temperature due to both Kr and Ar ion bombardment is plotted as a function of dose 
in Fig 6.3. The slope of the linear dependence gives the value of mixing rates (X^/$) 
as 20.6 nm'* and 5.4 for Kr and Ar ions respectively. It is to be noted that the 
thermal annealing of this bilayer samples also show (in Fig 6.4) a linear dependence 
of the square of the reacted region with annealing time[48] which suggests that the 
mi.xing process is like a dilfusion controlled process. The variation of mixing rate with 
the deposited energy density (Fd) at the interface is shown in Fig 6.5 and it is observed 
that the mixing rate is linearly proportional to Fd. Similar behavior has been reported 
for many medium Z metal/metal[25, 66, 76] and metal/Si systems[31, 75]. The slope 
of this line gives mixing efficiency (A^/ ^Fd) of 13.1 nrrF JkeV . This mixing parameter 
is a useful measure of the ion mixing since it normalizes the mixing rate with respect 
to deposited energy, thereby providing a convenient way to compare the results of 
different systems performed under different irradiation conditions. 

The comparison of ion-beam induced effects in Cu/Ge bilayer samples irradiated 
for a dose of 4 x 10^® Kr'^ jam? at different irradiation temperatures are shown by the 
RBS spectra of Fig 6.6. At sample temperatures above RT (300 K), mixing occurs at 
much higher rate. However, the mixing rate at 100 K occurs at almost similar rate as 
in the case of RT mixing. The mixing rates obtained at RT and 100 K do not differ by 
more than 20%. The temperature dependence of the mixing rates (normally referred 
to as Q curve) for Ar"*" ion is plotted in Fig 6.7a. This figure clearly shows that the 
mixing is nearly temperature independent up to 328 K with a steep increase at higher 
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Figure 6.2: Rutherford baokscattering spectra of Cu/ Ge bilayer thin film on quartz 
substrate for as-deposited condition (Q) room-temperature mixing (A) with 

1 MeV Kr ions at a dose of 8 x 10^® Krjcm^. The corresponding simulated spectra 
are shown by solid lines. 
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Table 6.1: Description of the layer structure used for the simulation of RBS spectra of 
as-deposited and ion mixed Cu/Ge bilayer samples shown in Figs 6.1 and 6.2. 


Sample 

description 

Layer structure used for simulation 

Layer Thickness Composition 

No. (nm) 

Cu/Ge bilayer 
film deposited 

on SiOz 

1 165 Ge 1.00 

2 105 Cu 1.00 

3 1000 Si 1.00 0 2.00 

After ion mixing 

at RT for a dose of 

6 X 10^5 Kr/cm^ 

1 140 Ge 1.00 

2 34 Cu 0.75 Ge 0.25 

3 SO Cu 1.00 

4 1000 Si 1.00 0 2.00 

After ion mixing 

at RT for a dose of 

8 X 10^® Krfcm? 

1 135 Ge 1.00 

2 30 Cu 0.75 Ge 0.25 

3 30 Cu 0.83 Ge 0.17 

1 54 Cu 1.00 

4 1000 Si 1.00 0 2.00 
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Time (min ) 


Figure 6.4: The squared thickness of the reacted layer of composition CusGe formed 
across the interface of Cu/Ge bilayer system after thermal annealing is shown as a 
funtion of annealing time for different temperatures (after ref 48). 
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temperatures. The critical temperature which marks the beginning oftb^ 

• 4 e temperature 

dependent mixing is around 320 K. It is to be noted that for this system r 

^ ’ , tne formation 


ove 400 K[47]. 


of the equilibrium phase after thermal treatment has been observed abcm^ 

However for Ar ions, mixing rate does not exhibits sharp temperature de^ pendence as in 
the case of Kr. The values of the mixing rates in case of Ar are lower foj^-^j. similar doses 
but the critical temperature does not change appreciably. The ^ activated 

contribution can be obtained by subtracting the smallest mixing rate ^ lowest 

irradiation temperature of the Q-curve from the remainder of the curve. ^ resulting 

plot for the temperature dependent contribution is shown in Fig 6.7b which the 

activation energy £4 of ion mixing of approximately 1 eV is determineci^^ rpj^-g 
close to the one reported for thermal diffusion[48]. 


6.3.2 Ni/Ge system 

Figure 6.8 shows the RBS spectrum of as-deposited Ni/Ge bilayer samp a^jQjjg 3 ^ 
spectrum taken after performing ion-beam mixing at RT at a dose of 8 x 7 ; ^ 0^5 /an^. 

The step in the Ni signal at channel number 374 in case of irradiated specr^^t. mm indicates 
the formation of a stoichiometric composition at the Ni/Ge interface. sjmulated 

spectra are shown by solid lines and the simulated structures are give:^^ Table 6 . 2 . 
The simulation analysis of the irradiated sample shows the formation of> ^ composition 
of NieeGcM similar to the multilayer Ni/Ge sample as described in 4 

analysis of ion mixed samples at other irradiation doses (at room gy. 

and Ar ions show formation of mixed layer of similar composition. This composition is 
belong to one of the stable compound phases, Ni^Ge in the equilibrium phase diagram 
(as shown in Fig. 4 . 27 of Chapter 4). Note that we have identified thi^ compound in 
case of ion-beam mixing of multilayer samples as well. 

The dose dependence of the squared thickness {X^) of the mixed laye^rjs produced by 
Kr and Ar ion bombardments at RT are shown in Fig 6.9 where for both t::;he cases linear 
variation is observed which is similar to the case of Cu/Ge system. F^^r comparison, 



Energy (MeV ) 


Si02/Cu /Ge (bilayer) 
o As deposited 

Ion mixed {4x lO'^KrVcm 2) 
X 100 K 


1-3 MeV He 
G =150° 


350 

Channel 


Figure 6.6: Rutherford backscattering spectra of Cu/Ge bilayer film on quartz substrate 
for as-deposited condition and after mixing for a dose of 4 x 10^® Kr/cm^ at 100 K, 
RT and 350 K. 
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Fl6u« 6,T: Wn« of At.h«3« pW 


Fig.6.6a. 
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Energy ( MeV) 



Figure 6.8: Rutherford backscattering spectra of Ni/Ge bilayer thin film on quartz 
substrate for cis-deposited condition and after room-temperature mixing with 1 MeV 
Kr ions at a dose of 8 x 10^® Krfcm^. The corresponding simulated spectra are shown 
by solid lines. 
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Table 6.2: Description of the layer structure used for the simulation of RBS spectra of 
before and after ion mixing in Ni/Ge bilayer sample shown in Fig 6.8. 


Sample 

description 

Layer structure used for simulation 

Layer Thickness Composition 

No. (nni) 

Ni/Ge bilayer 
film deposited 

on Si 02 

1 111 Ge 1.00 

2 60 Ni 1.00 

3 1000 Si 1.00 0 2.00 

After RT ion mixing 

for a dose of 

of 8 X 10^® Krfcm? 

1 82 Ge LOO 

2 46 Ni 0.66 Ge 0.34 

3 37 Ni 1.00 

4 1000 Si 1.00 0 2.00 
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Fig 6.10 shows the variation of square of the reacted layer thickness as a function 
of different annealing durations at 480 K. The values of mixing rate (X^/^) obtained 
from the studies of dose dependence are found to be 18.5 and 5.6 for Kr and 
Ar ions respectively. The dependence of mixing rate on the deposited energy density, 
Fd, at the interface is shown in Fig 6.11. The mixing rate varies linearly with Fd and 
mixing efficiency obtained from the slope of this curve is 11.3 nm^/keV. 

The RBS spectra of Fig 6.12 compare the effect of Kr irradiation at 400 K for a 
dose of 6 X 10^® Kr'^ (cnrF with the as-deposited one. In this case the amount of mixing 
is much higher than that obtained at RT or 100 K for the same dose. The Q curves for 
Kr ion irradiations are plotted in Fig 6.13a. The critical temperature for this system 
is 326 K and is shown by an arrow. Only the temperature dependent part of mixing is 
shown in Fig 6.13b from which the value of activation energy Ea of 0.3 eV is obtained. 


6,3.3 Co/Ge system 

The RBS spectra of Fig 6.14 show the as-deposited Co/Ge bilayer sample along with the 
spectrum taken after ion-beam mixing at a dose of 1 x 10^® AV^/cm^. The simulated 
spectrum (solid line) for the irradiated case is obtained using a composition of Coo^GtzA 
which corresponds to the Co 2 Ge phase (see Fig 4.28 of Chapter 4). The detail of the 
layer structure used for simulation of the as-deposited and irradiated samples are given 
in Table 6.3. The variations of X^ as a function of incident Kr and Ar dose are shown 
in Fig 6.15 whereas Fig 6.16 shows the variation of (Af^) as a function df annealing 
duration. This suggest that the mixing process is like a diffusion controlled process. 
The dependence of mixing rate on F^ is shown in Fig 6.17. The value of mixing 
rates at room temperature are found to be 7.7 nm^ and 1.6 ntn^ for Kr and Ar ions 
bombardment respectively. The slope of the line in Fig 6.17 gives mixing efficiency 
of 4 nm^/keV. The Q curve for Kr ion irradiations is plotted in Fig.l7. The critical 
temperature for this system is 347 K and is shown by arrow in Fig 6.1Sa. The activation 

energ}' Ea of 0.15 eV is obtained from the temperature dependent part as shown in 
Fig 6.18b. 



X^(10^^cm2) 
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Figure 6.9: Dose dependence of squared thickness of mixed layers of composition NkGe 
formed across the interface of Ni/Ge bilayer system irradiated at room temperature 
with 1 MeV Ar ions (Q) and Kr ions (#) respectively. 
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Figure 6.10: The squared thickness of the reacted layer of composition Ni^Ge formed 
across the interface of Ni/Ge bilayer system after thermal annealing is shown as a 
funtion of annealing time. 
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Figure 6.12: Rutherford backscattering spectra of Ni/Ge bilayer thin film on quartz 
substrate for as-deposited condition and after mixing for a dose of 6 x 10^® Kvjcn^ at 
410 K. 
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Figure 6.13: The vaxiation of mixing rate for Ni/Ge bilayer system as a function of 
reciprocal of substrate temperature for a dose of 6 x 10^’ A'’r/cm^ (b) Arrhenius plot 
of the mixing rate as a function of reciprocal of substrate temperature extracted from 
Fig.6.12a. 
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Table 6.3: Description of the layer structure used for the simulation, of RBS spectra of 
as-deposited and ion mixed Co/Ge bilayer samples shown in Fig 6.13. 


Sample 

description 

Layer structure used for simulation 

Layer Thickness Composition 

No. (nni) 

Co/Ge bilayer at 
film deposited 

i 

on Si 02 

1 107 Ge 1.00 

2 65 Co 1.00 

3 1000 Si 1.00 0 2.00 

After RT ion mixing 

for a dose of 

1 X 10^® Kr/av? 

1 94 Ge 1.00 

2 30 Co 0.66 Ge 0.34 

3 45 Co 1.00 

4 1300 Si 1.00 0 2.00 




171 


Energy (MeV) 


0-6 07 0-8 0-9 1-0 M 



Figure 6.14: Rutherford backscattering spectra of Co/Ge bilayer thin film on quartz 
substrate for as-deposited condition and after room-temperature mixing with 1 MeV 
Kr ions at a dose of 1 x 10^® A'r/cm^. The corresponding simulated spectra are shown 

by solid lines. 




Co/Ge ( Ion-beam mixing, RT ) 


15.0 

4>(10'^ions /cm 2) 


Figure 6.15: Dose dependence of squaxed thickness of mixed layers of composition 
C 02 G& formed across the interface of Co/Ge bilayer system irradiated at room tem- 
perature with 1 MeV Ar ions (Q) and Kr ions (•) respectively. 
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Figure 6.16; The squared thickness of the reacted layer of composition CotGe formed 
across the interface of Co/Ge bilayer system after thermal annealing is shown as a 
funtion of annealing time. 







175 




Figure 6.18: The variation of miring rate for Co/Ge bilayer system as a function of 
reciprocal of substrate temperature for a dose of 1 x 10'' A'r/cm' (b) Arrhenius plot 
of the miring rate as a function of reciprocal of substrate temperature extracted torn 

Fig.e.lTa. 
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6.4 Discussion 


6.4.1 Critical temperature 


At first we will discuss the nature of the mixed region on the basis the Q cur vq presented 
in Figs 6.5, 6.13 and 6.18. It has been described in Chapter 2 that, in general, there 
exists two critical temperatures Tc and Teq with Tc < Teq such that metastable phases 
are expected between Tc and Tgq and equilibrium phajses are expected beyond Teq. de 
Reus et ai.[51] have correlated these temperatures with the energies required for creat- 
ing a hole due to the smaller and the larger element in the compound. The smaller hole 
formation energy is related with Tc while the larger one with Teq. These temperatures 
are predicted from the following expressions ( Eqs.2.1 and 2.2 from Chapter 2) 


Tc = 4.0 X AHvsmall (6.1) 

and 

Teq = 3.8 X AHviarge (6.2) 

where the temperatures is given in Kelvin and the hole formation enthalpies in 
kJ/mole. The hole formation enthalpies A.Hv for an alloy Ai—xBx are calculated from 
Miedema s formula described in Chapter 2 (Eq.2.3 and 2.4). For several systems the 
difference between Tc and Teq is not so pronounced and in those cases there is only one 
transition temperature[51]. 

Cheng et al.[20] have correlated Tc with average cohesive energy of mixing elements. 
They have plotted experimentally determined Tc values against AHcoh 3^4 obtained a 
straight Hne whose slope is 100 K/eV. Thus from Eq. 2.5 of Chapter 2 we have 

Tc^lOOAif^oA (6.3) 

For the present case of Cu/Ge system (see Fig 6.5) the value of Tc is 320 K. Theo- 
retically, the models of de Reus et al. (Eq 6.1) and Cheng et al. (Eq 6.3) give Tc values 
of 253 (± 75) K and 373 (± 100) K respectively. According to the formulation of de 
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Reus et al. the value of Teq is 330 (± 75) K. The comparison of calculated Tc and Teg 
values with the experimental results is reasonably well and is presented in Table 6.4. 
The values of cohesive energy (AHcoh) for the present systems are given in Table 6.5. 
Similarly, the calculated values for Ni/Ge and Co/Ge systems are compaxed with the 
experimental one in Table 6.4. It is to be noted that in all the three cases the equilib- 
rium phases have been observed around Tc. It has been found that several silicides form 
around Tc which could be due to the higher reactivity of these materials as compared 
to the metal/metal systems. For example, in Co/Si system the experimental [89] and 
theoretical[51] values of Tc and Teq are both equal to 300 K. In this case the formation 
of equilibrium phase Co 2 Si is reported to occur at room temperature under Xe ion 
bombardment [89]. In fact there are many more systems where an equilibrium phases 
are formed around Tc- Several such examples are shown in Table 6.4. Also note that 
the value of Tc depends on the mass of the irradiating ion[23, 29]. The value of critical 
temperatures increases with decreasing mass of the irradiating ion. 

6.4.2 Mixing rate: Prediction of existing models 

In this section the experimental values of mixing rates and mixing efficiencies will 
be compared with the values calculated for Cu/Ge, Ni/Ge and Co/Ge systems using 
empirical formulae proposed by different researchers (described in Chapter 2) . The 
experimental values of mixing rates obtained at room temperature which is just be- 
low the critical temperature (see Table 6.4) for these systems. Each phenomenological 
model assumes certain mechanisms viz. cascade mixing, thermal spike etc. and usually 
expresses in the form of an empirical formula for the evaluation of mixing rate. Thus 
comparison of the experimentally obtained values with the ones calculated using these 
models can reveal the underlying mixing mechanism in the systems under consider- 
ation. It should be borne in mind that these models also assume certain conditions 
viz. unlimited supply of atoms of each species, size and shape of the collision cascade, 
thermodynamic parameters, mobility of the specif involved etc. for their validity. 
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Table 6.4; Compari© 9 n of experimentally obtained critical temperatures with the one 
calculated from different empjerical models 


Bilayer 

system 

Critical temperature 

Tc in K 

Eq. temperature 

Tc, in K 

Observed 

phase 

around 

RT or Tc 

Ref. 

Exp. 

Calculated from 

Exp. 

Cal. from 


Cheng 

et al. 

model 

Eq 6.3 

de Reus 

et al. 

model 

Eq 6.1 

de Reus 

et al. 

model 

Eq 6.2 

Cu-Ge 

320 

368 

253 

320 

330 

Cu3Ge+ 

[this work] 

Ni-Ge 

326 

414 

338 

326 

455 

Ni2Ge+ 

[this work] 

Ni-Si 

250 

435 

390 

340 

399 

Ni2Si+ 

[ref 51] 

Go-Ge 

347 

412 

371 

347 

383 

G oiGeT 

[this work] 

Co-Si 

300 

451 

356 

300 

410 

C02Si+ 

[ref 51] 

Pd-Si 

* 

426 

373 


363 

Pd2Si+ 

7 ? 

Ti-Si 

380 

[ 

474 

392 

390 

426 

TiSi2 

7 ? 

Pt-Si 

' * 

524 

416 

320 

429 

' Pt2Si+ 

77 

Pd-Ge 

! * 

387 

371 

* 

383 

Pd2Ge^ 

77 

Nb-Si 

550 

610 

425 

570 

490 

1 

* 

77 


* Not available 

+ Phase formed around RT or Tc 
- Composition 
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Ballistic model 

According to this model[50] the mixing rate is expressed as (Eq. 2.8 of Chapter 2) 


6a/= ^ (6.4) 

$ CEd 

where C is the average atomic density, Ed is displacement energy or threshold energj’ 
for Frenkel pair production, ^ = [4:MiM2/{Mi + and A = 10 A is the minimum 
separation distance for the formation of stable Frenkel pair[14]. The values of average 
atomic density C, displacement energy Ed and f for Cu/Ge, Ni/Ge and Co/Ge systems 
are given in in Table 6.5. The deposition energy Fd for 1 MeV Kr and Ar ions in these 
three systems are given in Table 6.6. The mixing rates at room temperature for Kr 
and Ar irradiation calculated from this model for Cu/Ge, Ni/Ge and Co/Ge systems 
are compared in Table 6.7. It is seen from this table that the calculated ballistic 
mixing rates are one order of magnitude lower than that of the experimental ones. 
The calculated mixing efficiencies have been compared with the experimental ones in 
Table 6.8. Since all the three systems have similar masses and displacement energies, 
the collision process cannot account for more than 20% variation[14] in the obser\ed 
mixing rate, whereas actual values varied by a factor of four for the systems under 
investigation. 


Influence of chemical force and thermal spike 

The low values of mixing rates predicted by the ballistic model suggest that the col- 
lisional effects alone are inadequate to produce the observed mixing and that there 
axe driving forces of different nature which are also contributing to mixing[14]. These 
driving forces are termed as chemical forces since they depend on the heat of mixing, 
cohesive energy etc. We recall from Chapter 2 that the ion mixing rate of a Hf/Ni 
bilayer system is significantly larger than that for a Hf/Ti system. The collisional 
parameters for these two systems are similar but the heats of mixing for equal atomic 
concentrations are -62 and 0 kJ/gm.at. respectively showing the effect of chemical 
forces. The chemical forces act when thermal spikes are produced. On the basis of 



180 


Table 6.5: The values of average atomic number Zavei atomic density C, displacement 
energy Ed-, heat of formation A.Hfon cohesive energy AHcoh of Cu/Ge, Ni/Ge, Co/Ge 
and Ni/Si systems. 


Bilayer 

system 

Zave 


C in 

[7] 

Ed 

eV 

[67]. 

Equilibrium 

phase 

-AHfor 

(kJ/gm.at.) 

[63] 

AHcoh 

(eV/at.) 

[73] 

AHcoh 

Cu/Ge 

30.5 

1.00 

0.075 

17 

CuzGe 

4.05 

3.67 

0.012 

Ni/Ge 

30 

0.99 

0.072 

12 

Ni2Ge 

12.3 

4.30 

0.022 

Co/Ge 

29.5 

0.99 

0.071 

19 

Co^Ge 

5.00 

4.26 

0.012 

Ni/Si 

21 

0.94 

0.078 

20 

Ni2Si 

16 

4.5 

0.036 



181 


Table 6.6: The values of deposited energy density Fd, critical energy for global 
spike and the critical energy Ec for spike formation. 


Bilayer 

Irradiating 

Fd 

FE 

Ec 

Ref. 

system 

ion 

(keV/nm) 

(keV/nm) 

eV 


Cu/Ge 

Ar 

0.4 

0.95 

113 

[this work] 


Kr 

1.6 

2.69 

113 

7 ? 

Ni/Ge 

Ar 

0.5 

0.95 

109 



Kr 

1.6 

2.68 

109 

77 

Co/Ge 

Ar 

0.4 

0.94 

105 

77 


Kr 

1.4 

2.64 

105 

7 ? 

Ni/Si 

Au 

4.8 

5.53 

47 

i [ref 31] 
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Table 6.7: Comparison between experimental and calculated mixing rates. 

The theoretical mixed rates calculated from ballistic model (Eq 6.4) model 

Borgesen et al. (Eq 6.7), local spike model of Boise (Eq 6.8) (^ 
and global spike model (^ j^jo) (Eq 6.9). 


Bilayer 

system 

Fd 

(keV/nm) 

\exp 
^ 1 

(nm^) 

“ 5 “ \hal 

(nm*) 

Eq 6.4 

[Borgesen 
$ hoc 

(nm^) 

Eq 6.7 

X^ i Boise 
^ hoc 

(nm^) 
Eq 6.8 

f u 

{nm‘^) 
Eq 6.9 

Cu/Ge 

0.4 

5.4 

0.06 

0.11 

0.17 

0.42 

3 ? 

1.6 

20.6 

0.25 

0.44 

0.70 

0.68 

Ni/Ge 

0.5 

1 

5.6 

0.11 

0.13 

0.20 

0.55 

33 

1.6 

18.5 

0.35 

0.40 

0.64 

0.58 

Co/Ge 

0.4 

1.6 

0.06 

0.09 

0.13 

0.03 

33 

1.4 

7.7 

0.22 

0.32 

0.47 

0.48 
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Table 6.8: Comparison between experimental |“'p) and calculated mixing efficien- 
cies from ballistic (Eq 6.4) model lia/)) Borgesen et al local (Eq 6.7) spike model 
and Boise local spike (Eq 6.8) model Global spike (Eq 6.9) 

niodel(^ U). 


Bilayer 

system 

Jexp 
<PFn 1 

(nm^/keV) 

<SFx) 

(nm^/keV) 

i Borgesen 
l/oc 

{nm^lkeV) 

\Bolse 
iPp iloc 

{nm^jkeV) 

(nmykeV) 

Cu/Ge 

13.1 

0.16 

0.20 

0.39 

0.73 

Ni/Ge 

11.3 

0.22 

0.25 

0.3 

0.73 

Co/Ge 

4.0 

0.15 

023 

0.23 

0.20 
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fractal geometry, Cheng[14] has shown that the value should be larger than 20 
for thermal spike to occur. He’ has demonstrated that the space filling cascades which 
embodies the idea of thermal spikes can be initiated if the recoil energies fall below the 
critical value Ec (in eV) given by the following equation 

E, = ( 6 . 5 ) 

where Zave is the average atomic number of the target. The calculated values of Ec for 
Cu/ Ge, Ni/ Ge and Co/ Ge systems are given in Table 6.6. These values are much higher 
than the displacement energies Ed (Table 6.5) for the atoms indicating the occurrence 
of thermal spikes in all these three systems. The thermal spikes are generally divided 
into two categories(i) those formed by overlapping subcascades (global spikes) and (ii) 
those formed by non-overlapping subcascades (local spikes). 

As mentioned in Chapter 2, the density of spikes, and the conditions of their overlap 
as assumed by Johnson et al.[19], depend on the deposited energy density Fd and on 
the mass and energy of the incident ion. The threshold deposition energy[24] density 
also called as critical energy F§' above which the subcascades begin to overlap can be 
estimated from the formula (Chapter 3 Eq 2.15) given below. 



ZiZt 

Z9-23 ^ ^ 0.23 


[ 


m,- 


m,- d- Mt 


] 


( 6 . 6 ) 


Here, the constant k,n is of the order of 25 x IQ-^ ev-nm^[2i, 65]. The values of F^ 
calculated for metal/Ge systems under investigation are compared with the Fd values 


(for Kr and Ar ions) used in this work and shown in Table 6.6. The experimental 
values are less than the threshold values required for overlapping of local spikes. Thus, 
it IS expected that the mixing should occur under local spikes in the present cases. 

Borgesen et al.[22] and Bolse[24] have separately suggested on the basis of their 
senes of experimental results that the linear dependence of the mixing rate with Fd 
can be correlated with the formation of local spikes. Our experimental results also 
show a linear dependence of mixing rate with Fd (see Figs 6.5, 6.11, 6.17) suggesting 
the formation of local spikes in this systems. 
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6.4.3 , Mixing in non-overlapping subcascade 

By combining the effects of cohesive energy and heat of mixing under non-overlapping 
spherical spike condition and using regular solution approximation the mixing rate 

i 

from the model of Borgesen et al.[22] can be written as (chapter 2 Eq 2.13) 
y 2 L ' 71.5 p 

^ I Borgesen ^ 

where AHmix and AHcoh are the heat of mixing and cohesive energy respectively. 
Here Kb = 0.034 nm-eV^/^ and K 2 = 27.4. Later Bolse[24] have shown that these 
spherical local spikes overlap and form local cylindrical spike. The mixing rate then 
becomes (chapter 2 Eq 2.14) 


5 AM rnixl 

^ + 6^^^ ak: 


(6.T) 


£! KsZI^Fe 




1 




( 6 . 8 ) 


where Kb = 0.0163 nwP eV and A '2 = 41 are empirical constants obtained by fitting 
the mixing rate of metal/metal systems[76]. Here (in Eq 6.7 or Eq 6.8) first factor is 
due to the non-overlapping subcascade and the second one is the Darken term which 
indicates the energ}- gained by forming AB pairs after Treely moving’ A and B atoms 
generated in the collision cascade[104]. The calculated mixing rate and efficiencies for 
all the systems are compared with the experimental results in Table 6.7 and Table 6.8. 
The values of AH„ix AHcoh, C, and Fd are given in Tables 6.5 and 6.6. Thus by 
taking into account the heats of formation, cohesive energies and deposited energies 
we find that the mixing rates calculated from the local spike models are higher than 
those obtained using the ballistic approach. However, these values are still about an 
order of magnitude lower than the experimental mixing rates. 


6.4.4 Mixing under Overlapping subcascade 

Johnaon et al-llS] proposed this model for calculating mi.clng rates under regular so- 
lution approximation and is known as “global spike model". According to this model 
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mixing rate is given by (Chapter 2 Eq 2.12) 

^!i r, . 

* [* + ' 2 Aif„A J 

where A'l = 0.00175 nm and = 27 are empirical constants. Although the above 
formula is not applicable for present systems since the spikes are expected to be of non- 
overlapping nature', we have calculated the mixing rates for the sake of comparison. 
These values are given in Table 6.7 and are not very different as compared to the ones 
calculated from the local spike models. 



6.5 Need for a new model 

The inadequacy of all the above mentioned models indicates that there are some other 
factors which enhance mixing. There are many metal silicide[31] systems whose mixing 
rates are much higher than the ones calculated using above models. For example, 
recently, Desimoni and Traverse[31] have reported that the occurs in Pd/Si 

system and the mixing rate varies linearly as a function of Fp. They have found that 
the mixing rates calculated using the Johnson and Cheng model are 3 to 8 times lower 
than the measured mixing rates. All these results together with the results presented 
here strongly suggest that there are some other factors which contribute to the mixing 
process leading to compound formation. It may be noted that metal/Si and metal/Ge 
systems are strongly reactive. A formalism has been presented in the next section 
to calculate the mixing rates in such strongly reactive systems which is based on a 
mechanism in which the mobile species is driven to the growing interface where it gets 
consumed during compound formation along with simultaneous structural relaxation. 
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6.5.1 Model for compound formation under local spike mix- 
ing 

It is now well recognized that thermodynamics has to be taken into account in de- 
scribing ion-beam mixing and is normally manifested through heat of mixing. Johnson 
et al.[19] and Cheng[14] have shown that under regular solution approximation the 
diffusion coefficient can be expressed as 


D = Dq {1 — 2AHmixlkBT) (6.10) 

where Do is a constant, ks is Boltzmann constant, AHmix is the heat of mixing and 
T is the local temperature. The term within the bracket is responsible for enhancing 
the diffusion coefficient and is known as Darken term. 

However, in those cases where a chemical reaction also takes place at the inter- 
face of growing mixed layer, the reaction rate may strongly influence the concentration 
gradients and thus can become the dominant source of thermodynamically controlled 
modification of diffusion coeflficient. This effect of reaction kinetics on mass trans- 
port coefficients is well known in the analysis of compound forming thermal diffusion 
couples[105-107]. Desimoni and Traverse[31] extended this analogy for thermally syn- 
thesized compounds for explaining the ion-beam mixing results of silicides. 

This is best illustrated (see Fig 6.19) by considering a bilayer system A/B, where 
each layer is an infinite source of corresponding atoms; and for simplicity let only 
one atom, say A, be the moving species during compound formation. Assume that the 
solubility of A and B is very low but they are highly reactive resulting into a compound 
AaBb. On irradiation, initially the rate of the growth of the layer A^Bb is limited by the 
rate of reaction. However, at later times the supply of the mobile species to the reacting 
interface gets limited by its diffusion through the mixed layer(31]. If the reaction rate 
is high, then each atom arriving at the interface is consumed in the formation of the 
compound, and hence for the steady state growth, the flux (Ji) of mobile atoms due 
to diffusion must equal to the flux ( Ja) participating in the formation of the compound 
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Figure 6.19: Bilayer sample configuration of two elements A and B where compound 
AaBb grows across the interface under ion bombardment. 

at the interface. Assuming a small layer of thickness X, this equality of steady state 
flux Fs can be written as, 


F. = = kC, (6.11) 

where Ca is atomic density of layer A and Ci is the concentration of A atoms at 
the growing interface, k is the reaction rate and D is the diffusion coefficient. Note 
that the reaction rate coefficient determines ’’reaction kinetics" which in turn influences 
transport kinetics ' in a diffusion controlled process. The relationship between kinetics 
and thermodynamics in phase formation by ion-beam mixing has been discussed by 
Miotello and Kelly [26]. 

The time dependence of thickness X of the mixed layer can be expressed as steady 
state flux divided by the concentration of atom A in the compound, 

dt 


= $ 


II 

Na 

dX 


d^ 


( 6 . 12 ) 
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where is the conceiitratioa of A in the compound AaBb with atomic density 
C i.e., A'. 4 =[a/(a+b)]C. The second equality is obtained by replacing time by dose 
assuming flux $ to be constant. 

Using Eqs 6.11 and 6.12, elinxinating Ci and solving for X with the initial condition, 
X (t=0) ~ 0 one obtains, 


{NA^k)X^ +- {2 DNa^)X - 2DkCA^ = 0 
This gives two solutions for X of which the positive one is expressed as 


(6.13) 




I, ^ 


1/2 


1 


(6.14) 


For diffusion limited process and high reaction rate k, X becomes independent of k 
giving rise to 


X 


D 

k 



CA±]f_] 

NAiD^j 


1/21 


Ca ^ $ 

^ ■= 


(6.15) 


Note that though k does not explicitly appear in this equation, it is responsible for 
enhancement of D to some De// by a factor 2(7.4 /AU, whose origin can be attributed 
to the chemical driving force for formation of a compound[31]. In the approximation 
that there is no terminal solid solution, this is the only enhancement factor different 


from the Darken term in Eq 6.19. It is straight forward to extend Eq 6.15 to the case 


when both atomic species axe mobi le and is given by, 



[Na Nb. 


4 


(6.16) 


It is interesting to note that the tenhancement factor appears as a ratio of concentra- 
tion of mobile species to its concentration in the compound and is independent of the 
deposited energj’[31]. However, a nler irradiation conditions the diffusion coefficient is 
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separately dependent on deposited energy jFd- This dependence is known to vary as 
Fq in case of overlapping subcascades[14] and only as Fd in case of non-overlapping 
subcascades [22, 24]. It has already been mentioned that in case of non-overlapping sub- 
cascades/local spikes, under regular solution approximation, the mixing rate ^ 
for the later case is given by (Eq 6.8) 


$ 


\Bolse 

\loc 


= Kb 




l+Ki 


^Hcoh - 


(6.17) 


The first factor is due to the spike effect and the second one is due to the effect of 
chemical driving force. In case of compound formation, this would be further enhanced 
by the factor in the square bracket in Eq 6.16. However, in the absence of terminal 
solution, the Darken term can be simply replaced in the above equation by the com- 
pound formation enhancement factor in Eq 6.16 and the mixing rate ^ would be 
given by 


V2 

{com Tir 

IT I/AC - 


+ 


Cb 


2FdZI^ 


coh 


(6.18) 


[al{a-\-b) hl{a + h)\ 
where Kg is the proportionality constant. When atom A (B) is the only moving 
species during phase formation then the first (second) term within the bracket will be 
retained. 

The quantitative success of the above formalism in predicting the experimental 
mixing rates for metaJ/Ge systems presented here as well as the mixing rate of Ni/Si 
system reported in the iiterature[31] is shown in the following subsection. 


6.5.2 Comparison with experimental results 

To calculate the mixing rates, the values of compound formation factor and thermal 
spike factors and the proportionality constant A", are required. To determine the 
compound formation factor in Eq 6.18 for the present systems, it is essential to know 
whether only one or both the species are mobile during the formation and growth of 
the germanide phases. In order to investigate this we have studied the movement of 
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atomic species across the interface by marker experiments the results of which have 
been presented in the previous chapter. In Chapter 5 it has been shown that in case 
of Cu/Ge systems, both Cu and Ge axe mobile but Cu is the dominant mobile species 
whereas in case of Ni/Ge system Ni and Ge both are equally mobile. However, in case 
of Co/Ge system Ge is the only mobile species. 

The values of mixing rates calculated using Eq 6.18 for the metal/Ge systems are 
tabulated and compared with the experimental values in Table 6.9. The parameters 
used for these calculation are given in Tables 6.5 and 6.6. The calculated values are 
very close to the experimental mixing rates indicating excellent predictionabilty of the 
present model. The calculated mixing efficiencies for these systems are also quite close 
to the calculated values which are compared in Table 6.10. 

For the sake of completeness, we shall now compare the experimental mixing rate 
with the calculated ones from the modified relation of the ballistic mixing rates. Ac- 
cording to Desimoni and Traverse the ballistic mixing rate ^ vrritten as 

(Chapter 2 Eq 2.16) 




com 

bal — 


c. 


+ 


Cb 


[a/{a-h6) b/ia + b) 




(6.19) 


The mixing rates obtained from this modified expression of ballistic model are com- 
pared with the experimental one in Table 6.9. This shows that mi.xing rates calculated 
from the modified ballistic model are one order of magnitude less than the experimental 
values. Now we shall calculate mixing rates using the modified global spike model[311 
where mixing rate is proportional to Although we have seen a Hnear dependence 
of mixing rate on Fd we would like to compare the predictions of this model with our 
experimental values. The mixing rate is given by ( Eq 2.17 of Chapter 2) 


$ 


com_ E' 
glo — -''•1 


Ca 


■ + 


Cb 


2Fl 




( 6 . 20 ) 


[a/(a -f b) bf{a + b) 

The values calculated using the above relationship axe also given in Table 6.9. The 
mixing efficiencies calculated from these two models are compared with the expenmen- 
tal one in Table 6.10. It is seen from both the tables that the model presented in this 



192 


Table 6.9: Comparison between experimental and calculated mixing rates. 

The theoretical mixed rates |f°^) calculated from our model (Eq 6.18) agree well 
with the experimentally obtained mixing rates. Note the significant departures of the 
rates calculated from modified ballistic (Eq 6.19) model and global spike 

(f IX) (Eq 6.20). 


Bilayer 

system 

phase 

formed 

Moving 

species 

Fd 

(keV / nm) 

X'^ lexp 
0 1 

[nrn^) 

\com 
^ iloc 

{nm^) 

Eq 6.18 

* 

[corn 
^ I6a/ 

(nm^) 
Eq 6.19 

# 

XI Icom 
$ l^/o 

(nm'*) 
Eq 6.20 

# 

Cu/Ge 

CuzGe 

both 

0.4 

5.4 

5.36 

0.49 

0.25 


1 


1.6 

20.6 

21.46 

1.96 

4.0 

Ni/Ge 

A'isGe 

both 

0.5 

5.4 

4.14 

0.78 

0.25 



77 

1.6 

18.5 

13.25 

2.49 

2.51 

Co/Ge 

CoiG^ 

Gft 

0.4 

1.6 

1.7 

0.21 

0.09 


?? 

77 

1.5 

7.73 

6.5 

0.79 

1.27 

Ni/Si 

-Ai 

Ni 

4.8 

11.8 

8.92- 

2.04 

9.27 


* Model proposed in this work. 

# Model proposed in ref 31. 
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Table 6.10: Comparison between experimental |«p) and calculated mixing effi- 
ciencies. The theoretical mixing efficiency |«-) calculated from our model (Eq 
6.18) agree well with the experimentally obtained mixing efficiency. The mixing effi- 
ciencies from modified ballistic (Eq 6.19) model and global spike (|r) 

model (Eq 6.20) are also shown for comparison. 


Bilayer 

system 

JCL |«P 
^Fd I 

{nm^fkeV) 

\com 

4>Fd I/oc 

{nm^/keV) 

* 

|com 

<iFD 

(nm^/keV) 

# 

\com 
iFo <9io 

{nm^fkeV) 

# 

Cu/Ge 

13.1 

13.4 

1.23 

1.56 

Ni/Ge 

11.3 

9.3 

1.56 

1.01 

Co/Ge 

4.0 

4.35 

1.08 

1.09 

Ni/Si 

2.46 

1.86 

0.42 

1.93 


* Calculated from the proposed model. 

Calculated from the model presented in ref 31. 
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work provides the best agreement with the experimental results. 

The formalism developed above is applied to calculate the mixing rate for Ni/Si 
system for which the experimental data is available in the literature[31]. Desimoni 
and Traverse have done mixing experiments on Ni/Si system to produce nickel silicide 
{Ni 2 Si) phase and explained the experimentally obtained mixing rate by the modified 
global spike model for compound formation. Since, Ni/Si system consists of medium 
Zave elements, it is expected that the model proposed here is appropriate to explain 
the mixing rates in this case as well. Here, Ni is the only moving species during the 
nickel silicide phase {Ni 2 Si) formation hence only the first term in the compound 
formation factor (Eq 6.18) is to be considered. The calculated mixing rate from our 
model and experimental mixing rate at low temperature are in close agreement as 
shown in Table 6.9. The experimental and calculated mixing efficiencies are compared 
in Table 6.10. 

It is to be noted that we have used the same average value of the proportionality 
constant A® = 0.95 x 10~® nm^ — keV (Eq 6.18) in the present model for calculating 
the mixing rates for all the cases. A more general value of of 0.99 x 10“® nm^ — keV 
has been obtained from the plot (see Fig 6.20) of the experimental mixing rate versus 
the product of the compound formation factor and thermal spike factor (the right 
hand side of Eq 6.18). In general, these results suggest that this new model should 
be able to explain the mixing rates for many other reactive medium Zave systems for 
which compound formations have been reported in the literature. Unfortunately, the 
information regarding mixing rates for compound formation for these systems are not 
available. The knowledge of mixing rates will be helpful for further validation of the 
proposed model. 

Before concluding, it may be mentioned again that the growth of mixed layers at rel- 
atively low temperatures for highly reactive systems, such as silicides and germanides, 
must be viewed differently from those systems where a solid solution or amorphous 
phase is produced under ion irradiation prior to the formation of equilibrium com- 
pounds, In the later case long range atomic migration or overcoming of barrier to 
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Figure 6.20: Determination of average value of the proportionality constant A, used 
in Eq 6. IS form the experimental mixing rates. The best lit gives a value 0.99 x 10 

keV — nm^. 
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nucleation is necessary. However, for reactive systems where the diffusive flux of the 
reacting species is controlled by the reaction rate k, such barriers are either absent 
or negligible so that the system can undergo structural relaxation to produce stable 
equilibrium compound. The role of ion irradiation is to provide the required flux of 
reacting species for incorporation into the compound at the reacting interface. Hence 
models such as the one presented here and of Desimoni and Traverse[31] or Miotello 
and Kelly [26], would be of general applicability for systems with high reaction rate 
coefficient and low solid solubility. 


6.6 Conclusions 

The main conclusions of this chapter are, 

(1) For the bilayer structures of all the three Cu/Ge, Ni/Ge and Co/Ge systems, 
irradiation with Ar and Kr ions leads to the formation of CuzGe, Ni^Ge and CozGe 
phases respectively. The square of the thickness of ion mixed layers is linearly propor- 
tional to the irradiation dose of Kr and Ar ions and the mixing rates too exhibit linear 
dependence on deposited energy density Fd- The thermal annealing experiments also 
exhibit similar compound formation and show linear variation of squared thickness of 
the reacted layer with the annealing duration. 

(2) The room-temperature mixing rates for Ar ion irradiation in case of Cu/Ge, 
Ni/Ge and Co/Ge systems axe found to be 5.4 nm^, 5.6 nm^ and 1.6 nm* respectively. 
Similarly for Kr ion irradiations the mixing rates are 20.6 nm* 18.5 nm'* and 7.7 nm‘^ 
respectively. The mixing efficiencies are independent of the deposited energy density 
in each case and the experimental values for Cu/Ge, Ni/Ge and Co/Ge are obtained 
as 13.1 nm^/keV, 10 nm^/keV and 5 nm^JkeV respectively. 

(3) The temperature dependence studies of ion mixing in all the three systems 
show that the critical temperatures Tc, above which the mi.xing becomes temperature 
dependent, are just above the room temperature and are given by 320, 326 and 347 K 
for Cu/Ge, Ni/Ge and Co/Ge systems respectively. There is no distinction between 
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critical and equivalent temperature for crystallization. Compound formation in all the 
cases occurred around the critical temperature Tc- 

(4) The calculated mixing rate from aU existing models are much lower as compared 
to the experimental values. 

(5) A new model has been proposed for medium Zave elements where compound 
formation occurs during ion-beam mixing under non-overlapping subcascades. In this 
model the mixing rate is enhanced by the compound formation factor in highly re- 
active systems, and it inherits the linear dependence on deposited energy density of 
non-overlapping subcascades. The predictions by this model are very close to the exper- 
imental results presented here and the model is also able to explain observed mixing 
rate of Ni/Si system reported in the literature. Hence it is demonstrated that this 
model would be applicable to the important class of silicide and germanide systems. 



Chapter 7 


Summary and Conclusions 


The discovery of low resistivity Cu^Ge phase, produced by solid state reaction at 
elevated temperatures, reported at the beginning of this decade has triggered a great 
interest in metal germajiide phases since these materials have potential application as 
contact and interconnect in modern semiconductor technology. In the present thesis the 
syntheses of CuzGe and other metal germanide phases in Ni/Ge and Co/Ge systems 
have been investigated using ion-beam mixing (directed energy deposition method) 
which has distinct advantages of spatial selectivity and low processing temperatures as 
compared to the conventional method of thermal treatment at elevated temperatures. 

MeV Kr and Ar ions in the dose range of (2 - 20) x 10^® ions f cm} were used to 
induce mixing in the mutltilayer and bilayer films of the metal/Ge systems deposited on 
quartz substrate. Also the temperature dependence of ion-beam mixing in the range of 
100-450 K was investigated. Syntheses of the germanide phases by thermal treatment 
in the range of 400-700 K were also studied for comparison. To study the atomic 
movement caused by irradiation or thermal treatment, ion beam mixing experiments 
on bilayer films embedded with a thin marker layer at the interface were performed. 

The compositions and thicknesses of all samples before and after ion mixing or ther- 
mal treatment were determined using Rutherford backscattering spectroscopy (RBS). 
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The 2-MeV Van de GraafF accelerator at Indian Institute of Technology, Kanpur, was 
used for both irradiation and RBS measurements. X-ray diffraction techniques (conven- 
tional and small angle) were utilized for the identification of different phases produced 
after ion mixing or thermal annealing. The surface topography of the ion mixed and 
heat treated samples was examined after various doses and annealing treatments using 
scanning electron microscopy to detect the changes in surface features. 

The information on the synthesis of the first and subsequent phases was obtained 
from experiments performed with both multilayer and bilayer samples. On the other 
hand, the values of mixing rates, mixing efficiencies and critical temperatures for the 
onset of thermally activated processes were obtained from the results of the bilayer 
experiments. The nature of the phases formed in these metal germanide systems was 
examined on the basis of existing empirical rules. The ion-beam mixing mechanisms 
have been discussed and a new phenomenological model based on chemical kinetics, 
applicable to reactive systems ( silicides, germanides etc.) has been proposed to explain 
experimentally determined values of the mixing rates and efficiencies. 

The important conclusions of these studies are summarized as : 

1. We have found that MeV heavy ion irradiation technique is found to be a very 
efficient low temperature process for producing stable phases in Cu/Ge. Ni/Ge 
and Co/Ge systems. In case of Cu/Ge system, at lower doses, the equilibrium 
phase CuzGe forms first, followed by the formation of Cu^Ge phase at higher 
doses. In case of Ni/Ge system, only WijGe phase and in case of Co/Ge only 
Co 2 Ge phase form at all doses. The “first phases” to form by thermal annealing 
in these s}-stems are Cu^Ge^ Ni 2 Ge and Co 2 Ge respectively. These results show 
that for metal/Ge systems, as in the case of metal/Si systems, the first phase 

which is produced by ion-beam mixing is identical to the one produced by thermal 
annealing. 


The resistivity values of CusGe, NliGe and Co^Gt layers produced by ion-beam 
mixing are 20, 50 and 74 /ifi - cm. respectively. The resistivities of the phases 
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produced by thermal annealing are 8.6, 23 and 50 /iO — cm respectively. These 
resistivity values are lower as compared to metal silicides. 

3. We have shown that during the formation of Cu^Ge phase by ion-beam mixing 
at room temperature, both Cu and Ge are mobile but Cu is the dominant moving 
species across the interface while for Ni/Ge system both Ni and Ge are equally 
mobile for the formation Ni^Ge phase. However, for Co/Ge system only Ge 
atoms are mobile to produce the GoaGe phase. 

4. The first and subsequent phase formations in case of Cu/Ge and Co/Ge systems 
are correctly predicted by utilizing the concept of effective heat of formation and 
mobility of the participating atoms. 

5. The squared thickness of the mixed region (compound) in all the three cases, 
obtained after ion-beam mixing of bilayers at room temperature show linear de- 
pendence on irradiated dose for both Kr and Ar irradiations. Mixing rates at 
room temperature also exhibit linear dependence on deposited energ}’ density 
Fd- The growth of the mixed layer is like a diffusion controlled process. The 
thermal annealing experiments also exhibit similar compound formation and show 
linear dependence of the squared thickness of the reacted layers as a function of 

annealing time. 

6. The temperature dependence studies of ion-beam mixing in all the three systems 
show that the critical temperatures Tc, above which the mixing becomes temper- 
ature dependent, are just above room temperature. The values axe about 320 
K, 326 K and 347 K for Cu/Ge, Ni/Ge and Co/Ge systems respectively. There 
is no distinction between critical and equivalent temperature for crystalhzation. 
Compound formation in all the cases occur around the critical temperature 

7. The room-temperature mixing rates for Ar ion irradiation in case of Cu/Ge^, 
Ni/Ge and Co/Ge systems are found to be about 5.4 nm'*, 5.6 nm'* and 1.6 nm 
respectively. Similarly for Kr ion irradiations the room temperature mixing rates 
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are about 20.6 nm^, 18.5 nm^ and 7.7 respectively. The mixing efficiencies 
are independent of the deposited energy density in each case and the experimental 
values for Cu/Ge, Ni/Ge and Co/Ge are obtained as 1.3.1 nm^j'ktV , 10 nm^jkeV 
and 5 nm^jkeV respectively. These results and those for silicides reported in 
literature show that mixing rates and efficiencies are much higher than most 
of the metal/metal systems. The enhancement of mixing is is attributed to 
the highly reactive nature of metal/semiconductor systems as compared to the 
metal/metal systems 

8. The calculated mixing rates from all existing models are found to be much lower as 
compared to the experimental values. A new model has been proposed for systems 
composed of highly reactive elements where compound formation occurs during 
ion-beam mixing under non-overlapping subcascades. In this model the mixing 
rate is enhanced by the compound formation factor and it inherits the linear 
dependence of deposited energy density of non-overlapping subcascades. The 
predictions from the model are very close to the experimental results presented 
in this thesis as well as it predicts the mixing rate of Ni/Si system reported in the 
literature. It is expected that this model would be useful for predicting mixing 
parameters of other important silicide and germanide systems. 


Scope of future work 

The present work has demonstrated that metal germanide phases can be synthesized 
by ion-beam mixing at room temperature and above in Cu/Ge, Ni/Ge and Co/Ge 
systems. Several other raetal/Ge systems such as Cr/Ge, Fe/Ge, Nb/Ge, Mn/Ge, 
Pt/Ge etc. can be investigated on similar lines. These systems have stable equilibrium 
phases which could be synthesized by ion-beam mixing and thermal annealing. It is 
expected that the ion-beam synthesis of these phases and the studies of their various 
properties viz. resistivity, electro-migration, diffusivity, oxygen affinity, atomic mobility 
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etc. would make these materials useful for technological applications. 

With regard to the understanding of the mixing process, especially in reactive 
systems ( or in other words easily compound forming systems) such as metal/Ge, 
metal/Si etc., the validity of the phenomenological model presented here and the models 
reported in the literature so far, has to be checked further by performing man}^ more 
systematic experiments on mixing rates, mixing efficiencies, critical temperature etc. 
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